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HAUPTAUFSÄTZE 


Rückkehrkanten in ebener kompressibler Potentialströmung. 


Von @. @Guderley in Braunschweig. 







Mitteilung aus dem Institut für Gasdynamik der Luftfahrtforschungsanstalt Hermann Göring. 


Wie bekannt ergeben sich in analytisch gefundenen Überschallströmungen 
häufig Rückkehrkanten. Es wird an einigen Beispielen mit Hilfe der 
Prandtl-Busemannschen Charakteristikenmethode gezeigt, was bei 
solchen Strömungen physikalisch zu erwarten ist. Je nach dem Verhalten 
des Druckes an der Rückkehrkante ist diese entweder das bekannte Kriterium 
für das Auftreten von Verdichtungsstößen oder ein Anzeichen, daß sich der 
Strömungsrerlauf oberhalb ganz oder teilweise abändert, ohne daß dabei 
Verdichtungsstöße auftreten müßten. 




















In letzter Zeit haben die Arbeiten von Ringleb') das Augenmerk von neuem auf die 
schon früher von Tollmien?) gefundene Erscheinung gelenkt, daß bei der analytischen Be- 
handlung der Potentialströmung eines kompressiblen Mediums Rückkehrkanten auftreten 
können. Diese Ringlebschen Arbeiten ergänzend hat Tollmien®) den Beweis geführt, 
daß in jedem Falle die Strömung vor der Rückkehrkante „zusammenbricht“. Mit dieser 
Feststellung ist die Frage nach der physikalischen Bedeutung der Rückkehrkanten natür- 
lich noch nieht beantwortet®). Einige Überlegungen, die in dieser Hinsicht nützlich sein 
können, bringt die vorliegende Arbeit. 

Das wesentliche Hilfsmittel hierbei ist das Prandtl-Busemannsche Charakteristiken- 
verfahren). Es ist besonders anschaulich, weil bei ihm die Machschen Wellen nicht nur 
als Charakteristiken, d. h. als Linien, längs deren sich eine Störung möglicherweise fort- 
pflanzen könnte, sondern unmittelbar als Träger einer wohlbestimmten Zustandsänderung 


!) F. Ringleb: Exakte Lösungen der Differentialgleichungen einer adiabatischen Gasströmung. Z. angew. 
Math. Mech. Bd. 20 (1940), S. 185 bis 198. F. Ringleb: Über die Differentialgleichung einer adiabatischen Gas- 
strömung und den Strömungsstoß. Deutsche Mathematik, Bd. 5 (1940), S. 377 bis 384. 

2) W. Tollmien: Zum Übergang von Untersehall- in Überschallströmungen. Z. angew. Math. Mech. Bd. 17 
(1937), S. 117 bis 136. 

») W.Tollmien: Grenzlinien adiabatischer Potentialströmungen. Z. angew. Math. Mech. Bd. 21 (1941), S. 140 bis 152. 

°) Hierauf weist schon Tollmien hin, a. a. O. (Fußnote 3), S. 152. 


5) A. Busemaun: Beitrag Gasdynamik im Handbuch der Experimentalphysik (Wien-Harms), Leipzig 1931, 
Bd. 4, 1. Teil $ 26. 
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aufgefaßt werden®). Zur Orientierung in unseren Bildern diene die Bemerkung, daß Ver- 
diehtungswellen durch ausgezogene, Verdünnungswellen durch gestrichelte Linien gekenn- 
zeichnet sind. Der „Sprung“, d.h. die Winkeländerung, die der Geschwindigkeitsvektor beim 
Überschreiten einer Machschen Welle erfährt, ist durchweg 2°. 


1. Rückkehrkanten mit Expansion. 


Nach ihren physikalischen Eigenschaften sind zwei Typen von Rückkehrkanten zu 
unterscheiden, nämlich Rückkehrkanten, bei denen eine Kompression und Rückkehrkanten, 
bei denen eine Expansion stattfindet. Wir wenden uns zunächst zu einer Strömung, in der 
eine Rückkehrkante mit Expan- 
sion auftritt (Bild 1), sie ent- 
spricht etwa dem Beispiel, das 
in Bild 13 der ersten Tollmien- 
schen Arbeit?) dargestellt ist. 
Die Anströmung geschieht von 
rechts. Das analytisch berech- 
nete Strömungsfeld erstreckt 
sich bis zur Linie AB. Wie 
bekannt, ist eine Rückkehrkante 
die Hüllkurve Machscher Wel- 
len. Bei Rückkehrkanten mit 
Expansion handelt es sich, wie 
unser Beispiel zeigt, um Ver- 
dünnungswellen, deren jede auf 
der vorhergehenden ihren An- 
fang nimmt. Eine solche Strö- 
mung ist durch keine am Rande 

Bild 1. Strömung, in der eine Rückkehrkante mit yorsunchmenden Maßnahmen 
Expansion auftritt, Fortsetzung dieser Strömung realisierbar; die stärkste Kon- 
durch ein Schaufelgitter. zentration von Verdünnungs- 
wellen, die man einer Strömung 
vom Rande her aufzwingen kann, ergibt sich nämlich, wenn man dort eine Ecke anbringt. 
Man erhält dann einen Fächer von Verdünnungswellen, die in einem endlichen Abstand vom 
Rande jedenfalls auch einen endlichen Abstand voneinander haben und deshalb im Innern 
der Strömung nicht mehr auseinander hervorgehen können. Die Strömung mit Rückkehr- 
kanten besitzt also eine stärkere Konzentration von Verdünnungswellen, als sie durch irgend 
eine Randbedingung zu erreichen wäre. Erst durch Maßnahmen innerhalb des Strömungs- 
feldes, z. B. indem man längs der Rückkehrkante Ablenkschaufeln anbringt, ist das analy- 
tisch gefundene Strömungsbild zu verwirklichen. In Bild 1 ist unterhalb der Linie AB das 
hierzu notwendige System von Ablenkschaufeln nebst dem zugehörigen Strömungsverlauf 
angegeben. Da das Strömungsbild mit endlichem Sprung aufgezeichnet worden ist, ergibt 
sich ein System von endlicher Schaufelzahl. Geht der Sprung gegen Null, so geht die 
Schaufelzahl nach Unendlich. Um zu zeigen, daß sich die Strömung mit Sicherheit fort- 
führen läßt, sind die Wände so geformt, daß schließlich ein Parallelstrahl entsteht. Die 
Wirkung der Umlenkschaufeln besteht darin, daß sie auf das Gebiet oberhalb der Rück- 
kehrkante Einfluß nehmen und dort die von der analytischen Lösung geforderte Konzentration 
von Verdünnungswellen erzeugen. 

Aufschlußreicher als diese Überlegungen ist natürlich die Lösung der bereits von 
Tollmien formulierten Aufgabe, „die tatsächliche Strömung innerhalb der Berandung her- 
zustellen, die bis zum Ort der Grenzlinie mit dem Rand der Grenzlinienströmung überein- 
stimmt“. Sie läßt sich mit Hilfe der Charakteristikenmethode angeben, wenn der Abstand 
der berandenden Stromlinien einen gewissen, im allgemeinen nicht verschwindenden Betrag 
nicht überschreitet. Es sei das Strömungsfeld vom Unterschall her bis zur kritischen Ge- 
schwindigkeit und außerdem die gesamte Berandung wie bei der analytischen Lösung ge- 
geben. Als Berandung ist das frühere linke Ufer in Bild 1 sowie die Stromlinie CD gewählt 
(Bild 2). 

Wendet man nun die Charakteristikenmethode an, so zeigt sich, daß die Strömung über- 
haupt nur bis zu der vom Endpunkt A der Berandung ausgehenden, stromabwärts verlaufen- 


\AALLLLLDLAD) \ 














6) Genauer gesagt: Wegen der besonderen Eigenschaften der ebenen Potentialströmung läßt sich den Machschen 
Wellen eine solehe Bedeutung erteilen, daß beim Überschreiten einer Welle eine Zustandsänderung eintritt, die nur 
von dem örtlichen Zustand und dem für Gebiete des Strömungsfeldes einheitlich zu wählenden „Sprung“, nicht aber 
von der Länge des durchlaufenen Linienelements abhängt. 








gt. 
m 
n 
Ir- 
nd 
3S- 
ly- 
as 
uf 
bt 
lie 
rt- 
)ie 
k- 
on 


on 
r- 
n- 
nd 
ag 
re- 
se 


ılt 


T- 
n- 
ıen 


ur 
yer 














Z. angew. Math. Mech. 


Bd.22 Nr.3 Juni 1942 Guderley, Rückkehrkanten in ebener kompressibler Potentialströmung _ 123 





den Machschen Welle A E bestimmt 
ist. Bis dahin stimmt die analytisch 
berechnete und die mit der Charakte- 
ristikenmethode gefundene Lösung 
überein. Uber das weitere, bisher 
noch unbestimmte Gebiet verfügt die 
analytische Lösung durch die Fest- 
setzung, daß dort die Linien gleichen 
Zustands, wie bisher, Kreise sein 
Bild 2. Strömung, die zwei Stromliniendes Sollen, im allgemeinen also durch 


ey z analytisch gefundenen Teils in Bild 1 als analvtische Fortsetzung In Wirk- 
S Berandung besitzt, aber keine Rückkehr ü BEE R Se. 
nun £ ee lichkeit ist dieses Gebiet jedoch durch 


die Randbedingungen bestimmt, die 

man unterhalb von A vorschreibt und 
mit deren Hilfe sich, wie oben gezeigt wurde, die analytische Fortsetzung nicht erreichen 
läßt. In Bild 2 ist in A zunächst eine Ecke angebracht und der weitere Verlauf der Ufer 
so geformt worden, daß sich schließlich ein Parallelstrahl ergibt. 

Das Charakteristikenverfahren ergibt natürlich auch für jede glatte Fortsetzung der 
Kontur eine mögliche Fortsetzung der Strömung; die Ecke ist angebracht worden, um einen 
möglichst großen Abstand der berandenden Stromlinien zu erhalten (vgl. den folgenden Absatz). 

Diese Konstruktion versagt, wenn die beiden Ränder so weit voneinander entfernt 
sind, daß sich dabei Unterschallgebiete ergeben. Wie in diesem Fall die Strömung aussieht, 
läßt sich nicht ohne weiteres angeben, insbesondere ist zunächst nicht entschieden, ob sich 
die Änderung der analytisch gefundenen Lösung auch auf die Unterschallzuströmung er- 
streckt. Es ist durchaus denkbar, daß sich nur innerhalb der Überschallströmung isolierte 
Unterschallgebiete bilden, während die Unterschallzuströmung erhalten bleibt. Dies ist 
sicherlich nicht mehr möglich, wenn die Machsche Welle, die den durch die Rand- 
bedingungen bestimmten Bereich begrenzt (in Bild 2 die Linie AE) nicht auf der Kontur, 
sondern schon vorher auf der Linie der kritischen Geschwindigkeit endet. In diesem Fall 
nimmt die Fortsetzung, die man der Berandung unterhalb von A erteilt, Einfluß auf die ge- 
samte Unterschallströmung, die dadurch in ihrem ganzen Verlauf geändert wird. Auf jeden 
Fall ist das Versagen des Charakteristikenverfahrens nur eine Grenze der Berechnungs- 
methode, die keinesfalls bedeutet, daß bei ihrem Überschreiten physikalische Überraschungen 
zu erwarten wären. 


2. Rückkehrkanten mit Kompression. 

Als Beispiel für Rückkehrkanten mit Kompression benutzen wir die zu Bild 1 gehörige 
„Rückkehrströmung* (Bild 3). Sie entspricht dem Bild 15 der ersten Tollmienschen Arbeit. 
Die Strömung möge aus dem Unendlichen (von oben her) kommen und auf die Rückkehr- 
kante zu führen. Die Rückkehrkante entsteht hier dadurch, daß Verdichtungswellen zu- 
sammenlaufen. Wollte man die Strömung bis zur Rückkehrkante verwirklichen, so müßte 
man ebenfalls Leitschaufeln benutzen, die hier die Aufgabe hätten, das Zusammenlaufen der 
Verdichtungswellen zu verhindern. Dieses Leitschaufelsystem und die sich darin einstellende 
Strömung ist das gleiche wie in Bild 1; in Bild 3 sind die Anfänge angedeutet. Wenn diese 
Schaufeln fehlen, laufen die Verdichtungswellen zu einem Verdichtungsstoß zusammen, den 
man, solange er schwach ist, ohne grundsätzlichen Fehler als Verdichtungswelle mit größerem 


Bll3. Strömung, bei der eine Rückkehrkante mit Kom- Bild 4. Strömung mit zwei Stromlinien des Bildes 3 als Be 
pression auftritt (Rückkehrströmung zu Bild 1). randung. An Stelle der Rieckkehrkante mit Kompression 
ist ein Verdiehtungsstoß getreten. 
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Sprung behandeln kann’). Dies ist in Bild 4 geschehen; der Verdichtungsstoß liegt etwas 
oberhalb der früheren Rückkehrkante. Auch hier gibt das Auftreten von Unterschallgebieten 
eine Grenze für die Durchführbarkeit des Verfahrens und beschränkt damit die Entfernung, 
die man den beiden Ufern geben darf. 


3. Anwendung auf das Ringlebsche Beispiel. 

Bei den Tollmie'nschen Beispielen wird die Rückkehrkante von jeder Stromlinie er- 
reicht‘). Einen anderen Strömungstyp, bei dem nur ein Teil der Stromlinien die Rückkehrkante 
trifft, zeigt das von Ringleb angegebene Beispiel, das im folgenden untersucht werden soll). 

Die Rückkehrkante weist hier Spitzen auf. Zwischen den beiden Ästen der Rückkehr- 
kante, die von einer solchen Spitze ausgehen, liegt ein Gebiet, das dreifach mit Zustands- 
größen überdeckt ist. Man kann es sich anschaulich vorstellen, indem man die Strömungs- 
ebene mit drei Blättern überdeckt denkt, deren jedes einfach mit Zustandsgrößen belegt ist. 
Längs der Rückkehrkanten sind diese Blätter zusammengeheftet. Es ist also so, als ob die 
Strömungsebene zwischen den Rückkehrkanten eine Falte aufweist. Diese Faltenbildung 
(Bild 5a) ist für diese Strömungen mit Rückkehrkanten typisch; betrachten wir nun das 
Ringlebsche Beispiel im einzelnen. 

Die Strömung möge von rechts unten aus dem Unendlichen kommen und nach rechts 
oben ins Unendliche führen. Zunächst soll die untere Hälfte der Strömungsebene ins Auge 
gefaßt werden (Bild 6). Von den beiden Ästen der Rückkehrkante ist nur derjenige einge- 
zeichnet, den die aus dem Unendlichen kommenden Stromlinien zuerst treffen (Kurve AB, 
die Spitze liegt in A). Diese Rückkehrkante AB ist eine Rückkehrkante mit Expansion. 
Zwischen zwei nicht allzuweit voneinander entfernten Stromlinien läßt sich, wie früher, eine 
Strömung angeben, die als Ränder diese ursprünglichen Stromlinien besitzt und keine Rück- 
kehrkante mehr aufweist. Dies ist in Bild 7 geschehen. Als Berandung sind dabei die 
strichpunktierten Stromlinien CD und EF in Bild 6 gewählt worden. Durch die Berandung 
zusammen mit den Anfangsbedingungen ist die Strömung nur bis zur Machschen Welle D@ 
bestimmt. Der weitere Verlauf hängt von der Fortsetzung ab, die man dem rechten Ufer 
über D hinaus erteilt. Sie wurde so gewählt, daß die Symmetrielinie 71 der ursprünglichen 
Strömung überall senkrecht durchströmt wird. Den weiteren Strömungsverlauf könnte man 
durch Spiegelung an der Linie IH erhalten. Der Abstand der Stromlinien, für den sich 
mögliche Strömungen mit Hilfe des Charakteristikenverfahrens angeben lassen, ist auch hier 
beschränkt. Ob sich beim Überschreiten dieses Abstands sofort die Unterschallströmung 
ändert, kann nicht entschieden werden. Sicher tritt dies ein, wenn an der äußeren Be- 
randung nirgends mehr Überschall auftritt. In diesem Fall ist die analytisch gefundene 
Strömung nur dann zu realisieren, wenn man das innere Ufer so weit in die Strömung 
hineinlegt, daß die Rückkehrkante ausgeschlossen bleibt. Ganz anders liegen die Verhält- 
nisse bei einer Rückkehrkante mit Kompression, wie sie sich in der oberen Hälfte der 
Strömungsebene ergibt (Bild 8). Hier kann man immer, unabhängig von der Lage der 
äußeren Berandung, den Verdichtungsstoß einzeichnen, der durch das Zusammenlaufen der 
Verdichtungswellen entsteht. Ist dieser genügend schwach, so entspricht die stromabwärts 
davon befindliche Strömung wiederum der analytischen Lösung; die Falte in Bild 5a wird 
einfach durch einen Sprung ersetzt (Bild 5b). Es ist hiernach zu vermuten, daß sich auch 
Beispiele finden lassen, wo der Verdichtungsstoß ins Unterschallgebiet führt. 

Einen Mangel besitzt dieser Lösungstyp allerdings noch, an der Stelle, wo der Ver- 
dichtungsstoß auf die Kontur trifft, weist diese eine Ecke auf. Das schränkt natürlich die 
Bedeutung solcher Strömungen ein. 

Bei diesen letzten Überlegungen ist angenommen worden, daß die Geschwindigkeits- 
verteilung, wie sie auf der Mittellinie 7 I herrscht, physikalisch realisierbar ist. Von der 
Zuströmseite stehen dem Schwierigkeiten im Weg, weil dort die innere Stromlinie so weit in 
die Strömung hineingelegt werden muß, daß die Rückkehrkante ausgeschlossen wird. Damit 
schneidet man jedoch wegen der Symmetrie der Strömung in bezug auf HI auch auf der 
Seite der Abströmung die Rückkehrkante und damit den Verdichtungsstoß aus. Um eine 
Strömung zu finden, die diese Schwierigkeit vermeidet, müßte man sich durch Kombination 
mehrerer der von Tschapligin stammenden Partikularlösungen der Differentialgleichung 
des Hodographen eine unsymmetrische Strömung mit Rückkehrkanten verschaffen. Dort 
läßt sich dann wahrscheinlich eine Berandung so wählen, daß die Rückkehrkante mit Expan- 
sion ausgeschlossen bleibt, während eine Rückkehrkante mit Kompression, die dann als Ver- 
dichtungsstoß zu deuten wäre, dem Strömungsbereich angehört. 

”,Busemanna.a. O0. $ 27. 

8) Auf einen Irrtum, der mir an dieser Stelle unterlaufen war, machte mich Herr Tollmien freundlicherweise 


aufmerksam. 
%) Z. angew. Math. Mech. a.a.0O. Bild 4 oder Deutsche Mathematik a.a.0O, 
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Bild 5a (links). Gestali der Strömungs 
ebene beim Auftreten von Rückkehr 
kanten im Ringlebschen Beispiel. 


Bild 5b (rechts). 
Beim Auftreten 
eines Verdich- 
tungsstoßes wird 
die Falte im 
Bild 5a durch 
einen Sprung er- 
setzt, 


u. 


Bild 6 (oben). Rückkehrkante mit Expansion in dem Ringlebschen Beispiel. 


Bild 7 (links). Strömung, die die Stromlinien CD und EF in Bild 6 als Berandung 
besitzt, jedoch keine Rückkehrkante aufweist, 


Bild 8 (rechts). Die 
Rückkehrkanten mit 
Kompression in dem 
Ringlebschen Bei- 
spiel sind durch einen 
Verdichtungsstoß er- 
setzt worden. 
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Zusammenfassung. 


Hinsichtlich ihrer physikalischen Bedeutung bestehen zwischen Rückkehrkanten mit 
Expansion und solchen mit Kompression wesentliche Unterschiede '). 

Strömungen, in denen Rückkehrkanten mit Expansion auftreten, sind durch die Ge- 
staltung ihrer Ufer allein in ihrem gesamten Verlaufe nicht zu verwirklichen. Erst durch 


Maßnahmen innerhalb des Strömungsfeldes — man müßte unterhalb der Rückkehrkante, aber 
jedenfalls nicht in dem analytisch gefundenen Bereich Leitschaufeln anbringen — kann’ die 


analytische Lösung erzwungen werden. Wo diese Maßnahmen fehlen, tritt eine Abänderung 
des Strömungsverlaufs schon oberhalb der Rückkehrkante ein. Wenn sich diese Änderung 
auf das Überschallgebiet beschränkt, läßt sie sich mit Hilfe der Charakteristikenmethode 
angeben. Die Strömungen, die man dabei erhält, weisen keinerlei physikalische Besonder- 
heiten auf. Die meisten Vermutungen über die Bedeutung von Rückkehrkanten, die das 
bisherige Schrifttum enthält, sind danach für Rückkehrkanten mit Expansion nicht zutreffend. 
Auch in den Fällen, wo die Änderung der analytisch berechneten Strömung nicht mehr an- 
gebbar ist, ist das Auftreten physikalischer Besonderheiten unwahrscheinlich. Rückkehr- 
kanten mit Kompression sind ein Zeichen für das Auftreten von Verdichtungsstößen. Bei 
geeigneten Beispielen läßt sich der Verdichtungsstoß näherungsweise sofort angeben. 350 














10) Die bisherigen Arbeiten machen wegen ihrer vorwiegend mathematischen Zielsetzung diesen physikalischen 
Unterschied nicht. 
























Eingrenzung von Verschiebungen mit Hilfe der Minimalsätze. 


Von Constantin Weber in Dresden. 








Zur Einführung in den nachfolgenden Aufsatz wird hier zunächst die dort 
allgemein entwickelte Methode an einem einfachen Beispiel (Eingrenzung 
der Verschiebung einer belasteten quadratischen Membrane) auseinander- 
gesetzt. 


Zu meinem Aufsatze, der diesem folgt, möchte ich hier ein Beispiel bringen. Um ein 
recht einfaches Problem zu erhalten, ist eine quadratische gleichmäßig belastete Membrane, 
Bild 1, gewählt, deren Verschiebung in der Mitte eingegrenzt werden 

Pe Be soll. Die sehr einfachen Näherungsansätze geben natürlich nur 


\ 


eine rohe Eingrenzung, zeigen aber das Wesentliche der Methode. 





I. Die Minimalsätze. 
Die elastischen Eigenschaften der Membrane seien derart, 
daß bei einer Verschiebung » (x, y) die Zerrungsenergie je Flächen- 
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Das Minimum der potentiellen Energie gibt: 
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J(w, p) = \ 15 
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(A357.1) 
Bild 1. Belastete Membrane. Setzen wir verschiedene Verschiebungen (x, y) ein, die den Rand- 
bedingungen genügen, so machen die wirklichen Verschiebungen 


den Ausdruck J(w,p) zum Minimum, das gleich J„= —\ w„pdf wird. Mit Hilfe der 
r . . . . . N 
Variationsrechnung erhält man die Eulersche Gleichung 

0° Wi 0° m 


2 


Jet ap rpm. ern. 
mit We Raı 

Das Variationsproblem läßt sich nach Friedrichs!) in ein zweites Variationsproblem 
transformieren; hierbei wird der Minimalwert des ersten Variationsproblems zum Maximal- 
wert des zweiten (dieses zweite Variationsproblem entspricht nach Umkehr der Vorzeichen 
dem Castiglianoschen Minimalprobleme). 


er 0. 


1) K. Friedrichs: „Ein Verfahren der Variationsrechnung .., .“*. Nachr. Ges. Wiss. Göttingen 1929. 
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Für die Membrane lautet das zweite Variationsproblem wie folgt: Zwei Funktionen o (x, %) 
und r(x,9) mögen die „Gleichgewichts“-Gleichung 








RE N 
tar» ZERERKTTEEETEIE 
befriedigen. Wir setzen: 
4 j 
K=\>-(0°+1)df=Min Be a ee 





m) 





Bei Wahl aller möglichen Funktionen o, r, die der Gl. (4) genügen, gibt die Lösung des 







ae = j “ . e Om Om 
Variationsproblemes (5) die wahren Funktionen o,, 7„; diese werden o,= N *, wa 7 
x 
und der Minimalwert des Variationsproblemes ist 
1 l[odwo\ , (dmwo\ 1{ 
Ku=\z|( “) d =Z|\pmwud N A EN a > A 
u | 2 hi) x + d Yy f > ) ji u f 
2) 0) 








II. Eingrenzung des Integralwertes der Variationsprobleme. 


Soll een w„df eingegrenzt werden, so nehmen wir erst für w einen 

h 
Näherungsansatz mit wx,nı=0, dann für o und r einen Näherungsansatz, der die „Gleich- 
Or 
%Y 






00 en 
gewichtsbedingung“ + +p=0 erfüllt. Mit Hilfe des ersten Ansatzes berechnen wir 





en en NT 





mit Hilfe des zweiten Ansatzes 






und hieraus 








BETLEMER . 22: ec 
m) 





Für die quadratische Membrane, die durch & +1, y + 1 begrenzt und durch p=1 be- 
lastet ist, nehmen wir z.B. 





w=c(1—x’)(1— y’), 
RED) mw. 


=—2cx(l —), 






=— 2cy(l— ) 






2 9y 
und finden 
+1 +1 
ar 
J= \ \ 3 eftkı1—yP’+4y(l — ar] —cll — a?) (1 7) dedy 

u 
18, 16 
ir A 










2 .) e" . > 
Für das optimale ZT: erhalten wir das kleinste J= 15. 
> fe 
“ : x Y \ ; SEN i 
Als Ansatz für o und r nehmen wir o=— 5, =—,, womit die Gleichgewichts- 






bedingung (4) erfüllt ist. 





Man erhält 





ı +1 


K= | Ss +mMdrdy=z. 


u 





Hieraus folgt die Eingrenzung 
’ +1 +1 
li 6 
- \ \ n„dxedy=s 18° 
a 
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18 





IN 





Diese Eingrenzung stellt nichts Neues dar und ist nur zum Verständnis des Folgenden voran- 
gestellt. 
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III. Eingrenzung der Verschiebung der Membranenmitte. 

Die gegebene Belastung sei p,, die wirkliche Verschiebung w,„:; eingegrenzt werden 
soll die Verschiebung der Membranenmitte. Wir nehmen eine zweite Belastung p, mit der 
Verschiebung w,,,, wobei die Gesamtlast \p,df=1 sei. Wir zeigen erst, wie der Ausdruck 
\p,w,wdf=\p,w,wdf eingegrenzt wird. An Stelle der Belastung p, nehmen wir etwas 
allgemeiner die Belastung Pp,. Für die Belastung p, wählen wir Näherungsansätze w, 
bzw. o,,t,, für die Belastung Pp, die Näherungsansätze Pw,, bzw. Po,, Pr,. 

Wir stellen die Ungleichung (7) erst für die Belastung p, + Pp,, dann für p, getrennt 
und für Pp, getrennt auf, ebenso die Ungleichung (8): 

1 Bu 
>\pı w,wdf P\p, wwdf—-5 r . (10), 


J(w,+Pw,,p, +Pp)Z 


i Lf j 
J(w,,P,) | 5 \p ww df EI EEE UNTER Pr re 1 


J(Pw,,Pp.,) > =P?:\p, (0, 


Ir 
K(o, +Po,,r,+Pr,) a \»: w,wdf+ P\p, %.udf +- 


K (o,,t,) \pwudf 


K(Po,, Pr,) 
Durch Addition der Ungleichungen (10), (14), (15) bzw. (11), (12), (13) erhält man 
P\p, w„AdAfzJ/(w,+Pw,,p +Pp)+Klo,tm)+K(Po,,Pr,). . . . (16), 
P\p,w,„d/= K(o, +Po,r,+PrT)+J(w,p)+J(Pw,,Pp.) - -» -» - MM. 


Durch jede der Ungleichung (16) und (17) ist \p,w,„df nach zwei Seiten eingegrenzt, da 
P sowohl positiv als auch negativ gewählt werden kann. 


Wählt man p, so, daß p, nur auf einer kleinen Fläche um den Mittelpunkt der Mem- 
brane von Null verschieden ist, so geht der eingegrenzte Wert über in 


\?.w,wdf>w,um \ 2 df=w,wm- 

(f) (/) 
%, wm Ist die wirkliche Verschiebung (Index w) durch die Belastung p, (Index 1) im Mem- 
branenmitte (Index m). 


Bei ungeeigneten Ansätzen von w, und o,,r, erhält man als Grenzen — und +». 
Um diesen Fall auszuschließen müssen die Näherungsansätze für w, bzw. o,,r, im Angriffs- 
punkt der Kraft 1, zu der die Belastung p, zusammenschrumpft, die richtigen Singularitäten 
aufweisen. 


Für unser Beispiel mit p, =1 machen wir folgende Ansätze 
für die Belastung p;: 
> 1 ‚2 1 2 
= 1-2), 


für die Kraft P, also für die Belastung Pp,: 
P 
Pe, = zo [In »® — In(1+x°”) —In(1-+ 9?) + In 2] 


> > y 
Px EEE: Py 
? Zar 


Po,= 


nr’ 








10), 
11). 
12), 


13), 


16). 
17). 


da 


'11- 


'Mm- 


ffs- 


en 
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Die Ansätze für die Kraft P haben im Punkte z=y=0 die gleichen Singularitäten wie bei 
einer Kreismembrane mit der Kraft P. 

Bei der Auswertung der Integrale sind die Anteile, die den Wert © geben, in Polar- 
koordinaten mit der unteren Grenze r =e ermittelt. 

Wir stellen uns hierzu vor, daß die Kraft P nicht in einem Punkte, sondern gleich- 
mäßig verteilt auf einer Kreislinie vom Halbmesser & angreift. Im Ansatze für Pw, ist dann 







“ RIES Z E . a - Lu Y 
fürr>e Inr’, fürr<e Ine’ zu setzen, ebenso im Ansatze für Po,, Pr, fürr >e m und 
3 3 2 : 





für r<e O und ®. 
Die betreffenden Anteile heben sich bei der weiteren Rechnung fort. 
Man erhält 








J(w,+Pw,,p, + Pp.): 
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*; E [In »?® — In(1 +.) —In(1!+y?) + In 2]) dedy >: MA 
> 2 2 / 1 3 
. ie! +5, [me+ 5 in?) 
P (2 rg 
5: 5 > /2 71 | P:? P: P? /z 1 
Tu (St) t7zme + \Ineospdp+ ze (5 ‚) 
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f) 
Die Werte J(w,,p,), J(Pw,.Pp,), K(o,,r,), K(Po,. Pr,) sind ebenfalls durch Gl. (18) 
und (19) durch die Glieder ohne P bzw. mit P? gegeben. 

Die Ungleichungen (16) und (17) geben 


3 P e P:? 
wm =jg 7 0,9285 + 7 
P? 


3 
+7 FREE 











— Pw er 












> Pr - 
Pn em 





- P und lassen den Querstrich fort, so erhalten wir 






Ersetzen wir in Ungleichung (20) P durch 


>. AP PER 10 27 22 
umSjgtz 85 +75 0,1427 . a 






Pn 









Ungleichungen (21) und (22) unterscheiden sich nur durch das mittlere Glied der rechten 
Seite. Allgemein lauten sie: 


Pwumzsa+tPb+P:ec. 






Wir teilen die Ungleichung durch |P| 


Re a P 
. ä Ip! 
|Pj"em> P| pjPt Pie. 






Für die rechte Seite wollen wir P,,: finden, für das sie den kleinsten Wert annimmt. Das 
> 


em Mi. i 
Glied ip b ist hierbei unveränderlich. Man erhält 










a 


pP? EEE Pop == % J ale. 
op 
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Für Ppr=+tyale folgt wum S2yac+b, 
für Pp=—YVale folgt wuom 2—2yacHtb. 


Jede der Ungleichungen (21) und (22) gibt Eingrenzungen nach beiden Seiten; aus (21) folgt 


1 a Er di u 
u 0,1778) Mom =—(1+0178), 


aus (22) folgt: 
1 In a > l ui 
“ (0,9285 — 0,1778) S m, wm 7 0, 9285 + 0,1778). 


Die gesamte Eingrenzung wird hier etwas schärfer als durch die Ungleichung (21) bzw. (22) 
allein. 


Die Eingrenzungsspanne durch jede der Ungleichungen allein ist 4yYac; hierin ist a 
die Eingrenzungsspanne der Ansätze für die Belastung p, und ce die Eingrenzungsspanne der 
Ansätze für die Belastung p,. Sind also die Ansätze für p, oder die Ansätze für p, die 
richtigen Lösungen, so wird w,om eingegrenzt durch Werte, die sich nicht unterscheiden, 
die also den wirklichen Wert geben. 


IV. Zusammenfassung. 


Die Verschiebung der Membranenmitte wird durch Näherungsansätze eingegrenzt. Die 
Eingrenzung gelingt durch geeignete Kombination der Ungleichungen, die aus den Minimal- 
sätzen folgen. Benötigt werden hierbei Näherungsansätze sowohl für die gegebene Belastung, 
als auch für eine Einzelkraft P in dem Punkte, dessen Verschiebung eingegrenzt werden soll. 
Durch geeignete Wahl von P= P,,ı erhält man die schärfste Eingrenzung. Sind die Näherungs- 
ansätze für die gegebene Belastung oder für die Einzelkraft P die wahren Lösungen, so 
erhält man den wirklichen Wert der Verschiebung. 357 


Fingrenzung von Verschiebungen und Zerrungen mit Hilfe 
der Minimalsätze. 


Von Constantin Weber in Dresden. 


Die Verschiebung eines Punktes eines beliebig belasteten Körpers wird durch 
geeignete Verbindung der beiden Minimalsätze der Elastizitätslehre ein- 
gegrenzt. Hierzu werden Näherungslösungen für die Verschiebungen und 
für die Spannungen gebraucht, sowohl für die gegebene Belastung, als auch 
für eine Einzelkraft in dem Punkte, dessen Verschiebung gesucht wird. Die 
Lösungen für die Einzelkraft müssen im Angriffspunkte die richtige Singu- 
larität besitzen. Ebenso können Linearkombinationen verschiedener Ver- 
schiebungen und folglich auch Dehnungen und Spannungen eingegrenzt 
werden. 


Ich zeigte in einer früheren Mitteilung‘), wie mit Hilfe der Minimalsätze der Elastizitäts- 
theorie der Weg einer Kraftgröße eingegrenzt werden kann, und in einer späteren Arbeit’), 
wie die Methode zu ergänzen ist, um die Verschiebung des Angriffspunktes einer zweiten 
Kraftgröße zu bestimmen. Ich nahm an, daß hiermit die praktische Verwendung der 
Minimalsätze in dieser Richtung erschöpft ist. Es ergab sich aber, daß noch weitere Fragen 
hiermit gelöst werden können. 


So zeige ich in dieser Arbeit, wie bei beliebiger Belastung und Stützung eines 
elastischen Körpers, für den das Hookesche Gesetz gilt, die Verschiebung eines beliebigen 
Punktes eingegrenzt werden kann. Weiter wird die Methode auch für die Zerrung (Dehnung 
oder Schiebung) in diesem Punkte abgeändert, so daß hiermit auch die Spannungen näherungs- 
weise bestimmt werden können. Wenn ich mich auf den zweidimensionalen ebenen Zerrungs- 
zustand beschränke, so geschieht es nur mit Rücksicht auf die Darlegung. 


1) Z. angew. Math. Mech. Bd. 11 (1931), S. 24 und 25, 
2) Z. angew. Math. Mech. Bd. 18 (1938), S. 375 bis 379. 
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I. Die Minimalsätze der Elastizitätstheorie. 


Ein zweidimensionaler Körper (ebenes Zerrungssystem mit e,—0) sei beliebig gestützt 
und belastet, siehe z. B. Bild 1. Die Komponenten der Oberflächenkräfte auf das Ober- 
flächenteilchen do seien Zdo und Hdo. Wir nehmen ein beliebiges Verschiebungssystem 
u,v,w (w=0, das,=0) an, das mit den vorgegebenen Randverschiebungen im Einklange 
ist. Diese Randverschiebungen in den Stützflächen mögen gleich null sein. 

. i ur du 0% 
Aus den angenommenen Verschiebungen werden die Zerrungen = \_, =, 
0x ey 

du 0% 
Yey= , € 
sy 92 
energie 


=Yy,x=Yzy—0 bestimmt. Im Teilchen dedydz ist dann die Zerrungs- 


1 
adedyde=6 et) H Net + rudedyde 


aufgespeichert. Für den ganzen Körper erhält man durch Integration die Zerrungsenergie 
A=\\\adxdydz. Infolge der Randverschiebungen wird die potentielle Energie der 
1: 


äußeren Kräfte um den Betrag \(Zu+ Hv) do, verringert, wobei nur über diejenigen Ober- 
09 
B 

flächenelemente do=doy zu integrieren ist, für die die Belastungskräfte gegeben sind. 


Die wahren Verschiebungen seien u, tw, %o—=0 und die daraus zu berechnende 
Zerrungsenergie sei A„. Dann folgt aus dem Satze vom Minimum der potentiellen 
Energie’): 


A— \(Zu+ Hv)doy > Au — \(Zu.+ Hv.)d og. 


OR O1 


Für die wahren Verschiebungen ist A,. gleich der Hälfte der Verringerung der potentiellen 
r er . r ” ” (o a1 So . 
Energie der äußeren Kräfte, so daß man die Ungleichung in folgender Form schreiben kann: 


A— \(Zu+Hv)dog> 2 | (Zu.+ Hv.)dog . ee I 


OR 0 


Nunmehr nehmen wir ein beliebiges Spannungssystem 0%, 0,, Try; 0: = U(0x + 0,), 
x; —=Ty;,=0, das die Gleichgewichtsbedingungen erfüllt und mit den gegebenen Randkräften 
im Gleichgewichte ist. Im Teilchen dedydz ist dann die Spannungsenergie 


u 


1 
bdedyde=,G lot +0} +9) >; (tat +2.) dzdydz 


aufgespeichert. Für den ganzen Körper erhält man 


B=\\\bdxdyda. 
E 
Da die vorgegebenen Randverschiebungen gleich null sind, so ist die Spannungsenergie gleich- 
zeitig die Ergänzungsenergie des Systemes. 


Die wahren Spannungen seien mit dem Index w gekennzeichnet, ebenso die wahre 
Spannungsenergie, die infolge des Hookeschen Gesetzes gleich der Zerrungsenergie ist: 
€, 
Be = Au= 5 \(Zu»+Hv,) dog. 


_ 


og 
Aus dem Satze vom Minimum der Ergänzungsenergie (Satz von Castigliano) folgt: 


B=Bu 
oder 


a u 
B25\(Zue+Hv.)dog a ET A 


OR 
3) E. Trefftz: Handb. d. Physik, Bd. VI, Kap. 2. Math. Elastizitätslehre, S. 666 u.f. — C, Weber: Z. angew. 
Math. Mech. Bd. 18 (1938), S. 375 bis 379. Z. angew. Math. Mech. Bd. 21 (1941), S. 32 bis 42. 
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II. Eingrenzung der Verschiebung eines Randpunktes. 





IS y, 


Im Randpunkte O0 (Bild 1) soll die Verschiebungs- 
komponente in einer beliebigen Richtung infolge einer ge- / 
gebenen Belastung eingegrenzt werden. Wir legen durch 0 ZN“ 
ein Koordinatensystem, die x-Achse in Richtung der Ver- 7AmmmmnÖa« II 
schiebungskomponente, die y-Achse senkrecht dazu. Im / \ 

Punkte O0 bringen wir eine Kraft P an, die in Richtung “ pP 
der x-Achse wirkt. Bei gewissen Überlegungen denken wir \ 

uns die Kraft P auf die kleine Länge e verteilt, um un- \z -» 

endlich große Spannungen, Zerrungen und Verschiebungen Bild ı. Beispiel zur Eingrenzung einer 
zu vermeiden. Verschiebung. 

Für die gegebene Belastung und Stützung und für 
die Kraft P und derselben Stützung wollen wir erst für das Verschiebungssystem, dann für 
das Spannungssystem Näherungslösungen angeben. Die Näherungslösungen können auch einen 
oder mehrere Parameter enthalten, die noch frei wählbar sind. 

Die Näherungslösung für das Verschiebungssystem für die gegebene Belastung sei «,, 
v,,w,—=0; es sind hierbei nur die Randbedingungen in den Stützflächen zu beachten, d. h. 
es müssen dort entsprechend unseren Annahmen «,=», = sein. 

Für die Verschiebung des Angriffspunktes der Kraft P stellen wir auf folgendem Wege 
Näherungslösungen auf. Wirkt die Kraft P auf eine Halbebene, deren Begrenzung mit der 
Tangente an den Körper im Punkte 0 zusammenfällt, so erhält man mit e>0, wenn P=1 
gesetzt wird, den Verschiebungszustand: 























{1 a iu top ur 
oe B 19 2 Tai u)inr|, mit "=x” + y° 
ı’ ae +y ! Y 
ı I1 zy 1 te I 
u, 19 a vr ujarc i 
? naG|2@+y | sg AMOS . 
7 : . N N . e ? 
Das Verschiebungssystem folgt aus der Airyschen Spannungsfunktion F', —- — yarctg y 
% 4X 


In den Stützflächen erhalten wir Verschiebungen, die nicht null werden. Wir über- 
lagern darum ein weiteres beliebiges Verschiebungssystem «,, v,, das diese Randverschie- 
bungen wieder aufhebt. 

Für die Verschiebung des Angriffspunktes der Kraft P erhält man dann als Näherungs- 
lösung: 3 

up=P(u,+u), vp=Plv,-+v,). 

Jetzt wollen wir Näherungslösungen für das Spannungssystem, sowohl für die gegebene 
Belastung als auch für die Kraft P angeben. 

Die Näherungslösung für die gegebenen Belastungskräfte seien die Spannungen o;,, 0y,, 
Teyı, 9, =Ul0x, + 0y,), Trzı =Ty;ı—=0. Diese Spannungen müssen miteinander und mit den 
gegebenen Randkräften im Gleichgewichte sein. 

Für die Kraft P stellen wir wieder eine Näherungslösung auf, die aus zwei Teilen be- 
steht. Wirkt wie vorher auf die Halbebene die Kraft P, so erhält man mit e>0, wenn 
P=1 gesetzt wird, das Spannungssystem 


9) y3 
es E—- —— Te 
en a (ce +y?}’ 
2 ey? 
Oya — a 92 PIC 
a (x + y?)” 
a: #9 
Teya wi 


7 (+ a 
0,=U(0x,+ 0y,) , 
ge tys md. 
Die Spannungen folgen aus der angegebenen Spannungsfunktion F,. Die Randflächen, 


an denen keine Stützung erfolgt, sind dann allgemein nicht spannungsfrei. An der Fläche do, 
deren äußere Normale mit der x-Achse den Winkel «a bildet, müssen die Kräfte 


2 
-_ 


2,d0o=-—-- 
u 7T 





„ cos? p cos (p — a)do, 


> 


21 ; 
H,do=—- „ 608 p sin p cos (g — a)do 
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angreifen, damit Gleichgewicht vorhanden ist (9 und r — Polarkoordinaten, also z=rcosy, 
y=rsing). Da die Randflächen aber spannungsfrei sein sollen, fügen wir ein beliebiges 
Spannungssystem 9x3. Oys, Txy;s usw. hinzu, das mit den Kräften Z,do=—2,do, 
H,do=— H,do im Gleichgewichte ist. Beide Spannungssysteme bilden zusammen das ge- 
suchte Spannungssystem für die Kraft P=1, welches allgemein nur eine Näherungslösung 
sein wird. Für die Kraft P erhält man somit: 

%p=Pl(o,, +0%,), 

oyp=FPflo, +0y,)) 

Typ=Pfrzys+ Txys) 

usw. 

Nunmehr bilden wir für den Körper, der durch die gegebenen Kräfte Z,do, H/,do und 
durch die Kraft P belastet ist, die potentielle Energie. Als Verschiebungssystem nehmen wir 
die Näherungslösung v= u, +P(w, + u,), e=v, + P(v,+v,). Die Verschiebung in «-Richtung 
des Flächenstückes e in der Umgebung des Punktes O0, auf das die Kraft P verteilt ist, ist 
w=Uu,+ Plü,+ %,,). Hierin ist «,, die von uns gesuchte Größe, die den wahreı Wert der 
Verschiebung t,,„ nach einer Seite eingrenzen soll. 

Die Zerrungsenergie a enthält Glieder von der Gestalt 

(u, +Plw,+ - (2) u, du, , „m du, du, 


+2P 


dx dx dx dx 7 dx dx 


ıp: Ö "+ p: Ri) "+ 9p: du, dm, 
| 0x, 0x & 0x 0x 


Durch Ordnen der Glieder nach den Indizes zerfällt auch a in entsprechende Bestandteile: 
a=a,+2Pa,+2Pa,+P?’a,„+P?’a,+?2P?a,.: 
Durch Integration über den ganzen Körper erhält man: 
A=A,+2PA.+2PA,+P?A„+P?A,+2P?4,. 
Ungleichung (1) wird hiermit: 
Au+2PA,+2PA,+ P?A„+P?A,+2P?A, 
— \[&, (0, + P(u, + u,))+ H,(e,+ Pir, + v,))]) dog — Pu, + Pu, + Ws)] 
0 
1 , = > | H a 0 , | 1 > > 
43 9} \ [3, (ww+J4H (u, w +4; w)) T ‚vw I I (d,w+ ®; w))] O4 zi [uw +4H (Un ttosw)]- 
OR 
Wirkt die Kraft P=1 auf einen Punkt, so wird «,, und %,,„ unendlich. Um dieses 
zu vermeiden, denken wir uns P=1 auf eine kleine Strecke e verteilt. Verteilt man die 
Kraft gleichmäßig, so ist für «,, der Mittelwert der Verschiebungen in &-Richtung zu 
nehmen für alle Randpunkte, die zu e gehören. Man kann aber auch von der Vorstellung 
ausgehen, daß P=1 so auf die Strecke e verteilt ist, daß ein einheitliches «,, entsteht. 
Da der Verschiebungszustand «,, ®, die wahre Lösung des elastischen Problemes einer 
durch die Kraft P=1 belasteten Halbene ist, so wird 


Ass - | (2,0, + H,v)do — u, = — A, 
1) 


Das Integral ist über die ganze Oberfläche einschließlich der Stützflächen zu nehmen. 
. r . . ” o . . . P 
Dieses setzen wir in die linke Seite der Ungleichung ein, um «,,, das für &>0 unend- 
wird, zu beseitigen: 


Au+2PA,+2PA,+P?A,+2P’4,— P’B, 


) [3, (u, + P(u, + u,))+A,(w,+P(v,+v,))]dog + pP| (Z,u,+ H,v,)do 
0 v0 

Püw.+Pus)Z 

LT.n 
P2 \ [#; (wu+Plu,o+ %w))+ H, Wuw+tPwu+ v;w))] d OR 


O7 





1 
7 Pr P (u,, w +P (uw + Ua) ) 
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Nunmehr stellen wir die Ungleichung (2) auf, einmal für die gegebene Belastung, dann 
für die Kraft P. Die Spannungsenergien B werden entsprechend den Indizes der Spannungen 
jeweils mit zwei Indizes versehen. 


Für die gegebene Belastung erhält man: 
Ei; 
B. 25 \ (Z,uW,0+ H,v,w)dog . 
op 


Für die Kraft P erhält man 
P’(B„+B,+?2B,)Z Sp: (Us + sw) 
Durch Addition der Ungleichungen (3), (4) und (5) folgt: 
Au+t2PA,+2PA,+P?A,+2P’A,„+B.,+P’B,„+2P’B, 


_ | [&,(u, + P(u,+u,))+A, (we, + P(v,+ v,))]Jdog + P? | (Z,u,+ H,v,)do 
0% 0 
Pu, tw) SZ — . P \ [F,(w,0+ uw) + A, (du + dv w)] dog — an 


j v Mt 


Die rechte Seite besteht aus einem Integrale und dem Gliede —— Pu,,„- Beide Glieder 


sind gleich nach dem Maxwellschen Verschiebungsgesetze, das für die wahren Verschie- 
bungen gilt. Auf der rechten Seite steht folglich — Pu,,.. Teilt man die Ungleichung durch 
|\P|, da P sowohl positiv, wie auch negativ sein kann, so erhält man für die gesuchte Ver- 
schiebung: 
F., 1 . 
Yu pi jpjAut2PAa+2PA,+P?Au+2P*A,, 
* B.+ & B„+2 P’B,, \ [>; (u, +P(u, + ",)) .: H, w, +P(, + v,))] dog 


0 B 





+ P: (Zu, + H,o,)do— Pu + Pu) 


Alle Glieder der rechten Seite folgen durch Integration aus den gewählten Ansätzen für die 
Verschiebungen und Spannungen. Die Energien A,, und B,,, die unendlich werden, sind 
herausgefallen. Hiermit ist die Einführung dieser Lösungsanteile begründet. 


In allen Ansätzen können noch freiwählbare Parameter enthalten sein, die so zu wählen 
sind, daß die rechte Seite zum Extremum wird. 


Dasselbe gilt auch für die Größe P. Die rechte Seite enthält Glieder mit P”', P° 
und P*+'; durch Differentiation nach P erhält man eine Bestimmungsgleichung für P mit 
Gliedern, die P"? und P*° enthalten. Man erhält zwei Werte für P mit entgegengesetzten 
Vorzeichen. Setzt man diese Werte in Ungleichung (6) ein, so erhält man eine Eingrenzung 
von %,,;„ nach oben und nach unten. Wir werden im weiteren %,,„ noch nach einer zweiten 
Methode eingrenzen. 


Die zweite Methode der Eingrenzung besteht darin, daß wir erst die Ungleichung (2) 
für die gegebene Belastung und Kraft P zusammen aufstellen, dann die Ungleichung (1) für 
die Belastung allein und für die Kraft P allein. Durch Addition der so erhaltenen drei Un- 
gleichungen und Division durch |P| erhält man für «,, „ die Ungleichung: 

P 1 J 9p? 22 9p 

Uw IP| = 'P| Au +21 Ay + As+ B,, +27 B;; 

; (7). 
+2PB,+ P?B„+2P?B„,— \(Z,w+H,v,)dog+P?\(Z,u,+ H,v,)do— P!uss} 
RB 0 

Auch hier wird P zweckmäßig so bestimmt, daß die rechte Seite zum Extremum wird. 
Man erhält wieder zwei Werte für P, die sich nur durch ihr Vorzeichen unterscheiden. Un- 
gleichung (7) gibt hiermit ebenfalls eine Eingrenzung von %,,„ nach oben und unten. 


Von den so erhaltenen zwei Eingrenzungen nach oben und zwei Eingrenzungen nach 
unten wird man diejenigen wählen, die «,,„ am stärksten eingrenzen. 
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Ill. Eingrenzung der Dehnung in einem Randpunkte. 


Wir nehmen erst nach Bild 2 die benachbarten 
Randpunkte O (#«=0, y—=0) und O0’ («= —6, y=0) und 
bringen in O die Kraft P/ö in positiver «- Richtung 
in 0' die Kraft P/ö in negativer »-Richtung an. 

Im weiteren wird genau so vorgegangen, wie im Bi: 
vorherigen Abschnitte. Der Verschiebungszustand für u in 
die Halbebene wird für P=1: in einem Randpunkte. 


BE I. I] A (+ 6)? 
=, 2 2? +? d— a)iny I (+ Pzreer 


273 1 | » ey 23 2 )arctg Yy | 1 (2 +6) y t. ( 2 u) PR 
Tänal2erty (3 “SIT IT EC HF \eg Prjarcie 
r — + y, r? - (“e > Ö)? + y. 


Für die Kräfte P/ö sind diese Werte noch mit P zu multiplizieren. 
Die Spannungen werden 





- (1— u)In ‚) 

y 
KL E Ö | 
mit 


USW, 


2 (2° (e+0)° 

ER Ta —— 

x2 Ö zı | r* r' 4 } 
Die wahren Verschiebungen der Punkte O0 und 0’ in «-Richtung durch die gegebene 

Belastung sind «,,. und ur,"; die durch die Eingrenzung bestimmte Größe wird 


a Um Ugrıw P R 
ö IP| 
Das ist aber die durchschnittliche Dehnung zwischen 0’ uud 0. 
Läßt man ö>0 gehen, so erhält man für P=1 und für die Halbebene die Verschie- 
bungen 


nn 


Ei 2 x 
alt)? Mer y|’ 


1 y’ Y | 
2G |(x” + y?)’ ra + y? 


und die Spannungen 
A 


-3.20° 

7 ar 

Die durch die Eingrenzung bestimmte Größe wird dann 
P 


€ , 
xzou >| ? 
IPI 


2x 
— usw. 


die Dehnung in &-Richtung im Punkte 0. 


Da die Spannung in y-Richtung im Punkte O bekannt ist (sie sei o,,), so kann auch 
die Spannung 0x, berechnet werden: 


1 


Oxow = a (2G ex0w + 4 0yo) 


Hiermit ist auch o,, eingegrenzt. 


IV. Eingrenzung der Verschiebung eines inneren Punktes und der Dehnung und 
Spannung in diesem Punkte. 


Die Eingrenzung der Verschiebung in x-Richtung eines inneren Punktes wird in gleicher 
Weise wie vorher durchgeführt. Nur sind andere Ansätze für %,, d%, Und 0,,, Oyz, Try, ZU 
verwenden. Wirkt im Punkte =y=0 im Inneren eine Kraft P=1, so erhält man die 
Airysche Spannungsfunktion 


1 


1 ? 
F,= 3 |* Inr—y are  — 5 = 


2. Mr a 


Aus dieser Spannungsfunktion kann der Spannungszustand und Verschiebungszustand be- 
stimmt werden. 
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Um die Dehnung in «-Richtung zu bestimmen, bringen wir auf der x- Achse in den 
Punkten O0 mit «=0 und 0’ mit «=—ö die Kräfte P/ö wie vorher an und lassen ö gegen 
null gehen. Dann geht F‘, über in eine neue Spannungsfunktion F/, die aus F, nach Gl. (8) 
sich berechnen läßt: 
(9). 





Aus den Dehnungen lassen sich die Spannungen ermitteln: 


2(1—u) „_ 
0x = 1 > ee (Fr 


i 


U 


ext 1 == u eu . 





E u ; e , 
Um den Ausdruck e, +7, einzugrenzen, bringen wir 








e : 0) Are ar: PH 
nach Bild 3 im Punkte 0(-=5.7=0) die Kraft P/ö in .«- 
. . ‚ ) . > x ICER 
Richtung, im Punkte 0 („= - 5,4=0) die Kraft P/ö in - 76 


Bild 3. Zur Eingrenzung einer 


u Spannung in einem inneren Punkte. 


negativer «-Richtung, im Punkte 0” (x=0, y-5d = | 
a u 


die Kraft P/ö in positiver y-Richtung und im Punkte 0" (20, y= si = ö) die 


Kraft P/ö in negativer y-Richtung an. Läßt man dann ö>0 gehen, so erhält man den ge- 
suchten Ausdruck. Man kann aber auch e, und +, getrennt eingrenzen und dann den be- 
treffenden Ausdruck bilden. 322 





Über die Kompensation der Reibung durch zusätzliche 
Schwingbewegungen. 
Von @. Heinrich in Wien. 


Zur Klärung der Frage, wieweit es möglich ist, durch Überlagerung einer 
schwingenden Bewegung die Reibung zu kompensieren, wird der Einfluß 
solcher kleinen Schwingungen einer ebenen horizontalen Unterlage auf die 
als nicht geschwindigkeitsabhängig angenommene Reibungskraft, die auf 
einem sich auf ihr bewegenden, schweren Körper wirkt, untersucht, so- 
wohl für den Fall, daß sich der Körper in Richtung der Schwingungen, 
wie für den, daß er sich unter einem beliebigen Winkel gegen diese bewegt. 





Der Gedanke, die unvermeidliche Reibung, die bei der gleitenden Berührung zwischen 
zwei gegeneinander drückenden festen Körpern auftritt, durch zusätzliche Bewegungen der 
Gleitflächen in ihrer Wirkung möglichst auszuschalten, taucht meines Wissens nur in der 
Patentliteratur auf, ist aber bisher nicht zum Gegenstand einer wissenschaftlichen Unter- 
suchung gemacht worden. 

So wird in dem D.R. P. Nr. 349476') ein Erfindungsgedanke unter Schutz gestellt, der 
es ermöglichen soll, drelibare Teile derart zu lagern, daß das Reibungsmoment, das bei ihrer 
Bewegung überwunden werden muß, ein möglichst geringes ist. Die Einrichtung ist „dadurch 
gekennzeichnet, daß der gesamte auf dem drehbaren Teil lastende Druck auf zwei getrennte 
Lager verteilt ist, die gleiche Reibungsmomente erzeugen und derartig angeordnet sind, daß 
die stützenden Glieder der beiden Lager während der Benutzung des Instruments dauernd 
ineinander entgegengesetztem Sinn um ihre gemeinsame Achse gedreht werden können, und 
zwar mit je einer Winkelgeschwindigkeit gegenüber dem drehbaren Teil, die stets wesentlich 
größer ist als die größte auftretende Winkelgeschwindigkeit des drehbaren Teils“. Als An- 
wendungsgebiet ist vor allem an Meßinstrumente und feinmechanische Apparate gedacht. Das 
englische Patent Nr. 344239?) stellt den gleichen Erfindungsgedanken für die Lagerung von 
Kreiselgeräten unter Schutz. 

Durch die angegebene Maßnahme ist eine vollständige Ausschaltung der Reibungswirkung 
möglich, solange man in einem Geschwindigkeitsbereich arbeitet, in dem die Reibung sich 
praktisch nicht ändert und solange die beiden entgegengerichteten Reibungsmomente genau 


1) Firma Carl Zeiß in Jena. Einrichtung zur reibungsfreien Lagerung von drehbaren Teilen. (25. April 1925.) 


2) Improvements in Gyroseope Suspensions. (March 5, 1931.) 
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entgegengesetzt gleich sind. Es ist aber besonders die letzte Bedingung nicht immer zu er- 
füllen, auch wenn man, wie beim angeführten D. R. P., das Reibungsmoment des einen Lagers 
einstellbar macht. Wenn sich nämlich während des Betriebes die Laererdrücke und damit 
die Reibungsmomente ändern, wie dies z. B. bei Kreiselgeräten möglich ist, so ist eine dauernde 
Kompensation der Reibung auf dem angegebenen Weg nicht möglich. 

Es besteht aber die Möglichkeit, auf die stützenden Glieder der beiden Lager eine 
schwingende Bewegung zu übertragen’). Da die Reibungsmomente in den beiden Dreh- 
richtungen praktisch keine Verschiedenheit aufweisen, so ist zu erwarten, daß dadurch die 
Reibungswirkung weitgehend kompensiert werden kann. 

Es soll im folgenden die Frage nach der Kompensation der Reibung durch Überlagerung 
einer schwingenden Bewegung in verallgemeinerter Form untersucht werden. Auf einer hori- 
zontalen Ebene befinde sich ein schwerer Körper und es wirke zwischen Körper und Unter- 
lage, entgegen der Relativgeschwindigkeitsrichtung eine Reibungskraft, die in dem maßeeb- 
lichen Geschwindigkeitsbereich nieht geschwindigkeitsabhängig sein soll. Die Unterlage voll- 
führe eine geradlinige hin- und hergehende Bewegung und zwar sei die Zeit für den Hingane 
genau gleich der Zeit für den Rückgang. Reibungskraft und Frequenz der Schwingung seien 
so bemessen, daß der Körper trotz der schwingenden Unterlage keine merkliche hin- und 
herzehende Bewegung ausführt. Ist r die Schwingungsdauer der hin- und hergehenden Be- 
wegung der Unterlage, so ist, bei Gültigkeit des Coulombschen Gesetzes, wenn f den 
Reibungskoeffizienten und m die Masse des Körpers bedeuten: m-#—=f-m:g und daraus: 

: -f-g: 7°; x ist dann die Ablenkung der Masse m während eines Hinganges (oder Rück- 
ganges) in Richtung der Schwingbewegung, wenn zu Beginn die Geschwindigkeit der Masse 
gleich Null war. Es soll also f und r so gewählt werden, daß der Wert x unmerklich 
klein wird. 

Bewegt sich nun der Körper in Richtung der hin- und hergehenden Bewegung durch 
äußere Einwirkung mit der Geschwindigkeit vr, so entsteht die Frage, welche Wirkung hierbei 
die Reibung ausübt. Da der von der Reibung bei einem Hin- oder Rückgang übertragene 
Impuls nach Voraussetzung eine vernachlässigbar kleine Wirkung ausübt, kommt also für die 
Reibungswirkung nur der über eine längere Zeitdauer gemittelte Impuls in Betracht. Die 
diesem mittleren Impuls entsprechende Kraft soll als Scheinreibung bezeichnet werden. Sie 
ist demnach der Mittelwert der von der Unterlage ausgeübten Reibung, genommen über ein 
Zeitintervall, das groß ist im Vergleich zur Schwingungsdauer der hin- und hergehenden Be- 
wegung. 

Es sei W, der Augenbliekswert der Reibung und W, der Betrag der Scheinreibung. Die 
Schwingbewegung der Unterlage sei harmonisch, ihre Kreisfrequenz sei » und ihre Geschwin- 
digkeit zur Zeit f sei «. Es gilt dann: 


nn he Ss EA a ME 


Wird » in Richtung des positiven « positiv gezählt, so ist v,,=® - wu die Relativgeschwin- 
digkeit des Körpers in bezug auf die Unterlage. Nach der Definition der Scheinreibung 
gilt also: 


5) 
(2), 


wobei T eine Zeit bedeutet, die groß ist gegen die Schwingungsdauer r. Man kann T als 
ganzzahliges Vielfaches von 7 einführen, da es auf Zeitintervalle in der Größenordnung von 7 
nicht ankommt. Man kann also setzen: 

Ian 


nt EU E RN EN Fr En Be he ae EEE 
() 


wobei » eine große ganze Zahl bedeutet. Durch Einsetzen von (1) und (3) in (2) erhält man: 


u,-cosafl 
u di 


U, cos at 


Setzt man: 
ot > u De a a ae a he . (4), 


3, Vel. K. Killian: Verfahren zur Aufnahme von Luftbildern, Allg. Vermess.- Nachr. 51. Jahrg. (1939), S. 21 
bis 29 (siehe S. 27 u. >28). 
10 
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so kann man auch schreiben: 


cost 
2 
dy 


CoSs( 
! 


wenn man zur Abkürzung: 


setzt, wobei «, als wesentlich positiv zu betrachten ist. 

Es wird nun vorausgesetzt, daß die Änderungen von 'W;! und »* während der Inte- 
grationszeit T vernachlässigbar sind. Dann ist der Integrand von (5) eine periodische Funktion 
von g mit der Periode 2#. Der Wert von (5) ist also für jedes ganzzahlige » gleich. Man 
kann also etwa n=|1 setzen und erhält: 


Man kann nun zwei Fälle unterscheiden 


1. Fall: !v*| >1, oder nach (6): |» t,. Dann wechselt der Integrand von (5a) sein 
Vorzeichen im Integrationsintervall nicht. Man erhält also: W,;,—= F W,, wobei das obere Vor- 
zeichen für »* >0, das untere für v»®<0 gilt. Ist demnach die Geschwindigkeit des Körpers 
größer als die maximale Geschwindigkeit der Unterlage, so stimmt die Scheinreibung mit der 
wirklichen Reibung überein. 


2. Fall: |e*! <1t, also: |v’<u,. Nun ändert der Integrand von (5a) im Integrations- 
intervall das Vorzeichen, und zwar an der Stelle: „= arceos r*. Ferner kann man das Integral 
in (5a) folgendermaßen umformen: 

\ ©" CoSg 


'OS« 
er / 


da 
v 
Man kann also (5a) schreiben: 


arccos r* 


- W, | vo” j ’ COS « 
w | - pP + —dg 
e 7 \ e / R ’ cosgq u 


Die Integration liefert: 


I, | - arecos ‚ı N a ..M, 


wobei das Minuszeichen sowohl für »*>0 als auch für »*<0 gilt. Man sieht aus (7), daß 
die Scheinreibung nur von dem Geschwindigkeitsverhältnis (6) abhängt. Da insbesondere 
die Scheinreibung für kleine Körpergeschwindigkeiten (v* {) interessiert, soll (7) in eine 
Potenzreihe entwickelt werden. Man erhält: 

= 3 “ Bu dl = 

W, Pr Wil" + 6 (a) 
Aus (7a) erkennt man, daß für sehr kleine Körpergeschwindigkeiten die Scheinreibung deı 
Geschwindigkeit proportional ist. Sie verhält sich also wie eine Flüssigkeitsreibung. Für 
v*—=0 verschwindet sie vollständig, sie wirkt sich also, im Gegensatz zur Coulombschen 
Reibung in der Gleichgewichtslage nicht mehr störend aus. 

Deutet man W;, als Reibungsmoment, so läßt sich dieses Ergebnis in gleicher Weise 
auf das Scheinreibungsmoment einer schwingenden Lagerung übertragen. Ordnet man zwei 
Lager an, deren Schwingungen der Phase nach entgegengesetzt sind und deren Reibungs- 
momente W;, und W;, ungefähr übereinstimmen, so gilt auch hierfür die Gl. (7), wenn man 
IWi=|W,|#!W,,| setzt; überdies werden auf diese Weise die Wirkungen der einzelnen 
Lagerschwingungen weitgehend kompensiert. Für sehr langsame Bewegungen erscheint also 
die Reibungswirkung praktisch ausgeschaltet. 

Kehrt man wieder zur ursprünglichen Anordnung zurück, so ergibt sich die Frage, wie 
sich die Scheinreibung verhält, wenn die Geschwindigkeit » unter einem beliebigen Winkel 
gegen die positive «- Richtung geneigt ist. Auf das Lager angewendet bedentet dies, daß 
man auch axiale Verschiebungen zuläßt. 
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Legt man die -Achse eines ebenen Koordinatensystems in die positive «-Richtung, die 
yAchse senkrecht dazu, so sind die Komponenten der Relativbewegung » 


des Körpers 
gegen die Unterlage: 


Next u, Cosy | 
| er 522 2% : (5). 


rel 4 “u 


Die Komponenten der Scheinreibung sind dann gegeben durch: 


Wr, 





















































} 
Mit den Abkürzungen: 
) 
3 erhält man aus (9) und (10) durch Einsetzen von (8), (4) und (): 
r an 
- 1 i ®, COS ( 
W.= W a ee 
ea 2an \ , y(v; cosp)’ tv, j D , 
. | 2 
Wı,= w 3 ‚„dg he ie > > A: 
. 2an \ 2 V(v, cosp”’ try 
: 2 
Macht man wieder die Voraussetzung, daß die Änderungen von |W;|, v, und v), während 
der Zeit T vernachlässigbar klein bleiben, so lassen sich die Integrale in gleicher Art wie 
früher transformieren und man erhält: 
> r j COS ( 
Ny.= \- I — do... 2.2.2208. Ob, 
a 3 yie COS P) TrVy 
= Ww;|f "y 
Wın \ — ——de. . . Mm 534 
j x. DR } (v, cosy)”’—+ry 
u 5 
Die Integrale von (9b) und (10b) sind elliptisch. Ihre geschlossene Lösung soll hier nicht 
gesucht werden, da ja vor allem die Lösung für kleine »} und v,, interessiert. 
© 2. - . \ 3 . ° . . a 
W;x und W;, erscheinen nach (9b) und (10b) als Funktionen der beiden Variabeln v/, 
und »,. Sie haben daher die Form: 
n W, 
W,. RP | 00, 7) a ae? est MR 
7 
; W, 
" W,, SET ER RR 5% ° 
Ir T 
N Die Funktionen J, und J, sollen nun versuchsweise naclı Potenzreihen entwickelt werden, 
die nach den linearen Gliedern abgebrochen werden sollen: 
Ne? 
2 0J.10,0 0.J/x (0,0 
\ z L 0J,.V 0) 0,00) 
er I. EAN +, Tr tt 2: 0.2.2. 0), 
S- ( x [ Oy 
su (0,0 I „(0,0 
LE ‚03,0,0) ‚0J,0,0) 
N Je, t)=J, 0) + Fr, Apr +, Sr . (104). 
or, Ory 
S() 
Man erhält also aus (9b) und (Ye): 
ie r | a : 
j \ ( cosy)dyg i " COS ( 
el J.(0,0)= lim \ > a > lim \..: a dy. 
NE er. >) . (v, eosp”’—+r), >. vr cosgq| 
. u + u 


0, >» 
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Es ist demnach: 








J.(00)= lim (a—2arecose{)—=0 . . . | Eee 
vr >\ 
Ferner ist: 
Ö J; (v;, ",) \ "n d q 
Ork lex — cos’ + vi: 
u 
daher 


vy’+dg 





0.Jx (0,0) 





: h i rd gq 
= lim — - ur lim BEE! _ ec ARE 
.. > N IIe% -cosp’ te] 0 \ (+ cos’ gp)"> 
0,2 r 






Dies läßt sich umformen zu: 


4 







2) 












R 4 ryde« 
lim \ ker a - 
ya 0 j j gr 
+», 
Setzt man zur Abkürzung: 
1 B : 
[Fo ee 5 






und 
| 











so gilt für kleine v,: A’—=1 — v,,?, daher: 


. (43B), 





Es folgt also: 










9} 


\ 2k’’.dg 
(1 — k’sin’p)”: 


oJ. (0,0) 
Ork 









Setzt man zur Abkürzung*): 





Br .... 2.2200 





J.(00) .. ‚Cd m 
so ist also: “ =\ lim 2% °\ * - Wird A nach k’ entwickelt, so folgt nach Ausführung 
0 


Orr kr > 0 


der Integration: 





0 .Jx (0,0) - e | E k? sing cosy 


N) 2 K>0 k r k'* l 


12] 


wg 











mit der Abkürzung: 











\ l t EIER 
BifzB- tt .. 020.04... 
Es ist also: 
ER N N 
d vr „m 2k Br = » In | k’ > 3: le | 2. : 5 R r . (16). 


Das dritte Glied von (9d) ist nach (9b) und (9e) definiert durch: 










7 7 








Ö Pr PB (0,0) ji Ri) \ ex cosg ' li vy(v; cosy) d gq 
. ım “ ‚* —— 2 „.dy— im „* ®. u v “2]ü 
od, 0, >0 or, ; J (tx cos p) +try 2,0 : II: cosyp”’—+v, |": 
0,>) 0,20 
Daraus folgt: 
= 
oJ, (0,0) \ \ v,C0Sy:-de 
E ' im \ , ag Beh er Te 
or, 2,0. (eos’p+ rs?) 


















4) Vgl. Jahnke-Emde: Funktionentafeln, ?. Aufl., S. 124 ff., Leipzig 19933. 
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Das erste Glied von (10d) ergibt sich aus (10c) und (10b) zu: 







7 





7 D) 






j vyd« . b 2r,d 
J,09)= lim | „——— — = lim \ nn BE 
0,20 ö ira cos op)” vUy 2,0; cos’ pP + vy 
e,>V 






oder: 












) dg 


vil+r7 | sin’g 






J,,(0,0) — lim 


"y>NV ö 









man ?): 





Setzt 









de ö 
2 En In = ea 7: 
"sin? g 











i \ . 
J , (0,0) lım 2%’ |ln pt E’ı+... Bi et ri MR 
hr > (0 ı l 






Das zweite Glied von (10d) folgt aus (10e) und (10b) zu: 










0. , (0,0) d v,dg (r; — cos po)’ dg 
u.“ lim ns \ .. . 2 „.2 lim \ .* ET rar 7? 
Or, 00H ; irx cosy”’+ry >. OH cos pP}? + riFPre 






vy>NV v,>V 





9,00)  dg 


Ivy, \\eosq 





daher ist: Dieses Integral existiert nicht zufolge der Singularität des 















Die Funktion J,(r;,r,) ist also für kleine »; und »‘, in erster 





’n a 
Integranden für 9= 5. 
Näherung von vr, unabhängig, verhält sich aber für r,=0 nicht regulär. Sie besitzt, wie 
man aus obigem, nicht konvergenten Integral ersieht, an der Stelle »,—=0 eine logarithmische 
Singularität. Es muß also der Ansatz gelten: 






J,(w;,v,)=a-vy-In|v,„| +höhere Glieder. . . - . . 2.20) 





mit noch zu bestimmendem «a. Aus (20), (10b) und (10e) folgt: 







d./,, (0.0) i k oO L v,d«e 

‚\ = lim [a (in |v,|+D)]= lim |. r - Ri we; t 
0, >0 v0 Or, ; ] (r} cosy)”’—+try 

£ Tr 





or, 







v,>V 





Daraus ergibt sich für a der Wert: 






a 





! (vi — c0sy)’-dg 
a— lim \ Ki : - 81 „#2]8],] » 
v>0 ; (1+1n|r;|) ((v; cospY’ try | , 
v,>0 





oder, nach einfacher Umformung: 









. he 
a— lim \ RE 
Da ET ET 


3 


sin? \”2 
1+ =) 






Nach (13), (13a) und (13b) wird daraus: 






rt 






2) 


s—= lim \ cos’pdg 
w>o | (+ In |R’])- + RN] - (1 A® sin? p)9e| 
k >| ) 
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Mit der Abkürzung (14) erhält man: 


2 " cos?y . dy 


ZN, + In |R)- [HR fo \ 5 


Nun ist’): 


5 7 B) . " a: i 
Base q dy ; , SM «<osg 
2? 2 | 


0 
wobei E und K durch (15) und (18) gegeben sind. 
Führt man dies in (21) ein, so erhält man bei Verwendung von (12a) 


en K-E 
= 0, d-FFGH+m RN): |d+ KR 


oder mit Benutzung von (15) und (18): 
= 


Man erhält also aus (9e), (dd), (12), (16) und (17): 


W.,.= 


Aus (20) und (22) folgt: 


W,=+- :o-Inlos|+--- . . N 25, FR 


Aus (23) und (24) erkennt man, daß die Scheinreibung für kleine Körpergeschwindig- 
keiten in der «-Komponente nur von rz, in der 4-Komponente nur von r,, abhängt. Die Ge- 
setze dieser Abhängigkeit sind jedoch verschieden. Während Gl. (23) eine lineare Abhängig- 
keit darstellt und für kleine »; mit Gl. (7a) übereinstimmt, wo die Bewegung ausschließlich 
in der Schwinzungsriehtung erfolgt, ist nach Gl. (24) die Abhängigkeit nicht mehr linear. 
Es ist zwar lim W,,=0, d. h. die Reibung wird für beliebig kleine Geschwindigkeiten in 

vy„>V 
der y-Riehtung selbst beliebig klein. Doch das Anwachsen von W;,„ an der Stelle #,=0 ist 
dWı,| 2 


unendlich groß, da lim gr z W; lim |(A+In v,))!=». Immerhin ist es bedeutungs- 
. 0%, 7 


v,>\ v,>V 
voll, daß auch in der Richtung senkrecht zur Schwingbewegung die Scheinreibung mit gegen 
null gehender Querbewegung verschwindet, wenn dies auch nicht mit derselben Größenordnung 
geschieht, wie bei Bewegungen längs der Schwingungsrichtung. 


Zuletzt sei darauf hingewiesen, daß die resultierende Scheinreibung der Geschwindigkeit 
des Körpers im allgemeinen nicht entgegengesetzt gerichtet ist, sondern schräg zu ihr liegt. 
Sind «a und $ die Winkel der Körpergeschwindigkeit und der Scheinreibung gegen die .x-Achse, 
so gilt: 
W;,„ 
W;, 

0; -vy (In1vy, 1) 
ve +ryPIn|v, 


und », von null verschieden sind. 


"y 
tg a uud he 
z 


Daraus und aus (23) und (24) folgt aber: tg (P  «) ‚ das immer von nuli 


verschieden ist, wenn gleichzeitig r/ 

Die obigen Untersuchungen dürften vor allem bei der Lagerung von Kreiselgeräten prak- 
tische Bedeutung besitzen, da hierbei die der Theorie zugrunde liegenden Voraussetzungen gut 
erfüllt sind. 365 
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gut 
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Über eine Näherungslösung der gewöhnlichen linearen 
Differentialgleichung 1. Ordnung. 
Ein Beitrag zur Theorie der Nachlaufsteuerungen und verwandter Geräte, 


Von Hans 








Biückner ın Berlin. 





Für die Differentialgleichung y+-y-kla)=fi(xr) wird unter bestimmten 
Voraussetzungen über 7, insbesondere für kleine 7. eine Näherungslösung 
ron der Form u f(r + nm) angegeben, wobei u eine ausschließlich von } 
abhängige Funktion bedeutet. 











Einleitung. 


Mehr und mehr gewinnen in der Technik Fern- und Nachlaufsteuerungen an Bedeutung. 
Verwandt mit ihnen sind Apparate, die zum selbsttätigen Abgleich von gewissen automati- 
schen Recheneinrichtungen dienen. Die Mehrzalıl dieser Geräte ist in ihrer Wirkung be- 
stimmt durch eine inhomogene, gewöhnliche lineare Differentialgleichung 1. Ordnung. Dies 
sei am Beispiel einer Nachlaufsteuerung näher erläutert. 

Wir nehmen an, daß der Zeiger eines Meßinstrumentes den Weg f(t) als Funktion der 
Uhrzeit # beschreibe. Verlangt werde, daß ein zweiter Zeiger, eine Welle oder dergleichen 
kraftgebend einen Weg y(#) zurücklege, der dem Weg f(f} möglichst proportional sein soll. 
Technisch wird diese Aufgabe zumeist dadurch gelöst, daß y(f) als der Weg der Welle z. B. 
eines Elektromotors dargestellt wird, dessen Umlaufgeschwindigkeit einer ihm angelegten 
Steuerspannung proportional ist; die Steuerspannung wiederum wird der Differenz zwischen 
y und f verhältnisgleich gemacht. Demnach gilt 












/yıh) = Fit) — yıl) 
oder 


HOFER. ...  ° 






Hierin ist 4 eine durch die Dimensionierung der Nachlaufsteuerung gegebene positive Kon- 
stante. 

Wie man sieht, ist es auf dem eingeschlagenen Wege grundsätzlich nieht möglich, eine 
strenge Proportionalität bzw. Gleichheit zwischen f und y herzustellen. Doch kann man diesem 
Ziel um so näherkommen, je kleiner 4, je steiler also die Steuercharakteristik des Motors ge- 
macht wird. Daher begnügt man sich vielfach ‘mit der Feststellung, daß if) bis auf einen 
„Schleppfehler* mit fif) übereinstimme. 








Ist fit) = konst. #, also zeitproportional, so ist auch if), nachdem die Nachlaufsteuerung 
hinreichend lange in Betrieb ist, von einem Vorgang konstanter Änderungsgeschwindigkeit 
nicht mehr zu unterscheiden, und die Differenz f —- y, der Schleppfehler, ist eine konstante 
Größe. Auf diesem Umstand beruht eine Untersuchungsmethode von Nachlaufsteuerungen. 
Man gibt eine gleichförmige Bewegung fif}) in die Steuerung hinein und mißt dann den 
konstanten Schleppfehler; dies wird für verschiedene konstante Änderungsgeschwindigkeiten 
von fit) durchgeführt, und die Meßergebnisse geben einen gewissen Anhalt für das Verhalten 
der Steuerung auch im allgemeinen Falle. 

















Oitmals erscheint es aber notwendig, ganz unmittelbare und zugleich anschauliche Vor- 
stellungen über den zeitlichen Verlauf von y zu besitzen. Man könnte meinen, daß die exakte 
Lösung der Differentialgleichung (1), die bekanntlich durch Quadraturen lösbar ist, diese 
Forderung vollauf befriedigt. Dies ist aber nicht der Fall: die exakte Lösung hat namentlich 
für den Ingenieur den Nachteil, daß man ihr nieht ansehen kann, ob und wie sehr fıt) als 
dureh y(f) näherungsweise wiedergegeben gelten kann. 










Es ist daher das Ziel dieser Arbeit, für kleine Werte von / eine Näheruneslösunge der 
Gl. (1) aufzuzeigen, welche y in einen einfach zu überblickenden Zusammenhang mit f bringt. 
Diese Aufgabe soll aber nicht nur für die Differentialgleichung (1), sondern auch für den Fall 
behandelt werden, daß der Koeffizient A der Differentialgleichung eine variabele Zeitfunktion 
darstellt. Der Fall des variabelen Koeffizienten tritt bei dem vorstehend behandelten Beispiel 
dann auf, wenn die Erregung des Motors zeitlichen Änderungen unterliegt, wie das bekannt- 
lich bei Netzanschluß der Fall ist. 






Als Hauptergebnis dieser Arbeit ($-2) wird sich herausstellen, daß man der Näherungs- 
lösung die Form 


BEIDE ee ee. a a A 
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geben kann, wobei « eine ausschließlich von 4 abhängige Funktion von f ist und durch eine 
gewisse affıne Transformation aus / entsteht. Neben der Untersuchung dieser Lösung und 
ihres Fehlers beschreibt die Arbeit noch eine weitere Näherungslösung sowie eine Verall- 
gemeinerung der genannten Näherungslösung für den Fall A= konst. ($ 4). 

Ein Nebenergebnis, das selbständiges Interesse verdienen könnte, wäre vielleicht in den 
Hilfssätzen des $ 1 zu erblieken, welche sich ganz allgemein auf die Abschätzung des Fehlers 
einer Näherungslösung beziehen. 









$ 1. Allgemeine Bemerkungen zum Fehler einer Näherungslösung. 





Sämtliche Untersuchungen dieser Arbeit haben es ausschließlich mit reellen Größen 





zu tun. 
Es seien f(x) und 4 (x) zwei stetige Funktionen für a << b(a< b). Ferner sei 4 (x) 0. 
Dann besitzt die Differentialgleichung 






y+4(#)-y fi) a N de 0 N | 





für a x=_b gewiß das allgemeine Integral 





wei Bla (+10 


No 


E(s,2,)-ds) . . VRR © 








Hier ist 










Bia,2.)=1Ela,2)=0” Be ERSTE 





C und x, sind Integrationskonstanten, wobei die letztere dem Intervall <a,b > angehören 
muß. An sich ist eine dieser Konstanten überflüssig, doch ist es zweckmäßig, zunächst noch 
keine von ihnen festzulegen. Unter diesen Umständen besteht der folgende 






Hilfssatz 1: Man kann mindestens ein partikuläres Integral y der Differentialgleichung (3) 
finden, welches der Abschätzung 








Yy 9(a,b)-Max /f(®)| '), <<: <h 


genügt. Hierin ist 






8 ' E(a,b) - 1 : 
ns E(a,b) -+1| ) 







Beweis: Es ist, wenn in (4) C = gesetzt wird, 





f E(s, x) f(s): ds 


)(s) 






ya; 0, Po \ 


Xo 








in dem auf der rechten Seite stehenden Integral ist f(s) die einzige Funktion, die ihr Vor- 
zeichen ändern kann. Daher besteht die Abschätzung 







Re N. .- Safe ont I 


y(x;V, x.) = 
’ 4.(8) a<x<h 






Nun ist 





zu Recht. 





X 
| E(s,x)ds 
A (s) 






| an a Rt 






Xxo 






Nunmehr verfügen wir über =, in der Weise, daß die Gleichung 





1—E(x,, 4) 1+E(x,,b) 





erfüllt wird; daraus folgt, indem diese Gleichung mit F (a, x,) multipliziert wird, 





2-E(a,@x,)=1+E(a,b), 


l E(b,«a) E (a,b) l 
1+Kk(b, a) E(a,b) +1 






1--E(x,, a) (8). 








1) Offenbar existiert das Maximum wegen der Stetigkeit von f(x). 
®2) Der Faktor # läßt sich im Falle A = konst., f  konst. nieht verkleinern, ohne daß der Hilfssatz falsch wird. 
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Aus dieser Rechnung erkennt man auch, daß es eine Stelle x, der verlangten Art geben 
muß; denn einerseits verschwindet die linke Seite von (8) für x,= a, während sie für „,— b 
das Vorzeichen der rechten Seite und einen größeren Betrag als diese annimmt, und anderer- 
seits ist die linke Seite der Gl. (S) eine in x, stetige Funktion. Da die Funktion 1 E(x,. x) 
in .e monoton ist, was ohne weiteres aus (5) folgt, so gilt unter Berücksichtigung von (8) 
E(a.b) l 

E(x,.%) E(z,,a Ei(x,.b) I (a.b) Be, 
Coy A r,,d) ) ElabW+1 (a,b (NM). 

Aus (6), (7) und (9) ergibt sich jetzt der Hilfssatz 1. 
Es sei u (x) eine für a<.r<b erklärte, stetige und stetig differenzierbare Funktion. die 

wir in die linke Seite von (3) einsetzen. Das gibt 


u +/(la)W=g(x) ’ (10). 
Wir setzen z=y— u und ziehen (10) von (3) ab. 


z+4(x)-2"= f(x) - y(&) ’ . (1). 


Das allgemeine Integral von (11) hat die Form 


2 = yIx;,0,%,) - ul). 
Wendet man nunmehr Hilfssatz 1 auf (11) an, so gewinnt man den 


Hilfssatz 2: Man kann mindestens ein partikuläres Integral y der Differentialgleichung (3) 
finden, dessen Differenz gegen eine vorgegebene stetige und stetig differen- 
zierbare Funktion u (.r) der Abschätzung 

you (a,b): Max |f- y O(a,b)- Max |f—u- Als) w 
a<x<h a v<b 
genügt. 

Die im Hilfssatz 2 vorkommende Funktion f- uw — 4-w' bezeichnen wir im folgenden 
als das Defizit der Funktion « (x) bezüglich der Differentialgleichung (3). Sie stellt den 
Rest dar, bis auf den die Differentialgleichung durch u (=) befriedigt wird. Das Maximum 
des absoluten Betrages des eine stetige Funktion bildenden Defizits kürzen wir im folgenden 
ab mit 

I)(u;a,b)=Max|f—u-/4-w|, a 2<b, 


Das Defizit einer Näherungslösung kann zur Abschätzung ihres Fehlers benutzt werden. Um 
dies darzulegen, schätzen wir die Differenz zwischen einem vorgegebenen partikulären 
Integral y(#; €, x,) und der als Ersatz für dieses Integral genommenen Näherungslösung u (x) 
ab. Indem wir noch ein dem Hilfssatz 2 genügendes partikuläres Integral einführen und mit 
y(#;C’,#,) bezeichnen, dürfen wir schreiben: 

y(z;C,2,) - y(2;C’,x,)=Ei(x,,2)-(C 


y(2;C,x,) a =yle;l’, (a) + Eia,,.r) (ua ) C’)+ E(z,,2) - (C— u(&,)) (12). 


0 


y(x;C,x,) -uke)| S(+E(z,,x))- Alu;a,b)- (a,b) + E(x,,X): IE u(,) 


Damit haben wir den 
Hilfssatz 3 (Fehlersatz): Bedeutet (x) ein beliebiges partikuläres Integral der 
Differentialgleichung (3) und « (x) eine als Ersatz dieses Integrals vorgesehene 

stetig differenzierbare Funktion, so gilt die Abschätzung 
ya E(#,,2)-|yla,) - ul)! + Hu;a,b)- Ola,b)-(1+E(x 


| +) a<x<h. 

Hierin ist x, eine beliebige Stelle des Grundintervalls <a,b > und I(w; a,b) 

das Maximum des absoluten Betrages des Defizits von «(x#) im Grundintervall. 

Durch den Fehlersatz ist die Abschätzung des Fehlers (y — u) der Näherungslösung auf 

die Abschätzung der Differenz y— uw an einer einzigen Stelle und auf die Berechnung oder 

Abschätzung von I(w;a,b) zurückgeführt. Die Abschätzung des Fehlersatzes läßt sich noch 

verschärfen. Wir setzen Max (#,x,)=x,, Min (x,x,)=x,. Dann gilt, was an sich ohne 
weiteres aus dem Hilfssatz 3 folgt, 

ya) ul) E(x,.%): |y(a,) - u,))+ Ilu;a,,2,)- I(&,,2,)(1+E(x,.x)) : . (3). 


Nun ist auf Grund von (8) 


(I+E(&,,2)): 9 (#,, & 
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Somit 
y(z) - ul) E (x,,%) : |!y(&,) -ul&,)j+|1 Ei, Aa) > ir AR: 





Daß die Abschätzung (14) schärfer als die Abschätzung (12) ist, ergibt sich aus den Ungleiehungen 













(uw; x,,%,) I(u;a,b) und 9(&,..72,)< 9 (a,b), 








die zu verifizieren es keiner besonderen Schritte bedarf. 
Die praktische Bedeutung der Größe I(w;..) erkennt man besonders auch aus folgen- 





den Beispielen: 
I. Aa) >0; w(a)=yle); z a. Die Absehätzung (14) liefert 










ds 
ie) 





X 





(u; a,.x) I(u:a,b). 





y—u K(u;a,.x)-(1 exp| \ 






Stimmen also das partikuläre Integral und seine Näherung in den Anfangswerten überein, so 
liegt ihre Differenz dem Betrage nach unter I(w;a,b). 

2. f(x), u (x), A(#) seien für O< x <» definiert; ferner sei 0<A< A=konst. Wir setzen 
a=x,=0. Dann ist, ähnlich wie beim vorhergehenden Beispiel, 






ul) y(d) u(0) - EO,#) + IKu:V,.r). 





y() 












X 


Nun ist E(0,.«): exp | \ 





ds ‚ 5 . 4 
T ) e ' und somit E (0, x) eine mit wachsendem . monoton gegen 
)(s 








u 
Null abnelımende Funktion, so daß es eine Stelle x=x, geben wird, von der ab der erste 
Summand auf der rechten Seite der Fehlerabschätzung praktisch gestrichen werden kann. 
Näherungsweise können wir also setzen 










Yy(x) u(x) 9 FB. 





Sollte I(u:0,2) 0, d.h. « ein partikuläres Integral der Differentialgleichung (3) sein, so 
ist die letztere Abschätzung so aufzufassen, daß die Differenz y -u von der Stelle », ab 
praktisch annulliert werden kann. Dieses Beispiel trifft auf alle in der Einleitung genannten 
Geräte zu. Die Größe ., ist dabei als eine gewisse Anlaufzeit zu betrachten, nach deren 
Ablauf der Fehler der Näherungslösung der vorstehenden Ungleichung genügt. 

3. Die Differentialgleichung (3) sei, wenn f 0 gesetzt wird, in übersichtlicher und an- 
schaulicher Weise lösbar. Doch möge es wünschenswert sein, die als unübersichtlich an- 
genommene exakte Lösung des inhomogenen Falles durch eine Näherungslösung zu ersetzen. 
Hier genügt es, als Näherungslösung eine einzige Funktion u(x) zu finden, von der lediglich 
verlangt wird, daß sie schlechthin irgendein bestimmtes, aber nicht vorgegebenes, partikuläres 
Integral der Differentialgleichung (3) ersetzt. 

In diesem Falle kann man zur Fehlerabschätzung auf den Hilfssatz 2 zurückgreifen, 
und man erkennt, daß die Größe I(w;a,b) eine obere Schranke für den Fehler bildet. 

Der Fehlersatz und die vorstehenden Beispiele 1 bis 3 lassen es nunmehr gerechtfertigt 
erscheinen, bei den angekündigten Näherungslösungen die Fehleruntersuchung auf eine Be- 
trachtung des Detizits zu beschränken. 
















S$ 2. Die Näherungslösung und ihr Defizit. 






l. Es sei A (x) der im Grundintervall <a,b > als stetige Funktion definierte Koeffizient 
der Differentialgleichung (3). 
Wir ordnen jetzt der Funktion 4 (x) eine Funktion « (.) zu, die der Gleichung 







H(a,u)=4Ale+ u) + u= EEE U N N 






genügen soll. Diese Gleichung läßt sich nur für solche Wertepaare (.r, «) erfüllen, für welche 
die Zahl z=x+ u dem Grundintervall angehört. Man gewinnt eine Parameterdarstellung 
sämtlicher Lösungen («, x), indem man setzt 


i 







u = A (z) | 





(16) 







z=2—- u—=2z+i(e) ] 






und z das Grundintervall durchlaufen läßt. Durch die stetige Transformation & = 2 +4 (z) 
wird das Intervall <a,b > auf ein Intervalla<x<ß abgebildet, so daß also die durch (15) 
definierte Funktion (x) für alle Stellen des Intervalls <a, > und nur für diese erklärt ist. 
Indes muß z. (*) nicht notwendig eindeutig sein. Die Eindeutigkeit ist aber dann gesichert, 
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wenn die Intervalle < a,b > und <a, > eineindeutig vermöge der Transformation x —z + Az) 
aufeinander abgebildet werden. Dies letztere ist bekanntlich der Fall, wenn i(z) stetig 
ix 


ö ( 
differenzierbar und d 


0 ist. 


Von nun ab setzen wir die stetige Differenzierbarkeit von Ä sowie das Bestehen deı 
Ungleichung 
l 


voraus, wodurch die Eindeutigkeit von «(.#) gesichert ist. Mit dieser Voraussetzung wird 
jetzt auch die in x und « stetige Funktion H (#, «) nach jeder der beiden Variabelen stetig 
differenzierbar, und es gilt 

oH 


Q 14 


0, 


wenn 2—+ «=2z im Grundintervall liegt. Nunmehr liefert der bekannte Satz über implizite 
Funktionen das Ergebnis, daß die Funktion «(#) in ihrem Definitionsintervall <a,$> als 
Lösung der Gl. (15) sowohl stetig, wie auch stetig differenzierbar ist. Wir haben 


-(14 u’ (&)) 
oder 


(17), 
insbesondere also gilt 
(17). 


Ist die Funktion A (z) zweimal stetig differenzierbar, so trifft dies auch für « (.«) zu, und es ist 


(15). 


u(lz+4(z)) +4 (z t i (19), 


welche wegen (17’) nach 4 eindeutig im Intervall a,b auflösbar ist Zwischen A und u 
besteht also eine gewisse Reziprozität. Man gewinnt «a aus 4 nach dem gleichen Verfahren, 
welches 4 aus u ableitet. 

Die die Reziprozität begründende Darstellung (16) läßt erkennen, daß nicht nur durch 
2=2+4(z), sondern auch durch z=.r.—+ u (x) dieselbe eineindeutige Abbildung der Inter- 


> 


valle <a,b> und <a,$ > aufeinander hergestellt wird. 
Die Endpunkte des Intervalls <«a,#> werden von den beiden Zahlen @a+%4(a) und 
b+4(b) gebildet. ‚Welche dieser beiden Zahlen den Punkt @ und welche den Punkt ? liefert, 
hängt von dem Vorzeichen von 1+4'(z) ab. Ist 1+4’(z2)>0, also 4’ > - 1, so bestehen die 
Entsprechungen 
a=a+sla); P=b+i(b); 


denn die die Abbildung zwischen den Intervallen <a,b > und <a, > vermittelnde Funktion 
2=2-+4(z) nimmt mit wachsendem Argument monoton zu. Nimmt sie dagegen monoton ab, 
also 7’: l, so vertauschen sich die Entsprechungen. 

Wir setzen von jetzt ab noch A (x)=F0 voraus. Dann können die Intervalle <a, b 
und <a, > keinen gemeinsamen inneren Punkt besitzen, wenn 4’ (.r) I ist. In der Tat 
würden aus der Existenz eines solchen Punktes die einander widersprechenden Ungleichungen 


oder 


oder 
i(b) 0 
hervorgehen. 

Damit also die beiden Intervalle <a,b> und <a, > mehrere gemeinsame Punkte 
haben, muß die Bedingung 4’ (x) > —- 1 notwendig erfüllt sein. Um eine hinreichende Be- 
dingung zu finden, nehmen wir, um etwas Bestimmtes vor Augen zu haben, A (#) als positii 
an. Dann muß sein 

d« (A « bh: p . $ > > . a s . (20), 
und hieraus folgt 
a a+ila)<b £ } a A ET. (21). 

Die Bedingung (21) ist also notwendig für die Existenz eines gemeinsamen Teilintervalles 

von <a,b> und <a,?>. Sie ist aber auch hinreichend, und man erkennt, daß das nieht 
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nur aus einem Punkt bestehende Intervall <a=a-+-4(a), b> sämtliche Stellen enthält, an 
denen 4 (x) und u (x) gemeinsam definiert sind. 
Im Falle A<0 tritt an die Stelle von (21) die Bedingung 






a< pP h +4(b) EIS UP RE re» a Re A) 





dafür, daß die Intervalle <a,b> und <a,ß > ein gemeinsames Teilintervall enthalten. 
Der Definition von «u nach (15) liegt ein einfacher geometrischer Sachverhalt zugrunde 
der an Hand von Bild 1 erläutert sei. Hierin stellen X 
d und Y die Achsen eines rechtwinkligen kartesischen Ko- 
ordinatensystems dar, in welchem die Kurve E mit den 
Endpunkten @, und 9, die Funktion Y=4(X) wieder- 
eibt. Legt man nun durch den Punkt X = x der Abszisse 
die Gerade ® (x) mit der Gleichung 



























Y=—-XA-+r, 





so ist die Frage nach der eindeutigen Lösbarkeit von 
(15) mit der anderen äquivalent, ob und wann © (.) 
einen und nur einen Schnittpunkt mit E hat: denn einem Lösungselement (z,.) von (15) 
Punkt P mit den Koordinaten 


Y: A, X 2 u 










können wir einen 


















eineindeutig zuordnen, wobei P sowohl auf © (#) als auch auf E liegt. Somit stellt «u (x), 
soweit eine solche Lösung eindeutig vorhanden ist, die negative Ordinate des Schnittpunktes 
von E mit © dar. Nach einer anderen Interpretation ist «.(@) dem Abschnitt von ® pro- 
portional, der von der X-Achse und von E aus dieser Geraden herausgeschnitten wird. Wählt 
man daher ein affines Koordinatensystem X, Z, welches von der X-Achse einerseits und der 
Geraden © (0) als Z-Achse andererseits gebildet wird, so ist die Gleichung der Kurve & bis 
auf einen konstanten Faktor Z=u(X). Somit entsteht « aus 4 durch eine affıne Trans- 
formation. 

Die Funktionen 4 (x) und «u (.) sollen wegen der zwischen ihnen bestehenden Reziprozität 
als miteinander gekoppelte oder als Koppelfunktionen bezeichnet werden. Ihre wichtigsten 
Eigenschaften, wie sie sich aus vorstehendem ergeben, fassen wir zusammen in dem 














Kopplungssatz: 







Es sei A (z) eine im Intervall a=_ 2 b(a<b) zweimal stetig differenzierbare 
Funktion. Ferner sei 7’ (2) l. Dann ist die Gleichung 
(a +) +u=V 






in einem Intervall «==. eindeutig nach .« auflösbar. Die Intervalle 








a=z<b und a<.r<p werden durch die Transformationen #=2+4 (z) 
und z=x+ u(x) eineindeutig aufeinander abgebildet. Es ist u(r) eine 
zweimal stetig differenzierbare Funktion, es gilt 1. (#) =+ l und es besteht 






die zur Definitionsgleichung von „ analoge Relation 






u@+/)+i=0, 







welche in a<z=<.b eindeutig nach 4 auflösbar ist. Unter der Voraussetzung 
+ >0,a-+4(a)<b sind die Koppelfunktionen im Intervall <a=a+ti(a),b 
unter der Annahme A<0,a<b+4(b) im Intervall <a,b-+i(b) > zugleich 






definiert. 






2. Unter allen Voraussetzungen des Kopplungssatzes ist 

a rn a MR 
eine für a<x<P definierte stetige Funktion. Sie ist die angekündigte Näherungs- 
lösung. Da jedoch eine Näherungslösung nur insoweit von Interesse ist, als sie auch im 
Grundintervall oder einem Teil davon erklärt ist, beschränken wir die Variabele = der 
Funktion « (x) auf den gemeinsamen Teil <a,,b, > der Intervalle <a,b> und <a,ß 
Wie dieser Teil im einzelnen gelegen ist, steht im Kopplungssatz. Für a, <x<b, ist auch 
u (@-4+4(x)) definiert, denn es gehört #+4(x) dem Intervall <a,$#> an. Das gleiche gilt 
auch für w(e+9-4), 0<9<1. Dies vorausgeschickt, schätzen wir jetzt I(u;a,,b,) ab 
unter der Voraussetzung, daß f zweimal stetig differenzierbar sei. Dann ist auch « (x) zwei. 
mal stetig differenzierbar, und wir haben 











u 2 E 2: RR 2 
u(2c+4(e)) =u (8) +4(a) - w (x) + = wW(z+9-.4), 0O<d<t. 
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Ferner ist wegen (19) 


“(ce +4(e)) = fe +i+u(le +] = fe), 
somit 


ee . Ar{B) .,, z1= 
ul) +4: u (x) f(x) : (ve +%%) 


I(u;a,d5)<—>. Max 4’(z)- Max |w” («) 


ns a, <xı<sb, TH 


Weiter eilt wegen (17) und (18) 


u (af (at -( 


| : 
I(u;a,5)<>- Max 4(&)- Max 
a; \ bi a > 


Für den Spezialfall A = konst. erhalten wir 


I(u;a,,b,) . 4° Max If" do) te a: 
“ a=<ı=y 
Somit ist in diesem Falle für kleine |A. der Fehler eine ganze Größenordnung kleiner als '} 
Man wird aber in vielen Fällen erwarten dürfen, daß ähnlich günstige Verhältnisse auch bei 
variabelem 4 auftreten, da / mit seinem Quadrat in die Abschätzungen (24) und (25) eingeht. 
In $ 4 wird ein Beispiel hierfür angegeben werden. 


Die Näherungslösung «= f(x-+ u) entsteht aus der Funktion f(x) durch eine Verschiebung 
des Arguments «. In diesem Zusammenhang sei an den Spezialfall erinnert, daß die Funktion 
f(x) bei konstantem 4 eine Sinusschwingung beschreibt. Bekanntlich existiert dann als parti- 
kuläres Integral der Differentialgleiehung (1) wieder eine Sinusschwingung, die gegen die 
erste phasenverschoben ist. Um eine solche Phasenverschiebung handelt es sich auch bei 
unserer Näherungslösung, wenn man x als die Phase auffaßt. Aus diesem Grunde trägt die 
folgende Aufstellung der wichtigsten Merkmale der Näherungslösung die Bezeichnung 


Phasenverschiebungssatz (l. Art): Die im Intervall <a,b > definierte Differential- 
gleichung y+4(x)-y’ = f(x) besitzt unter sämtlichen Voraussetzungen des 
Fehler- und des Kopplungssatzes die Näherungslösung 
ul) f(x + 1) 2 
worin . die mit 4 gekoppelte Funktion bedeutet. Die Näherungslösung entsteht 
aus f(x) durch eine Verschiebung der Phase x um die Koppelfunktion (x). 
Im gemeinsamen Definitionsintervall <a,,b, > der Koppelfunktionen besteht 
die Abschätzung 
2” (3) “ 
I(u;a,,b,) Max a’ Max |e”() 


"hun Fu 2 u ft Zu 
wobei in u’ (x) die ersten beiden Ableitungen von / und f eingehen, voraus- 
gesetzt, das f’ und f” stetig existieren. 


Würden wir statt der Differentialgleichung (3) die neue Gleichung 
y+ıu-y=f(«) 
lösen wollen, so bekommt man auf analoge Weise die Näherungslösung 
" f(« +4), 
und es gilt eine zu (24) und (25)sanaloge Fehlerabschätzung. Wir haben damit den 
Satz der verbundenen Differentialgleichungen: Sind in den Differentialgleichungen 
y+iy=fla), y+tuy=f(#) 


die Funktionen A und z. miteinander gekoppelt und zugleich in einem ge- 
meinsamen Intervall erklärt, so sind die Gleichungen in dem Sinne mit- 
einander verbunden, daß u=f(#“—+ u) eine Näherungslösung der ersten, 
"—f(=+4) eine Näherungslösung der zweiten Differentialgleichung im ge- 
meinsamen Definitionsintervall der Koppelfunktionen ist. 
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3. Die im Phasenverschiebungssatz angegebene Funktion « ist nicht die einzige, mit 
der man das in der Einleitung gesteckte Ziel erreichen kann. Eine weitere Funktion dieser 
Art ergibt sich als Lösung der Differentialgleichung 


ut+4-w 5, En i ls (28) 


Ein Teil ihrer Lösungen läßt sich durch den Ansatz 


1) 


.)Jds, a AR Fa A 


erfassen. Gemäß den dem Hilfssatz 1 vorausgegangenen Ableitungen kann man .«, im Inneren 
des Grundintervalls so wählen, daß « (.) seinem Betrage nach nicht größer als das Maximum 
von | (x)| ist. Da u(#) stetig ist und an der Stelle #—=.r, verschwindet, muß ferner ein 


0 
Teilintervall <a,,b, des Grundintervalls existieren, für welches «+ u(#) nur Werte aus 
dem Grundintervall annimmt. Somit ist dann u=f(# + u) für a, <x<b, als Näherungs- 


lösung wohldetiniert. Für das Defizit erhält man 
u +4. W"— f=f(r- (Hu): f’(e 
fi u f (2 1) f(«) 


ni” 
9-0), 0<#, 
Somit 
I(u:a,,b,) FT EEE EAN ne TR 


197 1) +.) 


= Mi vb; a xv=h 


Diese Abschätzung erscheint einfacher und besser als die Abschätzung (25). Indes be- 
steht bei der zweiten Näherungslösung der Nachteil, daß « gemäß (28) nicht so einfach be- 
stimmt werden kann wie die zu 4 gekoppelte Funktion. Trotzdem ist aber die zweite Näherungs- 
lösung immer noch einfacher zu überblicken als die exakte Lösung der Differentialgleichung (3), 
da « ausschließlich von 4 abhängt. Dieses Ergebnis fassen wir zusammen im 


Phasenverschiebungssatz (2. Art): Die im Intervall <a.b > definierte Differential- 
gleichung „+4 -y = f(#) mit fals zweimal stetig differenzierbarer und A als 
stetiger, von Null verschiedener Funktion besitzt die Näherungslösung 
w—f(#—+ u), wobei «#« durch die Differentialgleichung u -+4- wu i defi- 
niert und so gewählt ist, daß sein Betrag das Maximum des Betrages 
von / nicht überschreitet, und daß ferner wu=f(x-+u) in einem Teil 
a,b, des Grundintervalls definiert ist. Es gilt 

4’(z) 
I(u;a,,b,) Max =. ur If (a) 
HMEX= bh 


bi u a X 


Die Näherungslösung des vorstehenden Satzes ist unbestimmt, weil « nur bis auf den 
Umstand festgelegt sein soll, daß sein Betrag das Maximum von / nicht überschreitet und 
daß u (#) für einen gewissen Teil des Grundintervalls erklärt ist. Man wird deshalb « von 
Fall zu Fall geeignet wählen müssen. Z. B. dürfte es zweckmäßig sein, im Falle A =konst. 
u ) zu setzen. Dann nimmt die Näherungslösung die Form u —=f(xr — 4) an. Sie entsteht 
aus der Funktion f(.) dureh eine konstante Verschiebung der Phase. Die von ihrem Defizit 
abgeleitete Größe wurde schon angegeben. Im übrigen liefert auch der Phasenverschiebungs- 
satz 1. Art diese Lösung, und zwar nur diese. 


S 3. Beispiele. 


l. Nachlaufsteuerung. Es bedeute x die Zeit und 4 (“) eine positive Funktion 
von der Dimension einer Zeit. Somit ist «(.) eine negative Funktion mit der Dimension 
einer Zeit. |. (.), entspricht einem vor dem Zeitpunkt = angenommenen, also veralteten Wert 
von 4. Die Näherungslösung «= f(.r + u) gibt solche veralteten Werte von f(x) wieder, die 
schon zur Zeit « — || angenommen wurden. Somit wirkt die Nachlaufsteuerung, 
wennihr Verhalten dureh die Näherung riehtig beschrieben wird, im Sinne 
einer Verzögerung der ihr zugeführten Funktion. Je kleiner die Verzögerung, 
desto genauer die Nachlaufsteuerung. Ist A konstant, so auch die Verzögerung. Pendelt A 
in nicht zu großem Ausmaß sowie mit kleiner Geschwindigkeit und Beschleunigung um einen 
mittleren Wert, so gilt dies auch für die Verzögerung. Zur Diskussion des Fehlers werde ein 
praktisches Beispiel vorgelegt: 


3) Auf diese Definition hat mieh auch Herr M. Krafft besonders hingewiesen. 
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konst. VOL see, 0 f(x) 6000 E (inheiten 
f(x) 160 E/see., f: (.r) SO E!sec*. 


Die Nachlaufsteuerung werde zur Zeit a 0,01 sec in Betrieb genommen, so daß die 
Näherungslösung # = f(# — 0,01) für „> 0 definiert ist. Wir benutzen nun die im 2. Beispiel 
des $ 1 angegebene Fehlerabschätzung. 


y(r) u(.r) y(0) u(0) - E(O.#) + Ilu:V, x) 


1 
y(0) RE a Max 
f 


da 


x 
yla)- ul) 6O00-e " "-0,0001.40) Einheiten. 


Spätestens eine Sekunde nach Einschalten der Steuerung ist, grob gerechnet, der erste 
Summand der rechten Seite verschwindend klein gegenüber dem zweiten und der Fehler deı 
Näherungslösung kleiner als 0,004 Einheiten 

2. Elektrischer Derivator. Wir behandeln die in Bild 2 beschriebene 
Kondensatorschaltung, bei welcher ein Kondensator und ein Widerstand in 
Reihe liegen. Liegt an der Schaltung die Spannung 4(f), so gilt für die am 
Widerstande liegende Spannung V die Differentialgleicehung V+T:-V-g-T. 

Hier bedeutet "= WC eine Zeitkonstante Für V ergibt sich näherungsweise 
VI T:g(t - T). Für kleine 7 mißt V also den Differentialquotienten von 4, 
aber nicht den aktuellen, sondern den um die Spanne T veralteten. Die Schaltung 
ist als Derivator brauchbar. Die Fehlerabschätzung der Näherungslösung erfolgt 
wie beim 1. Beispiel. Zu beachten ist, daß 5 (N in die Abschätzung eingeht. 


Differentialgleichung mit variabelem Koeffizienten Es sei von 


-2 näherungsweise irgendein partikuläres Integral zu finden. 


Ferner 


also auch 4’ | I. Da ferner 4 zweimal stetig differenzierbar ist, sind somit alle Voraus- 
setzungen für den Phasenverschiebungssatz 1. Art gegeben. Die Gleichung A (“+ u) 

ist für #« > 1 eindeutig nach u auflösbar. Es ist u(#)=1/#. Damit erhalten wir für « 

die Näherungslösung «=.r +1/.r. Zur Abschätzung von I (w:..) benutzen wir (24) aus $ 
Es ist 

1 ö ” 
> Max 4#’(a)- Max |ju” («) 
- mx hy x } 


I(#;a,,b,) 


Es sei a >2 und 5b soviel größer als a, daß a<b+i(b) gilt (dies läßt sich wegen 
lim 4(@=)—=0 bei hinreichend großem b erreichen). Dann ist a,=a, a=a+tila) Wir 
—>% 
haben dann 
j 7 (a) 
I(u;a,b,) 


2 3 
(a +ıla)) | A 


») 


Für a=4 ist I(n;4,.b,)< 0,0014; dabei hat (4) den verhältnismäßig großen Wert A (4) 


nen Wir vertauschen jetzt die Rollen von 4 und « auf Grund des Satzes der verbundenen 
ein Differentialgleichungen. Also lautet die Gleichung y-+y'/#— x, und die Näherungslösung ist 


u 54 r Il. Auch hier erhält man I(u:4,b, 0014. 
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$ 4. Verallgemeinerung der Näherungslösung für den Spezialfall = konst. 








- 


Nimmt man in der Gleichung 







y+4-y—=f(x) die Substitution 







En ii r = konst. Er a . &D) 





vor, so gewinnt man die neue Gleichung 
fi Br f(«) ach . 


IE: I+r-4 












Sie hat nach dem Phasenverschiebungssatz 1. Art die Näherungslösung 











“A sl« h. 


Durch Umkehrung der Substitution (31) folgt daraus als Näherungslösung für die ursprüng- 


1 -+ ri) 






liche Gleichung 









r4 






I+r4/ 






4 e 





u (ir) flx N 


ii Fri 





Ist diese Näherungslösung im Teil <.a,,b, > des Grundintervalls definiert und ist ferner die 
Funktion f(.«) zweimal stetig differenzierbar, so sind folgende Abschätzungen und Umfor- 






mungen gestattet. 














Max f («) 


a=zsXI=s 


I(u:a,,b,) 





—_ 
4 
_ 
— 
be 


h, 


Indem wir noch die Abkürzung 







DZ 





einführen, können wir (33) auf die Form bringen 







I(u:a,,b,) Max e".(1—r DEZE m)+ | Pl m)! fi) ... 4. 


da sXx—b, 





Hierbei kann man noch setzen 





P p ’ m’ [24 ( ( 
f)=fle - m) +mn fie m) +5 (a —9,m), 0<9,<1, 












rm 
e 





lI+rm-+r?’: m’. e!:""12, B<» <1, 











(1i—rm)-e"—1= rm’ Il + (rm — 1): e!:"" 12%, 





rm 


me m—=m’-Ir+r’-m:e:"" 12) 





Dann geht (34) über in 












v,m)-5 





Verkleinert man m seinem Betrage nach, so wächst das Intervall < a,,b, >, in welchem 
die Näherungslösung zugleich auch für das Grundintervall definiert ist. Strebt m gegen Null, 
so konvergiert <.a,,b, > gegen <a,b >. Aus der vorstehenden Abschätzung folgt: 





——. Alt:a,,B,) eh vr 
lim uud Max | f@)2 -r fo) +2]. 
= 





m>V0 m a= b 

Die verallgemeinerte Lösung kann besser sein als die Lösung der Phasenverschiebungs- 
sätze. Man erkennt dies z. B. aus dem Spezialfall f(.) = e'*, in welehem die verallgemeinerte 
Lösung auch ein partikuläres Integral der Differentialgleichung (3) im exakten Sinne ist. 
Durch passende Wahl der Konstanten r hat man die Möglichkeit, den Fehler möglichst klein 
zu machen. 

Im übrigen könnte man auch die verallgemeinerte Näherungslösung auf den Fall eines 
variabelen 4 übertragen. Doch dürfte dann eine zu komplizierte Form entstehen. 349 
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Synthese ganzer rationaler Funktionen. 
Von H. Rosenhamer in Västeräs. 


Zur Abdämpfung von Übergangserscheinungen pflegt man Grenzfälle zu 
konstruieren, die einer verhältnismäßig einfachen mathematischen Behand- 
lung zugänglich sind, worauf man die Beiwerte der vorgelegten Aufgabe 
derart abändert, daß der wirkliche Vorgang in ein (rebiet positiver Dämp- 
fung hineingeschoben wird. 

Man kann nun aber auch umgekehrt von einer erwünschten und dem vor- 
gelegten Falle angepaßten Übergangsform ausgehen und dafür sorgen, daf 
die hierfür erforderlichen Bedingungen tatsächlich erfüllt werden. Bei 
stetiger Übergangsform setzt dieses Verfahren vor allem die Kenntnis des 
Bildungsgesetzes ganzer rationaler Funktionen bei gegebenen Realkom- 
ponenten der im allgemeinen komplexen Nullstellen voraus. 


Bei der Untersuchung stetig verlaufender Übergangserscheinungen gelingt es meist, die 
sich abspielenden Vorgänge durch eine lineare Differentialgleichung mit konstanten, reellen 
Beiwerten zu beschreiben, sei es, daß das zu untersuchende System von Hause aus linear 
ist, sei es, daß man die Untersuchung auf so kleine Störungen des Dauergleichgewichtes be- 
schränken kann, daß sich die auftretenden Kräfte als lineare Funktionen dieser Störungen 
darstellen lassen („Metlıode der kleinen Schwingungen“). Das Integral einer solchen Differential- 
gleichung setzt sich additiv aus einem Dauergliede, verursacht durch die aufgeprägten Kräfte, 
und einer der Ordnung der Differentialgleichung entsprechenden Anzahl von zeitlichen 
Exponentialfunktionen als Übergangsgliedern zusammen. Soll sich ein den aufgeprägten 
Kräften entsprechender Dauerzustand einstellen, so müssen die Exponentialfunktionen ab- 
klingen, die Realteile der Exponenten also alle negativ sein. Da sich die Exponenten als 
Produkte aus der Zeit und den Wurzeln der der Differentialgleichung zugehörigen algebrai- 
schen Gleichung, der Hauptgleichung, berechnen, müssen demnach sämtliche reelle Wurzel- 
komponenten der Hauptgleichung negativ sein. 

Ein allgemeines Kriterium, das stets ermöglicht für eine algebraische Gleichung be- 
liebigen Grades anzugeben, ob deren sämtliche reelle Wurzelkomponenten negativ sind, hat 
Hurwitz auf Anregung durch Stodola ausgearbeitet‘). Das Verfahren, welches man als 
„Methode der Dämpfungsprüfung“ bezeichnen kann, läßt jedoch keine Schlüsse auf den Be- 
trag der Wurzeln oder darauf zu, in welcher Richtung dieser bei einer Veränderung der 
Gleichungsbeiwerte beeinflußt wird. 

Einen Schritt weiter führt die „Methode der Dämpfungserzeugung“, als deren Vertreter 
die „Methode selbsterregter Schwingungen“ genannt werde’). Sie ermittelt für einen geeignet 
herausgegriffenen Teil des gegebenen Systems die Kraft ®,, welche er zur Ausführung von 
ungedämpften Schwingungen benötigt; diese Schwingungen beeinflussen ihrerseits den rest- 
lichen Systemteil und veranlassen ihn zur Ausübung einer Kraft ®, zurück auf den ersten 


sy 


I} . > Dn . a 
Teil. Je nach der Seite, auf der man beim Durchlaufen der Kurve v in der komplexen Ebene 


S 


1 
den Punkt (1, j - 0), der ja Gleichheit der beiden Kräfte und damit dem wirklichen Schwingungs- 
zustande des Gesamtsystems entspricht, liegen hat, arbeitet dieses gedämpft oder ungedämpft. 
Darüber hinaus läßt sich aus dem Aufbau von ®, und ®, recht schön der Einfluß der ver- 
schiedenen Systemgrößen auf das Zustandekommen einer Dämpfung erkennen. Über den 
Betrag der Dämpfung gibt aber auch diese Methode keinen Aufschluß. 


Was man letzten Endes wünscht, ist, eine unter Berücksichtigung der jeweils gegebenen 
Verhältnisse als am vorteilhaftesten erkannte Übergangsform genau vorschreiben zu können; 
mittels möglichst geringer Veränderungen an den Beiwerten der vorliegenden Differential- 
gleichung mit Hilfe von zusätzlichen Beiwerten, die ja im allgemeinen doch aufzuwendende 
Zusatzdämpfungsapparatur bedeuten, muß mit anderen Worten eine möglichst wirksame 
Dämpfung erzielt werden. Man wird allerdings nur in den seltensten Fällen eine viel 
kräftigere Dämpfung erzwingen, als sie durch die reziproke Zeitkonstante der größten Maschine 
des Systems gegeben ist, da die sonst für die Abdämpfung erforderliche Kraftreserve ein 
Vielfaches der für den Dauerbetrieb erforderlichen Kraft betragen würde; auch die Rück- 
sicht auf die zulässigen Beanspruchungen, beispielsweise bei Spannungsregelung die Rücksicht 
auf Isolationsfestigkeit, zwingt zu dieser Begrenzung. Weiterhin ist es wenig vorteilhaft, 

1) A. Hurwitz: Math. Ann. Bd. 46 (1895), S. 273; — A. Stodola: Schweiz. Bauztg. Bd. 22 (189), S. 113; — 
W. Hort: Technische Schwingungslehre,’ Berlin 1922, S. 164. 

?) A. Leonhard: Die selbsttätige Regelung in der Elektrotechnik, Berlin 1940, S. 59. 
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große Unterschiede in den Dämpfungsexponenten zuzulassen; deren Summe ist nämlich, wie 
später (erste der Gln. (8a)) ersichtlich wird, gegeben, und der Dauerzustand wird doch nicht 
eher erreicht, als bis das langsamste Übergangsglied verschwunden ist. 
Bei dieser „Methode der Dämpfungswahl“ liegt somit die Aufgabe vor, die Beiwerte 
Bi er (1) 
einer ganzen rationalen Funktion n-ten Grades 
ga)=a,+x-0a, +x’-a,+---+2"7'-a._,t2". . ....0.2) 


zu bestimmen, deren reelle Nullstellenkomponenten gegeben und insbesondere alle einander 
gleich sind, wobei zusätzlich die Bedingung zu erfüllen ist, daß alle Beiwerte reell sein 
müssen; bezüglich der imaginären Nullstellenkomponenten bestehen keine besonderen Wünsche. 


Sind die Nullstellen der Funktion (2) 
er RR (3), 
so kann (2) bekanntlich auch in der Form 
mm ee) 2) a)... arte 
geschrieben werden und Ausrechnung sowie Vergleich mit (2) gibt dann 
— = &n +4,-,= En, Ans Un, | 


we An-a = &n, '' (— IP —=E, | 
x n 


(5); 


hierbei bezeichne das Symbol x„ die Summe derjenigen Nullstellenprodukte, die durch die 
“ 
möglichen Kombinationen der n Nullstellen zur #-ten Klasse ohne Wiederholung dargestellt 


n I az n i 
werden, so daß die Indizes Bi den Kombinationsanzahlen # entsprechen. Es bezeichne also 


BER + +: +0 
=2,%,+2.,% +: +, +%,%,+°° +&%n-ı%n 
2 
=. 0,0%. +: +, +9. %t °F on .n-ı%n 





= dt En 


Nullstellen sind im allgemeinen komplex und es sei 


x. =v,tjm, 
,=u,+tjm, 


In = +Jj Wy 


Setzt man (7) in (5) ein, so ergibt die Ausrechnung mit Hilfe des Schlusses von n—1 auf n 
folgende zwei Gleichungsgruppen, von denen die erste durch Vergleich der reellen Glieder 
mit den Beiwerten der ganzen rationalen Funktion (2) entstanden ist, die zweite durch Null- 
setzen der imaginären Glieder: 


— IdAyn-, = Un Wy 
8 

+ 44-3, = 0, }n — OU Wn 
2 0 0 2 


— Ay-3 = U MW — Un Wy 
3 1 2 


\x use . , , ’ 
(— 1 an_-x = Una — Un Nut Un mn +: 
x 0 “”—2 2 a | 4 





n —ı ‘ D) , 
1", ZU Wn— dy My t Un Wu — + 
ns n-ı 3 n-ı 4 








) 


); 
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Vu w n 
0 1 





Un Wy ==& 
1 1 





Un Wu — Un My —0 
Bi 1 0 3 





On Wa — Un Mt Un Wan — +: =0 
ı es 3 win 5 


_ ku 












Un Wu — Un wn+ Un Wy +: =0 
n-ı 1 n-s 3 n—5 5 












Das Produkt v, w,„ ist wiederum symbolisch zu verstehen und bedeutet die Summe aller hin- 
2 4 

sichtlich der Indizes wiederholungsfreien Produkte aus den möglichen Kombinationen von n 

Elementen ® zur (<—A)-ten Klasse mit den möglichen Kombinationen von n Elementen w 

zur 4-ten Klasse; insbesondere bedeutet hierbei 







„enm=l | 
0 0 


9). 
mon=0 | RT 


<o <o 






Die Anzahl der Summanden des Produktes v„ m, beträgt 


a 4 


ETF 


Einige einfache Beispiele mögen die Berechnung von »„ w„ verdeutlichen: 


ih % 


















v,n,=v,+v,+Pv, 
1 0 






vYw=vV,W,%,+0,%,%W, +%,W,@, 






1 2 
vn, =v,0,%+v,v,%w,+v,v,W, 
2 1 









vn, =v, (m,w, ww, +n,w)+tv,(w,w, +w,w, + mw, w,) 
1 2 






+7, (ww, tw, w,tw,w,)+tv, (w, w, + w, w, + w, w,) 


vw,—v,v,W,W, + v,v0,w,@w, + v,uw,@, + v,v,w,w, + VV,Ww,W, + VV,Ww,m,. 
2 2 














Bisher ist unberücksichtigt geblieben, daß die komplexen Nullstellen stets konjugiert 
auftreten®). Beachtet man diesen Umstand, so läßt sich 









Wn-ı 


0 (10) 


= —-n., MN, '* m= 















schreiben, wobei für w„ der obere oder der untere Wert zu wählen ist, je nachdem n gerad: 
zahlig oder ungeradzahlig ist. Die in den Gln. (8a) auftretenden Produkte ®, w„ bekommen 
damit folgendes Aussehen: PTRERE 






Un Mn =+ du 
“ 0 % 





PER 2 : 2 
Un Wa = — Un-3Wı — Un-z a W5 —''' 
3 =>, _ iu 






+ Un—ima4 (v, > v,) (vw, E v,) w,w, + Un-1,256 (v; RT, v,) (v, Pr v,) w,Ww, nahe 
“ “ 4 






(11a). 


, pe ‘ 2 „m? . a! 2 
Un Wy =+ Un-1,2,3,4 Wı Ws; un Un-1,2,5,6 Wı W; + ke 
E 4 Zu, “4 


“- 





. . . . 2 am 
Un- 1,21 3,4510 (d3 x v,)(v, De v,) w, mw, Ww, — 
4“ 













3) Serret-Scheffers: Lehrbuch der Differential- und Integralrechnung, Bd. 1, Leipzig und Berlin 1924, S. 598. 
13° 
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Die Produkte v„ w,„ der GlIn. (8b) werden dagegen: 
„-2 A 


U Wa = — Ya-3(d; —Y)W, — Un-.1ıld — 9); — 
as 1 “2 3 


u A . , „2 
Fr =+v, 120 (0 —-v,)w,n, +» 
J ” 


14 


®n 1,2» 3, 4, 5» 6 (v, — v,) (u - v,) (%, — v,) w,Ww,W, — 


. en . . . . 2 m? 
Yu NMy—=— Yy-,230 56 (dı — U) WWW, — 
“”—5 5 „6 





Die in den GIn. (11) vorkommenden Indizes n—a,ß,... sollen bedeuten, daß von den 
n Elementen 1,2, ...n die Elemente «a,ß, ... zu streichen sind; beispielsweise wird daher 
v,w,=—- (tv, tr, +V)wi (w, +0, +v,+v,)n; 
1 2 
(v, + V, + Vz; + v,) w;, 


vw =—- (vr, +0, v0, +90,%, +v,v, +v,v%, + vd, wi 
2 
(nv, +tr.,+t+e0, +90, +0,07, +0, 0)093 
9, +v,, +, vv, +,9%, +90, +20, dV) w; 
+(v, — v,)(w, —v,)w, w, + (v, — v,) (d, — v,) w, wm, 
+(v, —v,) (vd, — d,) mW, W,. 
Mit Hilfe der Gln. (8) und (11) lassen sich die Beiwerte einer ganzen rationalen Funk- 
tion n-ten Grades aus den im allgemeinen voneinander verschieden reellen und den paarweise 


entgegengesetzt gleichen imaginären Nullstellenkomponenten berechnen unter der Bedingung, 
daß die Beiwerte reelle Zahlen sein sollen. 


Eine bedeutende Vereinfachung ergibt sich, wenn man, wie oben angedeutet, alle reellen 
Nullstellenkomponenten gleich groß macht. Setzt man nämlich 


leur Yan Veen, 7 BE Er a EEE |. 


in das Gleichungssystem (11a) ein, so geht dieses über in 


[nn 
Un Wu= + | ) 
%” 


v x 


Ar. N — 
U N, =— v* , b, ( 
”—2 2 u 


he N - 
Un Wu =+ v”* ,( 
„a 4 x“ — 





wobei zur Abkürzung gesetzt ist 
2 2 
,=w, wm tn; +: +mwir_, 
1 


,=wW, =niny tmin; +. +mim tn; nm 4. +wr_,wi_, 
2 


mit n"=n bzw. n—1 im Falle, daß n geradzahlig bzw. ungeradzahlig ist. Das Gleichungs- 
system (11b) wird durch Einsetzen von (12) identisch Null, was bedeutet, daß die Beiwerte a 
stets reelle Zahlen werden, wenn man alle reellen Nullstellenkomponenten gleich groß macht. 
Nunmehr berechnen sich die Beiwerte folgenderweise: : 








end 
[57 
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ı— 2 
2 )+d: 


2 \ 
ve b, (" n ")+r ».(" ‘) 


(15). 


n .——_ u 
ran) te + 


x 








Die reelle Komponente sämtlicher Nullstellen wird durch den Beiwert mit der höchsten 
ÖOrdnungszahl bestimmt zu 
Ay-ı 


RE. Eee ee re 





mit diesem Dämpfungsfaktor verschwinden also alle Übergangsglieder. Nach Wahl von » 
können die Beiwerte a gemäß (15) berechnet werden; die Werte b kann man dabei beliebig 
wählen, wenn man nur die folgenden Bedingungen erfüllt: 


a ee er a ie 5 
wi; a a sa ee I 


(17a) ist eine unmittelbare Folgerung aus (17b) und nur mit Hinblick auf die praktische 
Anwendung besonders erwähnt. Zur Erfüllung der Forderung (17b) ist die Kenntnis der 
Bestimmungsgleichungen für w,,,,,.. notwendig. Die Bestimmungsgleichung für w},,,,,. kann 
aber sofort angeschrieben werden, wenn man beachtet, daß die b,,,,... nichts anderes als die 
Beiwerte einer ganzen rationalen Funktion mit den Nullstellen »?,,,... sind, wie ein formaler 
Vergleich von (14) und (5) zeigt: 
n' " 1 
bu (— 1)’+ mw! b, ‚| 1)? +. end, — nn tu" —=0. . . (18a). 


+) 





‚ bestimmt sich dann aus 


mw 


Ir 2». 








== %+ je ou „2 ... 18} 
+ Vmwt, + Yw}, RUHR RO u  - 
w, 


w, 
w, 


Hiermit sind alle Hilfsmittel gegeben, um die Beiwerte einer irgendwie zustande 
gekommenen Hauptgleichung durch geeignete Ergänzungen so gestalten zu können, daß eine 
wirksame Dämpfung erzielt wird. In der Natur der für die Beschaffung der Zusatzdämpfungs- 
glieder zur Verfügung stehenden Anordnungen liegt es, daß eine solche gewöhnlich mehrere 
Beiwerte erfaßt; dennoch kann es vorkommen, daß die Reichweite einer Zusatzdämpfungs- 
einheit den Forderungen nicht genügt und der Einsatz von zwei oder noch mehr Einheiten 
nötig wird. Regelungsanordnungen, die eine Beeinflussung des Beiwertes a„_,, also des 
Dämpfungsexponenten gestatten, werden üblicherweise als Schnellregler bezeichnet. Eine 
Heranziehung des Beiwertes a, oder allgemeiner von Beiwerten, die auch als Faktoren von 
Dauergliedern auftreten, zu Dämpfungszwecken muß vermieden werden, da sonst geringe Ab- 
gleichungsfehler starke Störungen des Sollzustandes verursachen können. 

An einem Zahlenbeispiel soll die Anwendung des Gesagten auf den Fall einer Funktion 
sechsten Grades gezeigt werden. Für die Übergangsdrehzahl habe sich bei einer Leonard- 
schaltung mit stetiger oder auch mit Mittelwertsregelung der Drehzahl ohne Vorhandensein 
von Zusatzdämpfung eine Differentialgleichung fünfter Ordnung ergeben. Die erforderlichen 
Zusatzdämpfungsspannungen sollen über Siebtransformatoren beschafft werden, womit Trans- 
formatoren gemeint sind, die, eingangsseitig an Mischspannungen angeschlossen, die Ände- 
rungen derselben übertragen, während die Gleichspannungsanteile lediglich die Sättigungs- 
zustände beeinflussen. Durch Einsetzen dieser Siebtransformatoren in das ursprüngliche 
System erhöht sich die Ordnung seiner Differentialgleichung um zwei Einheiten, voraus- 
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gesetzt, daß die Transformatoren alle dieselbe Hauptfeldzeitkonstante haben; sorgt man 
außerdem auch noch dafür, daß die Streufeldzeitkonstanten hinreichend klein sind, so 
kann man sich mit einer Ordnungserhöhung um eine Einheit begnügen. Die derart für die 
Übergangsdrehzahl y erhaltene endgültige Differentialgleichung laute 


ya+ty"”a+y”a,+y”a,+y’a+ty”’a+ty®=0 .. 0... (19a). 
ihre Hauptgleichung also 
u.trz20a, +20, +0, ra, +ea+x=0. .... 0.0... (19b) 
und die Beiwerte seien 
a, = 4INO er a,=3836-+ a, 5° 
a, = 1%62 " Fa a,=202-+a,s” re re 
a,= 37+a,85"* a,= 32+qa,S" 


hierbei entsprechen die Zahlenwerte den Verhältnissen ohne Zusatzdämpfung, während die 
Zusatzgrößen a, durch die Siebtransformatoren geliefert werden und nun geeignet zu be- 
messen sind. 


Zur Berechnung der für eine möglichst günstige Abdämpfung der Übergangsdrehzahl 
erforderlichen Beiwerte a stehen die nachfolgend zusammengestellten Beziehungen zur Ver- 
fügung. Die sechs ersten der Gln. (15) geben 


a, 
—ıVUV 


(20); 


b 


b, 


a + 


v® \ v? v* 





gemäß den ersten drei der GlIn. (14) bedeutet dabei 


=mwy=nitn;+ m} 
1 


‚en; =wWin; twin; tn; 


= mw; =w’nim) 
3 


die b und w sind gemäß den beiden Gln. (17) an die Bedingungen gebunden 
b,.;5>0, reell. 
EEE © 
zur Berechnung von w;,,,, hat man gemäß Gl. (18a) die Beziehung 
+) +) =0 
und die w,,.,..., endlich berechnen sich gemäß den Gln. (18b) zu 


w, / 
=+yYn}, 


w,__ en - 
w, a 2 HUN | 


w, w, 

Da a, und a, hier fest gegeben sind und sich durch die Siebtransformatoren nicht be- 
einflussen lassen, wird man zunächst die sich aus diesem Umstande ergebenden Beziehungen 
ermitteln. Man erhält sie aus den beiden letzten der Gln. (20) zu 





YUM 
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2 Mi FR 12 u: 3 24 ’ 
> < u ae . . . . . . . . . . . +a 5 
re 
b a a IE 
en o un 1 Rah 1 o < “ 
Fer T + „ ee Ve a ar 


„rss b b h i 
damit sind yi und ze durch = ausgedrückt und die Gl. (23a) enthält jetzt nur mehr einen 
willkürlich veränderlichen Beiwert. 

Die Wurzeln der Gl. (23a) müssen gemäß (22a) positiv und reell sein. Hierfür geeignete 
Beiwerte sind im allgemeinen durch Probieren zu bestimmen. Bis hinauf zu Gleichungen 


a. f ; \ 
zehnten Grades in e entsprechend Differentialgleichungen zehnter und elfter Ordnung — 


lassen sich jedoch die in Betracht kommenden Gebiete durch Wurzelausziehung dadurch ein- 
i ; ah . mw? £ 
engen, daß man die Extremstellen berechnet und sich auf diejenigen Werte Re) beschränkt, 


für die alle Extremstellen reell sind; bis hinauf zu Gleichungen achten Grades kann man 
überdies die Gebiete reeller Nullstellen unmittelbar berechnen. Im vorliegenden Falle erhält 
man zunächst durch Differenzieren der Gl. (23a) und Nullsetzen 


‚2 / 
Wo 


m?\? 
- 2+| )-3=0 A a a ee inne A, 


v* v° v? v” 


und hieraus für die Extremstellen 


n® _ db, 1 ar 3-45 Be 3 ran) SE 





v A Tine ve. 3 oe 8 
die drei reellen Nullstellen sind somit sicher auf das Gebiet 


b 2 3 e 
„® >—3+ V 2 Pe en >70) 
beschränkt, wie aus (25b) unter Benutzung von (24a) folgt. Berücksichtigt man noch die Be- 
L ap i ; 
dingungen e >0 sowie (—v)>0, so findet man, daß der Dämpfungsexponent innerhalb 


des Bereiches 


'/a, Bi 


ee ADELS 


zu wählen ist; Einsetzung des Zahlenwertes gibt 


FREE 
u. 6 I 9,05 Ss . 


5 
0<-o<] 


: : b 
Die Grenzen, innerhalb deren Pr gewählt werden muß, ergeben sich aus (22a), (25c) 


und (24a,b) gemäß 





0 
2 \ 6, «€ 
= j - v® 2 Pe b, % - & ap 3 6 
Br ir v: 954 Ds De a a (26). 
u , 
mir > Me“ Men 





ü : b ; : i 
Um nun die für 2? zulässigen Werte genauer einzugrenzen, kann man für einige, inner- 


halb des Gebietes (25d) zu wählende Parameter (—v) sowie für unabhängig Veränderliche 





2? ‚ welche innerhalb des Gebietes (26) liegen, das Verhältnis 
—p° b . 
Frl) RR ge ae ar Be Pa rer Ber RS ME 
berechnen, wo 
En 
P=-(5)+ 5 a a 2 
. 3 EL - 
a=-(4-5) 2. 
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3 
) 2 . .ı 01° N r . x 
bedeuten; Fr >1 entspricht dann dem casus irredueibilis der Gl. (23a) und somit den Ge- 


a b i : a a 
bieten, in denen Fr zu wählen ist. Die Durchrechnung für (—v)=2 1,5 und 1s°' gibt drei 
Bereiche, die um so breiter sind, je kleiner (— v) ist. Ein berechneter Übergangsverlauf wird 
sich daher mit um so größerer Sicherheit herstellen lassen, je kleiner man (—v) wählt, das 
aber andererseits mit Rücksicht auf eine schnelle Abdämpfung wiederum möglichst groß sein 
soll. Als geeignetes Wertepaar kann man 

(-v)=15s"', 
nehmen und erhält hierzu gemäß (24a) und (24b) 
b, b. 
Z=114, — =3163. 
v* ne 


Die Wurzeln der Gl. (23a) berechnen sich schließlich in üblicher Weise zu 
w” N e. RR 
a 16,04, 55,30, 3,55. 

Zusammenfassend hat sich also ergeben: Eine Drehzahlübergangsform, welche nicht nur 
die gegebene Differentialgleichung (19a), sondern auch die an erstere zu stellenden Ansprüche 
befriedigt, wird erhalten, wenn man folgende Werte für die Zusatzbeiwerte herstellt: 

a.,=(15+75:.6+1141) 19 37= 78835°%, 
a.,= (20 +75-4): 1,5° — 386 = 6A: 
a,,— (15 + 75) - 1,5° — 202 — 08°, 


a,=6-15— 32 =— Be, 


Die Übergangsdrehzahl wird dann gleichmäßig mit der Zeitkonstanten 0,67 s abgedämpft 
und setzt sich aus drei Komponentenpaaren zusammen, die periodisch mit den Frequenzen 
fi. 210s"', f,,>218s' und f,.>0,5s"' verlaufen. 

Im Falle ganzer Funktionen zweiten und dritten Grades kann man sich auf sehr 
einfache Weise einen Überblick über sämtliche Nullstellen und insbesondere über die Lage 
aller Nullstellen mit gleichem Realteil verschaffen. 

Führt man zunächst in die Funktion zweiten Grades 

yg)=au, tra +... . ee ars, var Ge 
die Größen 

2 a, 
== Br; ei, 35 
+y ia| R +y|a 


ein, so kann man (28a) in der Form schreiben 


. 
x 
s 


(29) 


9 (x) 


a, 


u,x:P . (28b), 
wobei die oberen oder unteren Vorzeichen zu verwenden sind, je nachdem a, größer oder 
kleiner als Null ist. Die Nullstellen der Funktion (28b) sind 
y r P} 
. R a, JI@.N _ . 
E=utjr= "72V wıri. . Be * 
wo die beiden Vorzeichen unter der Wurzel sich wie vorher auf das Vorzeichen von a, be- 
ziehen. Reelle und imaginäre Wurzelkomponente u bzw.» sind in Bild 1 als Funktion von a, 
dargestellt. 


In die Funktion dritten Grades 


NEE TFTLERHE :..: 8 2 ern rn a ER 
kann man die Ausdrücke 


(32) 


. (31b). 





Z. angew. Math. Mech. € I 2 E 
Bd. Nr. 3 Juni 1942 Rosenhamer, Synthese ganzer rationaler Funktionen 





























aaa 


se 


a 
































| 
RR | 





Bild 1. Die Nullstellen <=«#+ jr der ganzen ratio- 


nalen Funktion zweiten Grades 15a, +5? (2=a,=0). 


Bild 2b (rechts). Kleinstwerte der reellen Nullstellen- 
komponenten « der ganzen rationalen Funktiondritten 
Grades 1+5a, +5?@3 +53 im negativen Grenzbereich. 
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Bild 2a. Größtwerte der reellen Nullstellen- Bild 2e. Imaginäre Nullstellenkomponenten +» 
komponenten « der ganzen rationalen Funktion der ganzen rationalen Funktion dritten Grades 
dritten Grades 1+5«, +5?aa3 +53 im negativen 1+5a, +5?a,+S53 im negativen Grenzbereich 
Grenzbereich. der reellen Nullstellenkomponenten. 
Die Nullstellenkomponenten u und » sind nur von a, und a, abhängig und lassen sich daher 
in einer Ebene darstellen, was in Bild 2a und 2b für « und zwar für dessen negativen 
Grenzbereich, der gemäß dem oben Gesagten hier hauptsächlich interessiert, geschehen ist. 
Der Bereich nur negativer Realkomponenten wird begrenzt durch die Kurve, auf der 
zwei der reellen Komponenten verschwinden, 4, =u,=0, und die gemäß den drei ersten 
Gln. (8) dargestellt wird durch die gleichseitige Hyperbel 
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1 
u = Er ee Ne see te 
2 


Alle drei reellen Nullstellenkomponenten sind gleich groß, u, = u,= u,= u, auf der Kurve 


A,\” 3 h 
=($).2+.=w-2-- TE PR BET EEE WERTEN 3 0. (SE, 


wie sich mit Hilfe der drei ersten Gln. (15) ergibt. Der Bereich, in dem alle drei Nullstellen 
reell sind, der also dem casus irredueibilis entspricht und in dem die Kurve (34) keine 
Gültigkeit mehr besitzt, hat als Grenze eine Kurve, die durch „= u, 1, ="r,=»v,=0 ge- 
kennzeichnet ist und durch die Gleichungen 


. (85a), 


. (35b) 


=- 
u 
beschrieben wird. 

Außer diesen Grenzlinien sind in Bild 2a auch die Schichtlinien für die größten u-Werte 
eingezeichnet, die also den am langsamsten verschwindenden Übergangsgliedern entsprechen 
in dem Falle, daß die kubische Gleichung die Dämpfungsexponenten einer Übergangserscheinung 
bestimmt, während Bild 2b die restlichen w-Werte enthält. Für die Berechnung dieser 
Schichtlinien im komplexen Gebiete erhält man mit «,=u, aus den ersten drei Gln. (8) 


1 
a,= 4, (a; + u) —7 EEE NT Tin 
3 


sowie 
. (36b). 


Man zeichnet demnach am einfachsten zunächst die Geraden (36a), welche bereits die u,- 
Linien darstellen; auf diesen trägt man dann gemäß (36b) die «,-Punkte ein und erhält durch 
Verbinden der Punkte gleicher Werte die «,-Linien. Im reellen Gebiete bestimmen sich die 
Linien auf Grund der Bedingungen », =r,—=»v,—=(0, was für jede der drei reellen Kompo- 
nenten 4,, 4, 4, eine Gleichung von der Form (36a) gibt, so daß die «-Linien des reellen 
Gebietes durch eine dreifache Schar von Geraden wiedergegeben werden. 

In Bild 2c sind noch die Schichtlinien der zugehörigen imaginären Nullstellenkompo- 
nenten aufgezeichnet. Sie bestimmen sich gemäß der Gleichung 


r 
a,= +) -— 0:27 en DEE a: 


ji 


4 
tz 


die man aus den ersten drei Gln. (8) mit «,= «, erhält und in der »» =—»v,=» gesetzt ist. 
Auch diese Schichtlinien lassen sich sehr einfach konstruieren, wenn man jeweils auf einer 
4,-Geraden die verschiedenen Werten » entsprechenden a,-Punkte gemäß Gl. (37) einträgt 
und dann die Punkte gleicher »-Werte auf den verschiedenen u,- Geraden miteinander 
verbindet. 

Die oben für die Dämpfung ausgewählten Nullstellen liegen im Falle der quadratischen 
Funktion auf der Geraden u,= u, des Bildes 1, im Falle der kubischen Funktion auf der 
Kurve „== u, der Bilder 2. 


Zusammenfassung. 


Der Zusammenhang zwischen gegebenen komplexen Nullstellen einer ganzen rationalen 
Funktion »n-ten Grades und ihren Beiwerten läßt sich unter Verwendung einer symbolischen 
Schreibweise in übersichtlicher Form darstellen. Besonders einfache Beziehungen erhält man, 
wenn sämtliche reellen Nullstellen als gleich groß vorausgesetzt werden. 

Dieser letztere Fall entspricht der günstigsten Abdämpfung von Übergangserscheinungen, 
die sich durch zusammenhängende lineare Differentialgleichungen mit konstanten Beiwerten 
beschreiben lassen, da ja die Summe der Dämpfungsexponenten hier gegeben ist, ihre Auf- 
teilung in gleich große Komponenten somit die schnellstmögliche Abdämpfung verspricht. 
Ein zufriedenstellender Verlauf der Übergangserscheinung kann derart bei mäßigem Rechen- 
aufwand genau vorausbestimmt werden. Die Anwendung dieses Verfahrens wird an einem 
Zahlenbeispiel für eine Differentialgleichung sechster Ordnung ausführlich behandelt. 

Zwei Bilder geben einen Überblick über die Lage der hier in Frage kommenden Null- 
stellen im Falle der beiden einfachsten Funktionen mit komplexen Nullstellen, der quadra- 
tischen und der kubischen. 327 
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KLEINE MITTEILUNGEN 


Die Differentialgleichung der Biegungs- 
schwingungen eines axial belasteten 
geraden Stabes, dessen Material dem 
Hookeschen Gesetz nicht folgt. 1. Ein pris- 
matischer homogener Stab schwinge in der verti- 
kalen Symmetrieebene unter der Wirkung einer 
axialen Kraft P, die als positiv zu betrachten ist, 
wenn sie auf den Stab Zug ausübt. Wir beziehen 
den Stab auf ein rechtwinkliges Achsenkreuz, 
dessen -Achse mit der Stabachse in der span- 
nungslosen Anfangslage zusammenfällt und dessen 
y-Achse ihre positive Richtung nach oben hat. 


YJ 











2. Wir betrachten das Gleichgewicht an einem 
Stabelement von der Länge dz, das durch zwei 
Querschnitte aus dem Stabe herausgetrennt wird. 
Die Bewegung dieses Elements kann als eine 
resultierende Bewegung von zwei einfachen Be- 
wegungen aufgefaßt werden, und zwar von einer 
Translationsbewegung in der vertikalen Richtung 
und von einer Drehbewegung um eine senkrecht 
zur Biegungsebene gerichtete Achse. Diese beiden 
Bewegungen werden anfangs gesondert für sich 
behandelt. 

3. An den Endquerschnitten des Stabelementes 


N 
wirken die Schubkräfte —Q und Q + Star. Wird 


die vertikale Komponente der etwa vorhandenen 
äußeren Kraft für die Längeneinheit mit Y bezeich- 
net, so ist der auf das Stabelement fallende Anteil 
dieser Kraft Ydx und die Resultierende sämtlicher 
Kräfte 


dz+Yaz. 


Diese ist mit dem d’ 
widerstand 


Alembertschen Trägheits- 


”yY 
dt? 
des Elementes gleichzusetzen, wo oe die Dichte des 
Stabes und F die Fläche des Querschnittes bedeu- 
ten. Man erhält somit nach Wegheben von dx als 
Differentialgleichung der reinen Translationsbewe- 
gung des Stabelementes 


oF dr 


‚”y_2Q 
dt? =), 


4. Um zur Differentialgleichung der Drehbewe- 
gung des Stabelementes zu gelangen, betrachten 
wir das Gleichgewicht der statischen Momente aller 
an dem Element angreifenden Kräfte in bezug auf 
den Schwerpunkt des Elementes. Die Biegungs- 
momente werden als positiv betrachtet, wenn sie 
in der Richtung von der positiven z-Achse nach 
der positiven y-Achse hindrehen. 

Die an den Endquerschnitten des Stabelementes 


wirkenden Schubkräfte — Q@ und Q +3 
c 


fern das Moment 


(1). 


pP lie 
A 


E02 re, 
Das Moment der entgegengesetzt gerichteten Längs- 
kräfte P ist 
HE) 
u. 7  ° 
dr 
Elementes angreifenden 
oM 
2 dz)der Normalspan- 


Enden des 
und IM- 
nungen ergeben das re en Moment 


aM __ 
3298. PR 


Die an den 


Momente M 


5. Wir nehmen an, daß die der Dehnung & des 
Materials entsprechende Normalspannung o durch 
eine stetige und eindeutige Formänderungsfunktion 


WET rn a ae 


bestimmt ist, ferner daß das Material diesem Ge- 
setz uneingeschränkt folgt. Unter der Ber- 
noullischen Annahme, daß die Ebenen, die vor 
der Verformung senkrecht zur Achse des Stabes 
liegen, auch nach der Verformung eben verbleiben, 
läßt sich für das Moment M der Normalspannungen 
o folgende Reihenentwicklung ableiten t) 

M = C‚@-7C a —c 


sat.» . . . 16) 


(7) 
und 
== 10) 


Ca 


= 


Hier bedeutet e, die Dehnung im Schwerpunkt 
des Stabquerschnittes, o den Krümmungsradius 
der elastischen Achse des Stabes, I das Trägheits- 
moment des Querschnittes in bezug auf eine zur 
Biegungsebene senkrechte, durch den Schwer- 
punkt des Stabquerschnittes gehende Achse und 
,=/mP®dAF, ,=|MdF... auf dieselbe Achse 
bezogene höhere Momente der Querschnittsfläche. 

> 


y ist die der mittleren Längsspannung o,,= F ent- 


sprechende Dehnung, die durch die Formänderungs- 
funktion (5) als o„,=0(y) bestimmt ist. 0’ (y) 
o’ n ...sind die Ableitungen von (5) nach ge, wenn 
Wir beschränken uns bei unserer Betrachtung auf 
die an der Fließgrenze des Materials oder unter- 
halb derselben vorhandenen Spannungen. Dann ist 
der Wert der Dehnung e, für alle Materialien so 
klein. daß. man ihn mit Eins verglichen ohne 
erheblichen Fehler weglassen und (7) in Form 


1 


0 


schreiben kann, wodurch (6) in 


M 


übergeht. 


1) Arvo Ylinen: Erweiterung der Bernoullischen 
Biegungstheorie auf den unelastischen Bereich. Annales 
Acad. Sceient. Fenn. Ser. A. Tom. LVII, Nr.7 (1941), S. 10. 
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6. Nach den Lehren der analytischen Geometrie 
eilt für die Krümmung der elastischen Achse des 
Stabes die Formel 

0? y 
1 02? 


0 Bi d Y 213], * 
f r 5 2) 


Wenn so kleine Ausbiegungen des Stabes voraus- 
daß 
vernachlässigt werden darf, so ist 
U .; 
oe 292° 


2 
Ö2 PRIEE u 
gesetzt werden, >.) der Einheit gegenüber 
c 


Wird dieser Wert in 
eingesetzt, so erhält man 


”y.,.,[d) 
M = c - es - 
e 'd2? vr. Or? 


die Reihenentwicklung 


e e y\ | 
s\y.r: T 


womit (4) in 


| (10) 
az] 


der Momente (2). (3) 


Ö a r ge, z (* 41 
dr| "dr? ®\d 2) »\d 22 


übergeht. Die Resultierende 
und (10) ist 


oY 

BR .. 

Qdx Ps „det “ 
| 5 

d ”r (N, ,_ a 

3) a5) - a de) 


Sie ist der kinetischen Reaktion der Verdrehung 
des Elementes gleichzusetzen, die sich als Produkt 
der Winkelbeschleunigung und des Massenträg- 
heitsmomentes des Elementes dr in bezug auf die 
Rotationsachse ergibt. Der Winkel, welchen die 
Tangente der elastischen Achse des Stabes mit 


i Rena ; | 
der x-Achse bildet, kann gleich -— gesetzt werden, 
Ö: 
wodurch sich als die Winkelbeschlennigung 


0? er es dv y 
at?\dz/) Mr dx 


ergibt. Das Massenträgheitsmoment des Stab- 
elementes ist o/dr, da die Dieke dx des Elementes 
im Vergleich mit der Stabhöhe als unendlich klein 
angenommen wird. Man erhält somit nach Weg- 
heben von dx als Differentialgleichung der Dreh- 
bewegung des Elementes 


03 y 0y, 0 
_ — 0 — > 
older ? Past) | 


„”y ”y\’, „(d 2) - | 
DE ol) tale 


Eliminiert man aus den Gleichungen (1) und (12) 
die Schubkraft Q, so ergibt sich als Differen- 
tialgleichung der Biegungsschwin- 
gungen eines axial belasteten ge- 
raden Stabes, dessen Material dem 
Formänderungsgesetz (5) folgt 


(12). 


Pr ey Ö 
dt? Piet;e 


0? 3 (92 y\? 2 y\® 
BC 
dr? "\dr?) dr? 


oF 





IR or y 
3 el mt .d 


Hierbei haben die Koeffizienten c die Werte (8). 


P=0, so verschwindet 
linken Seite und in den 
0 einzusetzen. 


Wenn die Längskraft 
das zweite Glied der 
Koeffizienten (8) ist y 

7. Tritt an Stelle der allgemeinen Formänderungs- 
funktion (5) das Hookesche Gesetz o—Ee, so 
ist 0’ (y)=E und die höheren Ableitungen 0” (y) 
—0",(y)= 0. Hieraus folgt, daß c, = IE, die 
Koeffizienten cz — C3 O und (13) sich auf die 
bekannte Differentialgleichung 


‚ay 


‚”yY P 0? y 
dr 


] dry 
[0] 
dt? dr? 


Trug, 
9 mt 
der Biegungsschwingungen eines axial belasteten 
eeraden Stabes reduziert. 

8. Das obige Verfahren zur Ableitung der Diffe- 
rentialgleichung der Biegungsschwingungen kann 
weiter für den Fall verallgemeinert werden, daß 
der Stab einen veränderlichen Querschnitt hat. 
Die Querschnittsfläche F des Stabes, die Momente /, 
I;. I,..., die der mittleren Längsspannung o,, 

> 
die Koeffizi- 


= — entsprechende Dehnung y und 


F 

enten € in (8) sind dann als Funktionen von 7 zu 
betrachten. Wir bezeichnen diese neuen Variabeln 
mit F (x), I(z), I (2)... r(2), cı(2), & (2)... 
Nimmt man dieselben Rechenoperationen wie im 
vorigen vor, so erhält man als Differentialgleichung 
eines axial belasteten zeraden Stabes mit ver- 
änderlichem Querschnitt 


D) 
dr? 


3? y\? 82 y\3 
& ca)? 2) +c;(&) ( %) 
2 s or? “ 0x? 


d° y iy 
( 2) - } 


Helsinki. 


0° 2 2 
Y_p®yı 
C 0x? 


| j( 04). 


Arvo Ylinen. 343 


Eine einfache geometrisch-mechanische 
Analogie. Betrachten wir die Aufgabe: 

Auf welche Art und Weise sollen » un- 
gleichen Massen Mm, Ms, mM; .. m, auf 
einem Kreis verteilt werden, so daß 
das MassenzentrumdiesesdemKreise 
3oeinzeschriebenen n-Eckes, mitin- 
homogenen materiellen Eceken, mit 
dem Kreiszentrum zusammenfalle? 

1. Hat man gleiche Massen, muß n offenbar minde- 
stens 2 sein — diametralentgegengesetzte Punkte; ist 
n gerade, so sind je zwei diametralentgegengesetzte 
Lagen stets zulässige Anordnungen. Der Bedin- 
gung entsprechen ebenfalls die regelmäßigen Viel- 
ecke: bei geradem n ist dies offensichtlich, da man 
wiederum zu diametralentgegengesetzten Lagen der 
Punkte kommt; bei ungeradem n folgt die Bestä- 
tigung aus dem bekannten Satz, daß die n Wurzeln 
der Einheit geometrisch durch die Ecken der regel- 
mäßigen, dem Einheitskreise eingeschriebenen, Viel- 
ecken repräsentiert werden. Sind die Ecken durch 
21, 22, 23, - „, 2, bestimmt, sit a+2+23-+..+ 
2,„—0, da die Summe der Wurzeln der Gleichung 
2" —1—0 Null ist. Immer noch bei Voraussetzung 
gleicher Massen ergeben sich neue, außer den schon 
benannten mögliche Anordnungen der n Punkte, 
wenn wir kombinieren: z. B. lassen sich 5 Punkte 
noch auf unendlich viele Weise verteilen, wenn 
van sich das Fünfeck aus einem gleichseitigen 
Dreieck und einem beliebigen Durchmesser des 
Kreises gebildet denkt; analog. das Sechseck, ent- 
weder aus Quadrat und Durchmesser, oder aber aus 
zwei gleichseitigen Dreiecken; das Siebeneck, aus 
Quadrat und gleichseitigem Dreieck oder aus einem 
schon gewonnenen allgemeinen Fünfeck, kom- 
biniert mit einem Durchmesser usf. 








14). 


343 
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uf 
aß 
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Wir stellen uns nun die Frage: welches sind 
dieBedingungenfürdieallgemeinsten 
Anordnungen dern (gleichen oder un- 
eleichen) Massenpunkte, die unserer 
Forderungentsprechen,undwiekann 
man diese allgemeinsten Vieleceken 
erhalten? 

2. Die analytischen Bedingungen der Aufgabe 
lauten: 


n " 
mi u, =I, my a—=6, . . ID, 

ki kl 
i+y4=1, (k=1233,..,mM. . . (2). 


Dieses System (1,2) von n + 2 Gleichungen mit 2 
Unbekannten :,., Y; (k=1, 2,.., rn) wird bestimmt 
nur wenn 2 —3 ist. In der Tat wird der eine 
Punkt beliebig auf dem Kreise gewählt. so besteht 
dieses System aus 4 Gl. mit 4 Unbekannten. Im 
allgemeinen haben wir mit n+1 Gl. mit 2 (nr —1) 
Unbekannten zu tun. Bei n—4 bekommt man 
5 Gl. mit 6 Unbekannten, so daß man oo! Vielecke 
hat. Wir verfügen nun über einen Freiheitsgrad. 
so daß wir noch einen Punkt beliebig wählen. Bei 
n—5 erhält man 6 Gl. mit 8 Unbekannten, d. h. 
5° Fünfecke. Wird ein zweiter Punkt fixiert, so 
wird das System auf 5 Gl. mit 6 Unbekannten redu 
ziert: wir verfügen also über noch einen Freiheits- 
erad, denn man gelangt erst nach dem Auswählen 
eines dritten beliebigen Punktes auf ein System 
von 4 Gl. mit 4 Unbekannten. Im allgemeinen, da 
durch die Festlegung jedes weiteren Punktes die 
Gleichungen um 1 und die Unbekannten um 2 ab- 
nehmen, sieht man, daß, damit man, für n>3, auf 
ein bestimmtes Gleichungssystem kommt, man n — 2 
Punkte fixieren muß. Und jetzt wird gefragt, ob 
liese n—2 Punkte ganz willkürlich auf dem 
Kreise &enommen werden können? Die Antwort 
auf diese Frage lautet verneinend. 
Setzen wir 
n—2 n— 2 
2 r Y r € 

mn=X, Zmumn=Y.. . (8; 

1 1 
in diesen Gleichungen sind 27,47; (k = 1,2,..,n — 2) 
bekannt, da sie sich auf die ersten n —2 beliebig 
genommenen Punkte M,. Ma, . .. M„_. beziehen. 
Auf Grund von (1) lassen sich die Koordinaten der 
letzten zwei (der gesuchten) Punkte M,„_ı, M, aus 


Mn nt Mn In=X, Mn-ı%n ıtMmuYn=Y (4) 
bestimmen: 


X — Mn-ı Tn-ı Y- mMn-ı Yn-ı 


„= DE (5); 
in My My 
X — My In Y— mıY . 
In-ı= . Yn-ı > . . (6) 
Ma_ı My -ı 


(5), (6) sollen den Gl. (2) genügen, woraus sich die 
Gleichungen 


Myu-ı My-ı 


El A er 

In )- (& = = (8) 

| 79 Ba In My My 
ergeben. Diese Gleichungen stellen zwei Kreise 
K,_ı, K„ dar, und zwar die, daß auf dem Kreis 
K„_, dessen Zentrum S,„_,ı die Koordinaten 

X } i m 

—— , hat und dessen Radius r„_, = “ 
M-ı Ma-ı Myn-1 
ist, der Punkt Man-ı (£n-ı, Yn-ı) liegt und auf 


dem Kreis X, mit Zentrum S,, ; ) und Radius 


Mn My 









Ma-ı 
„=, der 


andere Punkt My, (Zn, Yn)- 
My 


Andererseits müssen, der Forderung der Aufgabe 
nach, diese beiden Punkte M„_,. M,„ gleichzeitig 
auf dem gegebenen Kreise K liegen, was bedeutet, 
daß sie die Schnittpunkte von K,„_,, K„ mit 
K sind. Die gesuchte Bedingung für die Anord- 
nungen der n Punkte auf dem Kreise K wird also 
durch die Bedingung des Schneidens der Kreise (7), 
8) mit dem gegebenen Kreis gegeben; sie lautet: 
Mn |. Br 52210, 
ET 


My-ı Myu-ı My-ı 
wo U=-+YX?’+Y? der Abstand des Punktes 
RB, AE 7) \ und was dasselbe ist des Punktes 


n—?2 n—2 
a Y r . 
Ru-.| 2) mı 2, me ur) vom Koordinaten- 
u 4 


anfang O ist. Die obigen Ungleichheiten lassen 
sich einfacher folgendermaßen schreiben: 


Ma-ı — Mn| SU Mt Mu: - (9). 


3. Wir wollen zunächst den Fall gleicher Massen 
(= 1) betrachten. (9) gibt einfach | L 2; das heißt, 


\ > 
der Radiusvektor R),_s (X, Y), der der Resultante 


> 

der ersten n-- 2 Vektoren OM,. gleich ist, darf nicht 
größer als das Doppelte des Radius des gegebenen 
Kreises sein. Mit anderen Worten: das Ende des 


> 

Vektors R,_> muß innerhalb des Kreises € liegen, 
der zum Kreise K konzentrisch ist und den Radius 
2 hat. Die Gl. (7), (8) reduzieren sich auf den 
einzieen Kreis 


ld Bu ; L > Be 1 


welcher K in den gesuchten Punkten M„_,. M, 
schneidet. Damit die Aufgabe eine Lösung hat, 
darf das Zentrum von (10) höchstens auf dem Kreise 
€ liegen; in diesem Fall berührt er den gegebenen 
Kreis K. Dann besitzt das gesuchte Vieleck eine 
„Doppelecke“. Aus Bild 1 ersieht man die Konstruk- 
tion eines allgemeinen Fünfecks des erwünschten 
Typus. Bild 2 stelltden Fall unmöglicher Lösung dar. 

4. Ungleiche Massen. Diese können auch negativ 
sein. Jetzt schließen wir aus (9), daß, wenn man 
um den Anfangspunkt O zwei Kreise C’, C” mit den 
Radien |m„_, — m„| und |m„_ı + m,| beschreibt, 

> 


das Ende des Radiusvektors R,_s, im Falle mög- 
licher Lösung, nicht außerhalb des Ringes der bei- 
den konzentrischen Kreise C’, C” liegen darf. Die 
Sehnittpunkte von K mit (7) und (8) geben uns zwei 
richtige Lösungen, wenn die Punkte der Punkt- 
paare M„_-1,Mu: Mh-ı, Mn auf verschiedenen Sei- 
ten der Geraden liegen. Bild 3 zeigt ein Fünfeck 
dieser Art. Die Konstruktion folgt aus der obigen 
analytischen Betrachtung der Kreise (7), (8), aus der 
man sieht, wie die Mittelpunkte S,_,, S„ und die 
Radien r„_,, "„ gefunden sind (Radius von K—1). 
Im betrachteten Falle sind m,, Me, ms, ma, M, von 
vornherein gegeben und M,, M», M, beliebig aus 
dem Kreise K genommen, jedoch so, daß, sobald wir 
eine Lösung haben, und zwar ohne Doppelpunkte, 
die Bedingung (9) (wobei nur an einer Stelle ein 
Gleichheitszeichen zulässig ist) befriedigt wird. 
Kommen an den beiden Seiten Gleichheitszeichen, 
erhalten wir Bild 4, in welcher sich das Fünfeck auf 
ein Viereck reduziert hat; die vierte Ecke trägt 
jetzt zwei verschiedene Maßen. Keine Lösung er- 
eibt sich, wenn K„_,, K„ den Kreis X nicht treffen. 

Haben wir schließlich X=Y=0, bilden also die 
gewählten n—2 Punkte schon ein (na—2)-Eck der 
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Bild I. 


Bild 2. 
betrachteten Eigenschaft, so hat unsere Aufgabe, 
beim Falle ungleicher Massen, keine Lösung. In 
der Tat werden die Kreise (7), (8) 


> m 
&u-ı tNM-ı= ( . 
My, 


(11), 


.. 128) 


My 1 s 
My ) 


Be 
In tr Ya | 


Das sind zwei konzentrische Kreise um O, die K 
nicht schneiden können. Nur wenn die letzten zwei 
Massen gleich sind (die ersten n— 2 brauchen ja nicht 
gleich zu sein), läßt die Aufgabe eine Lösung zu; 
die beliebige diametralgegengesetzte Lage dieser 
zwei Punkte ergänzt die Konfiguration zur ge- 
wünschten Lösung. 

Damit sind alle möglichen allgemeinsten Lösungen 
des Problems erschöpft. 

Diese geometrische Aufgabe deckt sich mit fol- 
zendem mechanischen Problem: 

Auf welche Art und Weise sollen n 
verschiedene (schwere oder gewicht- 
lose Massen auf der Peripherie einer 
homogenen Scheibe, die auch als mas- 
selosangenommen werden kann, ver- 
teilt werden, so daß bei der Drehung 
dieser Scheibe um eine vertikale zur 
Scheibe senkrechten Rotationsachse, 
durch ihr Zentrum, die Drehachse 
keinen Druck erfährt? (Bild D). 

Nimmt man an, die Massen seien gewichtlos, so 
ist es gleichgültig, ob die Drehachse vertikal oder 


Bild 3. 


Bild 5. 


nicht ist; das Nichtvorhandensein der Gravitation 
verlangt ebenfalls nicht „Scheibe ohne Masse“, 
sondern sie muß nur homogen sein. In jedem Fall, 
betrachten wir die Scheibe und die n Massenpunkte 
als einen festen Körper, und zwar als einen solchen 
mit unstetig konzentrierten Massen in konstanten 
Zwischenabständen. Unter diesen Voraussetzungen 
ist es offenbar, daß auf jeden der n Punkte eine 
zentrifugale (inertiale) Kraft wirkt, welche durch O 
geht uud die Größe m; © R hat (w == Rotations- 
geschwindigkeit, R — Scheibenradius). Ist P der- 
jenige Scheibenpunkt, der mit dem Schwerpunkt 
(Massenpunkte) der n materiellen Punkte zusam- 
menfällt, so wird auf P die resultierende Zentri- 
fugalkraft — Mw?r wirken, wo M = Zm,, r—= OP 
ist. Diese Druckkraft würde die Drehachse der 
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Scheibe erfahren infolge der angenommenen festen 
Bindung der n Punkte mit Scheibe und Achse. 
Damit dieser Druck gleich Null sei. muß der 
Punkt P mit O zusammenfallen. Mechanisch drückt 
man diese Tatsache so aus: Die Rotationsachse 
bleibt frei, d. h. sie erfährt keine Druckkomponente, 
außer einer solchen in ihrer Richtung, wenn sie 
eine Hauptträgheitsachse in bezug auf O ist. Sie 
muß infolgedessen durch O gehen, da sie, wie ange- 
nommen, normal der Scheibe steht. Setzt man der 
Kürze halber ®’R — 1, was bedeutet, daß wir die 
Zentrifugalkräfte? gleich den entsprechenden Massen 
reduzieren, so wird es sich um eine derartige Ver- 
teilung der Massenpunkte M,. auf dem Scheibenrande 


> 
handeln, daß wir für die Summe der rn Vektoren m; 


> : 
mit Richtung OM,, Null erhalten. Dies ist aber iden- 
tisch mit der Bedingung (1.2). Man zerlegt im 
Falle gleicher Massen (Bild 1 und Bild 2) den Vektor 
> n—2 
-ORy = — 3 
> 1 
sein muß. da m. _, + m; =Z2, so in zwei Vektoren, 
daß ihre Enden auf den Kreis ÄK fallen. Bei unglei- 
> 
Bild 3 und 4) ist R’„_» in zwei 
> > 
von vornherein gerebene Vektoren m,_,. m 
zerlegen, was nur möglich ist. wenn die 
> > > 
Vektoren R’„_s. M,_ı. m,„ ein richtiges Dreieck 
bilden. Die betrachtete Aufgabe stellt folglich eine 
Veralleemeinerung des bekannten Stevinsprin- 
zips dar. 

Man sieht aus dieser doppelten Behandlung der 
Aufgabe, daß bei gewissen Problemen, wie bei dem 
untersuchten, das vereinigte, analytisch-vektorielle 
Verfahren sehr instruktiv sein kann. 

Budapest. Bl. Dolaptschiew. 355 


> 


> _>_, 
2,108, ,=- mj,), der 


chen Vektoren 


„. Zu 
drei 


Zur Berechnung der Eigenfrequenzen 
inhomogener Saiten. In einer in dieser Zeit- 
schrift erschienenen Arbeit!) behandelt A.Erdelyi 
das Problem der Annäherung der Eigenwerte der 
Differentialgleichung 


yY"+ig()y=V0 \ MN 
(Randbedingungen 4 (0)=y(1)=0) f ; 


wo 4 (2) eine im Intervalle 0 <r<< 1 postive, stetige 
und zweimal stetig differentiierbare Funktion be- 
deutet. Er wendet die Weinsteinsche Methode 
auf die Näherungsfunktionen 

x 
a ‚Vaw)du 
n=g +sinnn- 
\yaelı)du 
6 
Resultat ist die Ab- 


der Lösungen y,„ an. Sein 


schätzung 
N y 
1 i + 1 
\Ya(wWdu 
\o 


! i 
\rgmdz I | 


\ygawW)du' 
v 


wo S, eine berechenbare Konstante bedeutet und 


der Kürze halber 
bg’? 
16 g? 
gesetzt worden ist. 


1) Eigenfrequenzen inhemogener Saiten. Z. 


Math. Mech. Bd. 18 (1938), S. 177—185. 


angew. 


Anschließend daran beschäftirt er sich mit der 
Frage, in welchem Falle S,—0 ist, d. h. bei 
welcher Beschaffenheit der Funktion 
q(z) die Eigenwerte 4, der Differen- 
tialgleichung (l) mittels der Formel 
(2) exakt berechenbar sind? 

Er leitet ab. daß in diesem Falle r 
also q (x) die Differentialgleichung 


const ist, 


’ 


4gg' Ed 2 23) „©; 
befriedigen muß. Die 
durch die Gleichung 


. na ; 
Sams f FR 
\Yalu)du 


ausgedrückt. Für C=0 findet er die Lösungen 
dieser Gleichung, nämlich = (ax + b)”* und be- 
hauptet, daß im Falle C=F0 Lösungen von der Form 
q= (ax + J* existieren. Die allgemeine Inte- 
gration der Gleichung (3) ist ihm aber nicht ge- 
lungen. 

Der Gegenstand dieser Note ist die Ableitung 
aller partikulären Lösungen der Dif- 
ferentialgleichung (3). Zu diesem Zwecke 
machen wir die Substitution 


Eigenwerte werden dann 


0=3 


die die Gl. (3) in die Gestalt 


LT a I 6 u 


überführt. Wenn wir beide Seiten mit s’ multi- 


plizieren, erhalten wir 

s’ 

168’S"—=C —. 
5 


Diese Gleichung läßt sich schon einmal integrieren: 


2 
is 


m 


8 (s’)’ = 


Nach Multiplikation mit s? bekommen wir 


Wenn wir wiederum eine Veränderliche 


q - 1, ” p ++ ” * 
t—s’—q ” einführen, so wird unsere Gleichung 
von der Gestalt 


neue 


( 
-5+4t 


2 (it) = 


sein: in dieser Form läßt sie sich aber durch Tren- 
nen der Veränderlichen integrieren. Ihr Integral ist 
für AF0 


— (4), 
rt 


bzw. für 


==, - 
(£+B) 


worin B eine neue Integrationskonstante ist. 
Aus Gl. (4) erhalten wir Erdelyi’s Lösungen 
für den Fall C=0, wenn wir in ihr C gleich 0 


setzen, während Gl. (5) die von Erde&lyi gefun- 
denen Lösungen für den Fall € ==0 lieıert. 
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Es sei noch bemerkt, daß singuläre Lösungen 
unserer Differentialgleichung nicht existieren. Denn 
die notwendige Bedingung dafür, daß die im Koor- 
dinatensystem z, y, 2 durch die Gleichungen 


=r 
=q(z2, 4, B) 
= (2;4,B) 


definierten Kurven eine fokale Fläche mit endlichen 
y-Werten besäßen, ist nicht erfüllt. 
Budapest. E. Makai. 370 


Annäherung der logarithmischen Spi- 
rale durch Kreisbogen. Die logarithmische 
Spirale gilt in der Technik als eine der wichtigsten 
Kurven. Man hinterdreht nach ihr die Zähne von 
Fräsern und krümmt nach ihr bogenförmig geführte 
Blätter von Scheren. Fließfiguren bei statischen 
Festigkeitsversuchen von Rohren zeigen die Form 
einer logarithmischen Spirale. Weiter finden wir 
sie theoretisch wie praktisch bei vielen Bildern von 
Strömungsvorgängen z. B. Stromlinien, Leitblechen 
in Turbinen und Kreiselpumpen, Exponential- (Kurt 
von Sanden) Überfällen und Spritzrohren am Rund- 
herd von Erzaufbereitungen. 

Die Gleichung der logarithmischen 
Polarkoordinaten r, 9 lautet 


Spirale in 


FERN SE 5 (3), 
darin ist 7 der Winkel, unter dem die Spirale alle 
vom Mittelpunkt © ausgehenden Halbstrahlen 
schneidet, gemessen zwischen Halbstrahl und der 
positiven Fortschreitungsrichtung der Kurve. 

Im folgenden wird eine Näherungskonstruktion 
für die logarithmische Spirale gebracht, die rech- 
nungsmäßig wie konstruktiv in ihrer kreisbogen- 
förmigen Kurvenausführung t) viel vorteilhafter als 
1) ist und dabei beliebig genaue Resultate ge- 
stattet. 


Bild 1. 


Die Strecke OP, = a und der Tangentenwinkel r 
sind gegeben (s. Bild 1). Um den unter « verlaufen- 
den Halbstrahl r zu berechnen, ist zunächst nach 
dem Peripheriewinkelsatz 

ION P,= ET -0=y. (2) 
und 

{POP = 


IT P,M. Pı =na—-yr-o=a. 0). 


Nach dem Tangentenwinkelsatz ist XP,P,A im 
Punkt P, halb so groß wie der mit ihm auf dem 
1) Ein zeichnerisches Näherungsverfahren für die log- 


arithm. Spirale ist in der „Deutschen Mathematik“ Bd. 3 
(1938) S. 269 bis 272 gegeben. 


stehende Zentriwinkel P,M,P:. 


und {AP,O 


Bogen P,Pı 


i a 
daher ist =; 


ı=r1+7 A 3 


Mit Gleichung (4) sind alle Winkel im Dreieck 
OP,P, bestimmt und man kann jeden Fahrstrahl 
der logarithmischen Spirale beliebig genau be- 
rechnen, wenn man den Durchgang der Spirale im 
Punkt P, als ein Kreisbogenstück mit dem Mittel- 
punkt in M, ansieht. Ist ein unter <X« im Inter- 
vall 0—90° verlaufender Halbstrahl r einer log- 
arithmischen Spirale (aund rfest) zu berechnen, so 
kann man statt den Halbstrahl nach (1) zu be- 
rechnen, nach dem Sinussatz angewendet auf Drei- 
eck OP, P, Setzen 


sin ® 
sin y 


". 
s+- 


BB 5 r 
da nun aber und =t—-- ist, wird 


sint—+t 
va SEE © 


. a 
sint—> 
Gleichung (5) liefert für (1) nur grob zenäherte 
Resultate, zur besseren Annäherung wird (5) wie- 
derholt angewendet. Es wird dann 
E a ü ’ a % 
sint— sinT ” 
i 2n * 2n i 
lim — (6). 
n>& a 


e {gra- 


sind — 5, sint— 5, 

Darin werden r und a in Graden addiert und dann 
der Quotient der Sinuswerte gebildet. Je nach der 
verlangten Genauigkeit wird ein entsprechend 
großes n genommen. Auf diese Weise kann man 
nach (6) sehr bequem rechnen. In Zahlentafel 1 
sind die Abweichungen für die in Frage kommen- 
den Werte n angegeben. 


Zahlentafel 1. 





Abweichungen nach (6) in 0/9 


T =] =? e—s n"—1W n — Öl 








< 1,010,0 .. 10,2)0,0 .. 2,20| 0,0 ... 0,30 0,0 ... 0,074) 0,0 ... 0,00 





” a ’ i 
Wirda>,, so muß unterteilt und allgemein 


gesetzt werden 


nXx 


E > «a 
S sint+ Fr \ 
a(ectgre) za u 
i sin ı 


OP, 


2n 


zeichnerische Konstruktion der logarith- 
mischen Spirale lassen wir Halbstrahlen unter 
O, a, 2a,... verlaufen und berechnen darauf 
nach (6) ihre Länge oder tragen nach (4) alle 
Dreieckswinkel zeichnerisch an, wodurch die Kur- 
venpunkte Po Pi Paz... entstehen. Die Spirale 
soll nun zwischen P, und P,, P, undP,... stück- 
weise durch Kreisbögen ersetzt und somit be- 
quemer als bei punktweise durchgeführter Berech- 
nung mit dem Kurvenlineal gezeichnet werden. Da 
alle Dreiecke OP,P,, OP,Ps.... einander ähn- 
lich sind, braucht man nur die Angabe der Kon- 
struktion des Bogens P, P, inachen. 


Für die 
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Der dem Dreieck O P,P, umbeschriebene Kreis 
schneidet die auf P,P, nach der Seite von O hin 
errichtete Mittelsenkrechte in M,. Der Kreis um 
V, mit M,P,=M,P, als Radius wird als Näherung 
ler Spirale zwischen den Punkten P, und P, be 
nutzt. Für die Konstruktion beachtet man, daß die 
Kreisbögen P,P,. P,P. ... mit stetiger Tangente 
aneinander schließen sollen. Da ar, PP, 
ist XP, MıPı ıP,OP,. Der Mittelpunkt M, 
liegt also auf dem Kreis über P,P, als Sehne mit 
lem Peripheriewinkel «, 

Genauigekeitsnachweis: Bei der Be 
rechnung der logarithmischen Spirale nach Glei 
chung (1) und der hier neu abgeleiteten Gleichung 
6) erlangt man eine beliebig Genauiekeit 
durch Benutzung eines entsprechend hohen Ex 
ponenten. 

Eine hohe Genauigkeit in der zeichne- 
rischen Konstruktion erhält man, wenn die Drei 
ecke nach den errechneten Halbstrahlen gebildet 
werden. 

Wird die Kurve konstruktiv nach den Winkeln 
}) entwiekelt. so schneidet die Näherungskurve ( 
die Halbstrahlen OP, OP, . nieht unter dem 
Winkel r wie die Spirale S, sondern unter 2 @—-P, 
mithin ist 


hohe 


ebenso 


8) 


ein Maß für die Abweichungen dieses Bogenstückes. 
Die Größe «dieser Abweichung ist abhängig von 
wird der Fehler 


R 1 
dem gewählten a. Für « a8 


bis 2° eroß. Eine solche Genauigkeit genügt 
meistens für die Praxis. Wenn man «a noch kleiner 
wählt. so verringert sich die Abweichung zu- 
sehends. Bei dieser Konstruktion schneidet € nun 
S unter dem Winkel 4. Für 9 O liegt S deshalb 
von O aus gesehen außerhalb €, ebenso wegen der 
Ähnlichkeit für 9 @a+0. Daher muß für 
y =a—.( die Spirale von O aus gesehen inner- 
halb € liegen. Zwischen P, und P, müssen sich 
die Kurven also mindestens noch einmal schneiden. 


. } ah. 
Dieser Schnittpunkt liegt stets bei _ « und be- 


dinet dadurch eine sehr gute Annäherung. 


TO P, M, [9 


Beweis: Da +P, so 


folgt 
T a 


=+- 3% ET RR 


Bezeichnet R 
folet nach 


den Kreisradius M,P,=M,Pı. so 
(lem Sinussatz für Dreieck OM,P, 
a 


*OS 
os > 


R a 
sin (x 
und 


cosT 


e a\' 
sin (r ar 


Der Abstand des 
Punkte M, sei f. 


r? k? 


Spiralenpunktes r, 9 
Dann ist 


vom 


+ r? 


2krcos(t—- Y) . . 12). 


. 1 E ’ A 
Für g ‚ averschwindet das letzte Glied auf der 


rechten Seite. Es ist ferner nach (1) und (5) 
sin (? 4 


— für « 
A a\ ’ 
sın (ı 


woraus 


sin [7 


und mit (10) die Behauptung folgt. 
R ist auch ein mittlerer Wert der Krümmungs 
radien 
a 
sintT 


sin |T 


sin (t 


a ö 
> | sin I 


in den Punkten P, 
1 


>) 


und P,. und die Stellen 


q 0, a, «a sind im Intervall O g a die ein 


zigen, an denen S und € sich beim zeichnerischen 
Konstruktionsverfahren schneiden. 

Die bei der Berechnung von Leitblechen in Tur 
binen immer erforderliche Kurvenlänge b .Jer 
logarithmischen Spirale, die bekanntlich durch 
Integration bestimmt werden kann, wird nach (10 


Rn. 
16) 
180 

mit einer für die Praxis hinreichenden Genauigkeit. 
Die eingangs erwähnte und bisher in der Praxis 
schon bewährte Näherungskonstruktion durch Kreis 
bögen ist hier in eine sowohl in Rechnung wie 
Konstruktion beliebig genaue Kreisbogenkonstruk 
tion überzeführt worden. Durch den mit Formel (4 
gegebenen Fortfall der Zahlentafel ist die Be 
schränkung auf die dort angeführten Winkel auf 
gehoben. Diese können für jede Größe der logarith 
mischen Spirale ohne weiteres bestimmt werden. 
Vergleicht man die lögarithmische Spirale mit 
der hier abgeleiteten Näherung, so ergeben sieh der 
ersteren gegenüber beachtenswerte Vorteile: einmal 
kann ohne Logarithmen zerechnet werden, zum 
anderen ist bei logarithmischer Berechnung hier 
weniger Rechenarbeit erforderlich als nach (1). 
Hinzu kommt die nach (16) einfache näherungs 
weise Berechnung der Kurvenlänge und des In 

halts der logarithmischen Spirale. 


Berlin. H. Meincke. 342 


Eine Schieberanordnung für die Schlüssel- 
gleichung |; (7 («) + y (ß)) +h(h)=1; (a,B,y), 
wobei nur das Endresultat abzulesen ist. 
Die obige Schlüsselgleicehung I 
liebige „gesunde“ Funktionen 
neinen der Behandlung mit dem 
Widerstand entgegen, da sie 
/Zwischenresultats verlangt. 
keit auch Fehlerquelle 
beschriebene Anordnung von zwei Schiebern im 
rechten Winkel in Verbindung mit einer Kurven 
schablone vermieden. 
(a) ist auf der festen 
schiebers 1 angebracht, y 
Teil II (Zunge 1). Wird zur Abkürzung + y—u 
gesetzt. so ist die w-Skala mit der g-Skala zu- 
sammenfallend zu denken. f,(w) wird als Ordinate 
über u konstruiert. (Bei der praktischen Ausführung 
wird die Kurve als Blechschablone zeschnitten und 
an den festen Teil I des Schiebers 1 zeschraubt.) 
Mit dem Ende € der Zunge 1 (w-Skala) ist im 
rechten Winkel zu dieser der feste Teil IV des 
Schiebers 2 fest verbunden, so daß er sich bei Ver- 
schiebung der Zunge 1 (Bild 1) auf und ab bewegt. 


fl. $%, y be 
setzt im allge 
techenschieber 
Ablesung eines 
Unbequemlich 
wird durch die hier 


| )iese 


Skala (I) eines Rechen- 
(P) auf dem beweglichen 


12 
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Schieber 7 


I 
u-prw 
| 
| 
| 


4 (By. 


IV 





F 


EU 





Schieber 2 


Zunge 2 


Bild 1. Der besseren Übersicht wegen wurden die Skalen- 
träger im Bild mit größerem Abstand gezeichnet. In 
Wirkliehkeit schließen % und w Skala (in Bild I und II) 
ebenso /» und /3 Skala (im Bild III und IV) dieht anein 
ander (als Zunge und fester Teil von Schieber 1 bzw. ?), 
per punktierte u-Träger fällt mit I zusammen. 


sich die 
f,-Skala. 
und ge- 


\nf Zunge 2 des Schiebers 2 befindet 
f,„-Skala, auf dem festen Teil IV die 
Nötigenfalls können die Skalen eleich 
genläufig aufgetragen werden.) 

Die Rechnung geht so vor sich: $ (a) auf I auf 
suchen (im Beispiel=4APB), y(ß) auf II B’C) 
B auf I gegenüberstellen, dabei verschiebf sich der 
feste Teil IV des Schiebers 2 in die Lage C—F. 
\uf der gedachten u-Skala würde dann @ + y u 
abzulesen sein 1C), wo die Kurvenordinate 
D=[/,(u) ist. Jetzt wird die Zunge 2 einge 
schoben bis zur Berührung (Anschlag) mit der die 
Kurve f, darstellenden Blechschablone. fs (— DE) 
auf Zunge 2 gibt auf dem festen Teil IV des 
Schiebers 2 das Resultat /; (- CE). 

3ei der mechanischen Durchbildung des Systems 
könnte Zunge 2 zwangsläufig an der Kurve f, (u) 
entlang zeführt werden. wenn Zunge 9,1 eingestellt 
wird. Bei einer behelfsmäßigen Ausführung wiri 
man jedoch die ganze Anordnung auf Millimeter 
papier zeichnen, Zunge 1 und fester Teil des Schie 
werden durch einen rechten Winkel aus 
Millimeterpapier, Zunge 2 durch einen Streifen 
ebensolchen Papiers ersetzt. 


bers 2 


Von den zahlreichen Ausdrücken, die auf obige 
Schlüsselgleiehung zurückgeführt werden können, 
seien hier nur die Exponentialgesetze besonders 
erwähnt. 

Böblingen. H. 


Schminke. 356 


Bemerkung zum Buffonschen Nadel- 
problem. Bei diesem Problem handelt es sich 
bekanntlich um folgendes '); Man wirft auf eine 
ebene horizontale Platte, auf welcher Parallelen in 
eleichem Abstand «a gezogen sind, eine dünne 
Nadel von der Länge !<a und fragt nach der 
Wahrscheinlichkeit (dafür, daß hierbei die Nadel 
auf eine der Parallelen zu liegen kommt: vergleiche 
Bild 1. wo M«r,y) der Mittelpunkt der Nadel und 
d ihre Neigung geren (die .r-Achse ist 

a 
”) 


1) Vgl. z. B. E.Üzuber: Wahrscheinliehkeitsrechnung, 
Leipzig und Berlin 1938, 1. Bd. 5. Aufl., S. 98. 
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Bild 1. 


Das Problematische an dieser Aufgabe lieet nun 
in der Frage: wie soll man bzw. wie hat man die 
möglichen Nadellagen „abzuzählen”, oder anders 
gesprochen, welche Wahrscheinlichkeitsdichte ist 
zu nehmen? Offenbar geht in das Verhältnis der 
vünstigen zu den möglichen Fällen die y-Koordi 
nate nieht ein. Wählen wir dann mit Buffon als 
Parameter .r und #, und zählen gemäß 


\drd# 


ab, wählen also die Wahrscheinlichkeitsdiehte kon 
stant, so erhalten wir bekanntlich für (die gesuchte 
Wahrscheinlichkeit W 


So plausibel diese Abzählung auch erscheinen mag, 
so erscheint es zunächst nicht ausgeschlossen, ob 
es nicht auch andere, vielleicht ebenso plausible 
\bzählungen 


\ve(r.d)drdo) 


beziehungsweise, wenn wir vorübereehend die 


-Koortdinate mitnehmen 


ODE... TB 


gibt, wo ww (9) beziehungsweise w (7, 4,9) nicht 
konstant ist. Eine Begründung, warum also 4 
konstant, das heißt also gemäß [dxrd® bzw. 
(d.rdyd® abzuzählen ist, muß erbracht werden ?). 
Eine solehe Begründung scheint uns durch die 
folgende Bemerkung möglich, welche sinngemäß 
aueh für andere ähnliche Aufgaben über geome 
trische Wahrseheinlichkeiten gilt. Stellen wir näm 
lieh an das Maß (1) die natürliche Forderung ?). 
daß es für jeden Nadelbereich,. oder wie wir uns 
anschaulich ausdrücken wollen, für jeden „Nadel 
haufen“ 3 invariant gegenüber den Bewegungen 
der Ebene ist, so ist w notwendige konstant. 
Nehmen wir nämlieh durch eine solche Bewegung 
unseren Nadelhaufen ® mit, so entsteht aus ihm 

ein neuer Haufen B* vermöge der Zuordnung 
+ reosa+ ysin«a 


rsina—+yecosa 
a*, b*, a Konstante 


2), Vgl. hiermit auch die Note des Verf. „Bemerkungen 
zu Mises’ Behandlung des Nadelproblems von Buffon“, 
Deutsche Mathematik Bd. 6 (1941), Heft S. 108 bis 113, 
und speziell S. 112 Abs. 7, wo, um auf w = eonst zu schließen, 
die Zuflucht zu dem Prinzip vom mangelnden Grunde ge 
nommen werden mnßte. 


3) Hierüber und über die nun folgende kleine Rechnung 
siehe auch W. Blaschke: Vorlesungen über Integral 
geometrie, Erstes Heft. Leipzig und Berlin (195), 8. >. 
Übrigens ist die oben angegebene Abzählung vermöge 
dxdyd»% in der Sprache von Blaschke niehts anderes als 
die Dichte der ungerichteten Linieneleinente bzw. mit 2 
multipliziert die kinematische Dichte Poincares. 
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und unsere Invarianzforderung ergibt 


\wir,y.d)drdyd» 
a 


{ 


—= | w (r*, y*, 9*) dr* dy* d9*. 


Ak 
T 


[ransformiert man hier das rechte Integral veı 
möre (2) auf den alten Bereich ® zurück. und be 
ichtet, daß 


drdyd» dr*dy* d»* 











erhalten wir 
7,y.d9)drdydd 


w(a*.y*,d*)drdydd 


\ 
Ba 





und mithin, da ® beliebig gewählt werden kann. 


w(r,Yy,%) EEE 2 . Pr: 
Da man durch geeignete Wahl der Bewegung d. h. 
der Konstanten a*, b*, a jede Nadellage .r, y,® in 
eine feste .r*, y*, #* überführen kann, so folgt 
aus (3) in der Tat 


w (2.4, 0) const. 


Breslau. Ernst Mohr. 351. 


Bemerkung über ein von W. Höffding 
vorgeschlagenes, maßstabsinvariantes 
Korrelationsmaß. In Heft 6 (1941) S. 364-379 
dieser Zeitschrift habe ich das Korrelationsproblem 
der Statistik als Eigenwertproblem einer Fredholm 
schen Integralgleichung begründet und dabei auch 
lie Rolle der klassischen Korrelationsmaße, nämlich 
les Korrelationskoeffizienten r Gl. (5), der beiden 
Pearsonschen Korrelationsverhältnisse Kk?’,x bzw. 
k?xy@Gl. (7) und des Pearsonschen Kontingenzmaßes 
f? Gl. (8) im Rahmen dieser allgemeinen Theorie 
klargestellt. Nach Apschluß jener Arbeit wurden 
mir zwei neue, einschlägige Untersuchungen von 
Wassilij Höffding bekannt: 

1. Maßstabinvariante Korrelationstheorie. Schriften 
des Mathematischen Instituts und des Instituts für 
anrewandte Mathematik der Universität Berlin. 
Bd. 5. Heft 3. Berlin (1940) und als Fortsetzung 

2. Maßstabinvariante Korrelationsmaße für  dis- 
kontinuierliche Verteilungen, Archiv für mathema 
tische Wirtschafts- und Sozialforschung. Bd. 7. 
Heft 2 (1941). 

Es dürfte von Interesse sein, ebenso wie die ge 
nannten alten Korrelationsmaße auch die Höffling 
schen Vorschläge und Ergebnisse vom Standpunkt 
der Eigenwerttheorie aus zu betrachten. 

Das Hauptergebnis meiner Theorie ist die eng 
mit dem kleinsten Eigenwert einer Integralgleichung 
zusammenhängende Maximalkorrelation 


tz + 
K \ \ /@)gy@) wie.y)daedy 


u. ee) 


Maximum, 


die durch ihr Verschwinden (K -—- 0) mit Sicherheit 
anzeiet, wenn die beiden Veränderlichen z und y 
unabhängig voneinander sind, und die anderseits 
mit Sicherheit den Wert K—1 annimmt, wenn 
zwischen ihnen vollständige Abhängigkeit besteht. 
Entscheidend (und für die praktische Berechnung 











von K erschwerend) ist dabei, daß als Prüffunk 
tionen f(#) und 9(4) normierte Eigenfunktionen 
des kleinsten Eigenwerts der Integralgleichung zu 
verwenden sind. Gegenüber beliebiren Merkmal- 
transformationen X X (x), } Y (y) ist die Größ« 
k unempfindlich. 

Das Höffdinzsche maßstabsinvariante Korrela 
tionsmab, von ihm als normierter Korrelations 
koeffizient 2 bezeichnet, ist nun nichts anderes als 
ter Wert des obigen Doppelintegrals auf Grund der 
die Normierungs- und Orthogzonalisierunesbedin 
zungen (10) und (12) erfüllenden Prüffunktionen 







gay)=] 3(2\ WolY dy— 1). 


F. 


Für diese willkürlich gewählten Funktionen nimmt 
der obige Ausdruck im allgemeinen natürlich nicht 
seinen Größtwert an: daher ist 0? kK? 

Man kann den Vorgang bei Höffding weniger 
formal auch auf folgende Weise «deuten. Höffding 
führt eine Merkmaltransformation aus, derart. daß 
(lie beiden eindimensionalen Verteilungen w; (#) 
und ws(y) gleichförmige werden mit der Häufig 


keitsdichte Eins im Bereich zwischen ,„ und 


1 2 } 
+5. Wählen wir statt dessen, um besseren An 
schluß an die Gleichungen meines früheren Auf 


satzes zu gewinnen. das Intervall zwischen 1 
und + 1. so lautet die erforderliche Transformation 


” 
\ r) > \ m, r)« T l 
„U 
Y(y 2 \ w,(y)dy E. 
und dureh sie wird 
’ 1) 
w, (X): > | | N N | 
für | ‚sonst Null. 
1 ' 1 
,m-! | ISYs+i1 


Nach Ausführung dieser Transformation sind die 
Orthogonalpolynome (26) des Problems bis auf kon 
stante Normierungsfaktoren die Legendreschen Po- 
Iynome P,„ (2). Es ist fürn =0, 1, 2, 


YH(AM)=Y2n-+1P.(X), 
„an (Y)=Y?2n+1P,.(Y) 


Wir erhalten damit das methodisch bemerkens 
werte Ergebnis, daß man bei einer beliebigen geo- 
metrischen Verteilung w(z,y) durch Anwendung 
dieses Kunstgriffs stets mit einem einzigen be- 
stimmten Funktionensystem, z. B. dem der Legen- 
dreschen Kugelfunktionen auskommt. Bei arith 
metischen Verteiluneen allerdings zelingt diese 
Zurückführung auf ein einheitliches Funktionen- 
system nicht. 

Auf Grund der Transformation tritt an die Stelle 
der Verteilung w (x,y) eine neue =(X,Y), für die 


e(r,yJ)drdy=w(XY)drdY 


ist. Werden für das transformierte Variations 
problem 


12° 
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+1 
\ KM g(Y)w(X,Y)dXdY 
N 


K= — Maximum 


als einfachste Prüffunktionen f(X) und g(Y) die 
Yunktionen 
B 4 


\w(z)dz- ı). 


- & 


yı(X=V3P,(X) =V3 (? 


) 


/ NH 
NM) =V3 PM) =YV3 (? \w,(ydy ) 


verwendet, so entsteht «die Höffdingsche Größe .o 


+1 
\ y9,(Ky, (Y)Ww X, YıdXd/Y 


Lo x 
\ | \w.(r)dr 


3C 4 


\w,(v)dy -w(ze.y)drdy. 


“_e 


Die Vorzüge und Mängel von o als Korrelations 
maß liegen nun klar zutage. Höffdinze kommt es 
besonders auf die Invarianz gegenüber monotonen 
Merkmaltransformationen der Ausgangsverteilung 
an. Diese ist durch die beschriebene Transforma- 
tion von w, (7) und %2(y) auf Gleichverteilung ge- 
währleistet. Dazeren hat o mit dem Korrelations 
koeffizienten r die Mängel xzemein, daß aus r bzw. 
o—0 nicht auf Unabhängiekeit der beiden Merk- 
male geschlossen werden kann — worauf auch Höff 
(dinge ausdrücklich hinweist — und daß im Falle voll- 
ständiger Abhängiekeit o ebensowenig wie r stets 
!en Wert Eins annimmt. 

Während sicher  o 
als auch |o > r sein. Dies hängt davon ab, ob 
für die Ausgangsverteilune w(r, y) oder für die 
wransformierte Verteilung w (X, Y) die Regressions 
iinien besser dem Normalfall der Geragllinigkeit ent 
sprechen. Denken wir uns z. B. eine Ausgangsver 
teilung mit geraden, nahezu zusammenfallenden 
Rkegressionslinien und demgemäß einem Korre 
lationskoeffizienten r nahezu eleich Eins. In diesem 
Fall werden im allgemeinen nach der Transforma- 
tion die Regressionslinien für X und Y weniger 
zut geradlinig sein und daher o|< r. ausfallen. 
Diese Überlegung gibt aber zu starken Bedenken 
regen die Verwendung der Größe o als Korre- 
lationsmaß Anlaß: denn die im angenommenen Fall 
dureh r nachgewiesene große Abhängigkeit besteht 
zu Recht. auch wenn sie bei der Betrachtung von 2 
nieht in Erscheinung tritt. 

Die Ursache dieser Schwäche der Größe » liegt 
in der Sonderstellung einer bestimmten Verteilung, 
nämlich z. B. der Gleiehverteilunge für w, und ws 
innerhalb der Höffdingschen Theorie. So sehr dem 
von Höffding vertretenen Grundsatz, von einem 
idealen Korrelationsmaß Maßstabinvarianz zu for- 
mern, zugestimmt werden kann, ist doch festzu 
stellen, daß o die übrigen, an ein solches zu stellen- 
!en Anforderungen nicht erfüllt. Daß es möglich 
ist, diese Schwierigkeiten zu belieben, zeigt die 
Maximalkorrelation K, die ohne Preisgabe anderer 
Vorzüge auch invariant gegenüber Merkmaltrans- 
formationen ist. 

Veranlaßt vor allem durch die Unmöglichkeit, 
von 2— 0 auf Unabhängigkeit zu schließen, widmet 
Höffding seine Untersuchungen noch anderen mab- 
stabsinvarianten Korrelationsmaßen, nämlich dem 
Pearsonschen Kontingenzmaß und einer interessanten 


K ist, kann soxohll oe < r 


neuartigen Größe, dem von ihm als Abhängigkeits 
maß bezeichneten 22, in unserer Schreibweise. 


+ 


P—=W \ 


w,(z)w.(y)drdy 
o,()w,(y)drdy. 


Man kann nach der Transformation auf X und Y 
zunächst den Integranden nach den Polynomen 
f„(X) und w,ıY) entwickeln, denn nach Gl. (29) 


ist allgemein 


w(z.y) — w, (X) w; (Y) 
nn 
CGjpilr) yily)-w, (2) w;(Y) 
1 
r, und in diesem besonderen Fall 
l 
4 
I. 
1 u ' 


4 GjV2i 


w(X,Y) 


1-V2j-+1P;(X)-Pj(v), 


ml 

wobei fyı o ist. Darauf kann man die vorge 
schriebene Integration mit den Grenzen 1 und X 
bzw. Y durchführen mittels der für die Legendre- 
schen Polynome geltenden Integralformel 


” 


| P.D)42 


I 


an 


die Anwendung der Vollständig 
Darstellung für #2 quadra 
Koeffizienten ci; Anfang 
lautet 


und endlich führt 
keitsrelation zur 
tische Form der 
(dieser Entwicklung 


als 


Der 


+, |] 


Es ist mir nicht gegrlückt, P? etwa ähnlich wie 
das Pearsonsche Kontinzenzmaß mit den Größen 
der Fredholmschen Integralgleichung in Zusammen 
hang zu bringen, und daraus Beziehungen zwischen 
P und den klassischen Korrelationsmaßen odeı 
ier Maximalkorrelation zu gewinnen. Sofort an 
geben läßt sich dagegen wegen cu=0. die Ab- 
schätzung 


> 
od? o? 


8 0,7906 0 


oder &b 


lie auch von Höffding mitgeteilt worden ist. Von 
Interesse wäre vor allem die Frage, ob PK 
sein kann oder sein muß. Ob dies eintreten kann. 
wurde nicht entschieden. Daß dies nicht der Fail 
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sein muß, zeigen Beispiele, die Höflding in seiner 
zweiten Arbeit mitteilt. und bei denen ® kleiner 


als o und beide kleiner als der Korrelationskoefti- 
zient r auseefallen sind. während K selbstverständ- 
lich sogar größer als r ist. 

Zum Schluß sei noch erwähnt. daß bei der An- 
wendung der Höffdingschen Korrelationstheorie auf 
arithmetische Verteilungen grundsätzliche Schwierig- 
keiten auftreten, weil diskontinuierliche  Vertei- 





LUDWIG FISCHER, Die unabzählbare 
Menge, 2. Teil von „Grundlagen der Philosophie 
und der Mathematik“. V-+-57 8. Leipzig 1942, 
Felix Meiner Verlag. Preis brosch. 2,80 M. 

Der Verfasser hat sich zum Ziele gesetzt. die 
üblichen Nachweise für überabzählbare Mengen zu 
entkräften. Nun gibt es bekanntlich philosophische 
Einstellungen zur Mathematik. die Überabzählbarkeit 
ablehnen. Verschiedene mengentheoretische Opera 
tionen und Begriffsbildungen verlieren für diese 
Einstellungen — z. B. den Intuitionismus — in ent- 
scheidenden Punkten ihren Sinn. Der Verfasser 
unternimmt demgerenüber unter anderem den Ver- 
such. direkte Fehler in den Diagonalverfahren durch 
Überlegungen aufzuzeigen. die, soweit sie den An 
spruch mathematischer Argumentation erheben. auf 
Mißverständnissen beruhen. Dabei wird an einieen 
Stellen auf ein von Null verschiedenes Unendlich 
kleines (Abstandswert von .„Urlücken* im Konti- 
nuum u, a.) als „logisches Korrelat des Unendlieh- 
eroben“ zurückgegriffen. 

Arnolü Schmidt. 360 





Potsdam. 


Prof. Dott. DOMENICO ARGENTIERI, Ottiena 
Idustriale XXXI+ 687 S. m. 8 Fig. 4 Taf. 
ü. 3 Tabellen. Mailand 1942—XX, Verlage Ulrieo 
Hoepli. Preis geb. 48 Lire. 

Ein von einem ersichtlich großen Lehrtalent in 
klarer und lebhafter Darstellung zeschriebenes un: 
eut und übersichtlich ausgestattetes Werk von 
nahezu 700 Seiten in Taschenbuchformat. Währenid 
sich die deutschen Bücher meist auf die analytische 
Entwicklung beschränken und den Anfänger in 
einer Unzahl von Formeln hilflos sitzen lassen, gibi 
der Autor die wichtigsten Formeln der theore- 
tischen Dioptrik, der trigonometrischen Durch- 
rechnung und der Wellenoptik und Beugung. 
mit besonderer Berücksichtigung der Korrektur 
und Einstellung für beste Bildschärfe, ohne 
Herleitung an, zeiet aber dafür von Anfang an. 
wie man sie zur Synthese optischer Systeme ver 
wendet. Es wird nieht nur dem Praktiker dureh 
seinen Inhalt, sondern auch dem Theoretiker durch 
die Art der Darstellung gute Dienste tun. Besonders 
hervorgehoben sei. daß es dem Autor durch x 
schickte Substitutionen gelang, die langatmigen 
Seidelschen Bildfehlerausdrücke dünner Linsen für 
numerische Rechnungen zweckmäßig zu verein 
fachen. 

München. 


Staeble. 375 


P. RAETHJEN, Einführung in die Phy 
sik der Atmosphäre, Bd. 1. Statik und 
Thermodynamik. XI +126 Ss. m. 25 Abn. 
Leipzie und Berlin 1942, Verlag B. G. Teubner. 
Preis zeb, 8 M. 

‚Der Verfasser hat es sich zur Aufgabe gemacht. 
die physikalischen Grundlagen der Meteorologie 
elementar, aber exakt und gründlich darzustellen 
und hat dabei ein besonders in didaktischer Hin- 
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iuneen sich nieht in Gleichverteilunz verwandeln 
lassen. Diese Schwierigkeit, mit der sich Höffding 
in seiner zweiten Arbeit auseinandersetzt. besteht 
innerhalb der Eizenwerttheorie der Korrelation 
nieht; vielmehr tritt dort lediglich ein Ersatz der 
Integraleleichungen durch lineare  Gleichungs 
systeme ein. 


EBlineen. Hans Gebelein. 371 





sicht auseezeichnetes Lehrbuch zeschaffen. wie es 
bisher in der meteorologischen Literatur fehlte. 

Die Voraussetzungen an mathematischen und 
physikalischen Kenntnissen, auf denen sich das 
Lehrbuch aufbaut. sind möglichst gering gehalten. 
Es verlangt aber von dem Studenten physikalisches 
Denken und erzieht ihn dazu. Das Buch kann des 
halb nicht leicht «elesen, sondern nur durchgear 
beitet werden. Das ist bei der Fülle des Wissen- 
stoffes, der hier in kürzester Form gebracht wird, 
nicht anders zu erwarten. 

Das Studium des Buches wird dureh die außer 
ordentlich übersichtliche Einteilung des Textes 
sehr erleichtert. Fast in jedem Kapitel kehrt in 
enger logischer Verknüpfung die Reihenfolge 
„Experiment“. „Theorie“, „Satz“, „Übungsaufgabe“ 
wieder. Die Übungsaufzaben vertiefen und beleben 
das Studium und ermöglichen dem Studenten eine 
Selbstbeurteilung. 

Es ist sehr zu wünschen, daß dieses Lehrbuch, 
das den Bedürfnissen der Studenten so zlänzend 
angepaßt ist, weiteste Verbreitung findet. 

Prag. W. Findeisen. 354 


Dr. WOLFGANG RIEZLER, Dozent a. d. Univ. 
Bonn, Einführung in die Kernphysik. 
Meyers Kleine Handbücher, Bd. 6a.) 2. erweit. 
Aufl. 217 S. Leipzig 1942, Verlag Bibliographisches 
Institut A.-G. Preis geb. 2.60 M. 


Dr. WOLFGANG RIEZLER, Dozent a. d. Univ. 
Bonn. Tabellen und Tafeln zur Kern- 
physik. (Meyers Kleine Handbücher, Bd. 6b. 
Ergänzungsband zur Einführung in die Kernphysik. 
15 S. m. 7 Taf. Leipzie 1942, Verlag Bibliogra- 
phisches Institut A.-G. Preis geb, 2.60 M. 

Schon die erste Auflage fiel auf als kurze, klare, 
verständliche, nicht viel besondere Fachkenntnisse 
voraussetzende, dabei zuverlässige Einführung in 
die Kernphysik. Entsprechend den zahlreichen 
inzwischen erreichten Fortschritten ist die zweite 
Auflage ergänzt: der Abschnitt über die Kern- 
reaktionen ist umgearbeitet worden, natürlich ist 
auch die Kernspaltung eingefügt. Die Tabellen 
sind auf den neuesten Stand gebracht und viel 
umfangreicher geworden: besonders die Tabelle 
über die künstlichen radioaktiven Stoffe sei her- 
vorgehoben. Die Haupttabellen sind in dem 
zweiten Bändehen zusammengestellt: der wichtigste 
Inhalt der Tabellen ist in sieben mehrfarbigen 
Tafeln noch einmal dargestellt. In nach Protonen- 
zahl und Neutronenzahl fortschreitenden Schau- 
bildern sind die Isotopie- und Stabilitätsverhält- 
nisse aller bekannten Kerne, die künstlichen Um- 
wandlungen durch «a-Teilchen, Deuteronen. Pro- 
tonen, Neutronen, y-Strahlen und die Kernspaltun- 
een sehr übersichtlich und einprägsam aufgetragen. 
Die. Bändehen können zur Einführung, die Tabellen 
und Tafeln auch zum Nachschlagen für den for- 
schenden Kernphysiker selbst sehr empfohlen 
werden. 

Leipzig. F. Hund. 359 
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NACHRICHTEN 


Arbeitstagung der Gesellschaft für angewandte 
Mathematik und Mechanik in München. 


Die ursprünglich für den 15. und 16. Mai 1942 
geplante Arbeitstagung der GAMM mußte zu 
nächst abgesagt werden. Da sie aber als kriegs 
wiehtig angesehen wurde, fand sie dann doch in 
der Zeit vom 11. bis 13. Juni 1942 in München 
statt. Wie notwendig die Abhaltung einer solchen 
\rbeitstagung war, nachdem seit der letzten Zu 
sammenkunft in Braunschweig über zwei Jahre ver 
tlossen sind, zeigte einmal die große Zahl ler ange 
meldeten Vorträge, die es nötig machte, zweieinhalb 
Tage statt wie ursprünglich beabsichtigt zwei Tage 
für die Tagung anzusetzen. Obwohl eine Reihe 
von Herren, die Vorträge angekündigt hatten, im 
letzten Augenblick am Erscheinen verhindert 
waren. war diese Zeit noch reichlich ausgefüllt. 
Weiter eing das hervor sowohl aus der großen 
Zahl der Teilnehmer der von der Technischen 
Hochschule München zur Verfügung gestellte Hör 
saal mit 180 Sitzplätzen war bei vielen Vorträgen 
fast gefüllt wie aus den an die meisten Vor 
träge anschließenden lebhaften Aussprachen, aus 
denen das große Interesse an den behandelten 
Problemen hervorging. Nicht zu vergessen ist, dab 
einer der Hauptwerte einer solchen Tagung in der 
persönlichen Aussprache der Teilnehmer und in 
den wissenschaftlichen Erörterungen im kleineren 
Kreise liegt. 

Leider war Professor Prandtl, der 
regelmäßig die Tagungen der GAMM leitete, an 
der Teilnahme durch anderweitige starke Inan 
spruchnahme verhindert. An seiner Stelle eröffnete 
Professor Grammel die Tagung und dankte vor 
allem dem Geschäftsführer, Professor Weber, füı 
die dieses Mal besonders viel Arbeit erfordernde 
Vorbereitung der Tagung. 


SONnSst 


Die Vortragsfolge war die folgende: 
Donnerstag. den 11. Juni, vormittags: 

Margeuerre, DVL: Stabilität der 
schale veränderlicher Krümmung. 

Flügge, DVL: Ausbreitung von Biegungswellen 
in Stäben. } 

Weber, Dresden: Eingrenzung der Verschiebun- 
gen belasteter Membranen 

jilharz. Braunschweig: 
eindimensionale Bewegung. 
Donnerstag, den 11. Juni, nachmittags: 

Fromm, Danzig: Kräfte und Molekülzertrümme- 
rung in der Schwingmühle. 

Hartung, Danzig: Messungen des Schubmoduls 
zäher Flüssigkeiten. 

Pfriem. DVL: Technische Beispiele für vorteil- 
hafte Verwendung von Massenpunktkoordinaten. 

Mönch. München: Spannungsoptische Vorfüh- 
rungen. 

Freitag, den 12. Juni, vormittags: 

Collatz, Karlsruhe: Fehlerabschätzung beim 
Iterationsverfahren zur Auflösung linearer 
Gleichungssysteme. 

QWuade, Karlsruhe: Neuere Ergebnisse über tri- 
zonometrische Interpolation. 

v. Guerard, Darmstadt: Über Existenz und 
rechnung größter Biegeschwingungszahlen. 

L,öbell, München: Über ein Verzerrungesmaß für 
Flächenabbildungen. 

Neuber, Braunschweig: Ein 
numerischen Integration von 
chungen höherer Ordnung. 


Zylinder 


Über eine wzesteuerte 


Verfahren zur 
Differentialglei 


Richter. Aussige: Einfluß des 
schnittes auf die 
Wasserschloß. 
Freitag, den 12. Juni. nachmittags: 

Sehultz-G6runow, Aachen: Nichtstationäre 
eindimensionale Gasbewegungen. 

Sauer. Aachen: UCharakteristikenverfahren für 
die eindimensionale instationäre Gasströmune. 
Oswatitsch und Rothstein, Göttingen: Das 

Strömungsfeld in Lavaldüsen. 

Schubert. Peenemünde: Zur Theorie des statio 
nären Verdiehtunesstoßes. 

Schmieden, DVL: Kompressible Strömungen 
um ein Rotationsellipsoid (Rayleigh-Strömungen). 
Sonnabend. den 13. Juni. vormittags: 

l,yra, Göttingen: Zur Theorie des 

Aufwindes. 

Schmidt u, Schröder. 
sche Grenzschichtgleichung als asymptotische 
Näherunz der Navier-Stokesschen Differential 
gleichung bei unbegrenzt wachsender Reynolds- 
scher Zahl. 

Weissineer,. DVL: 
Traglinientheorie. 

Ginzel. Göttingen: 
eines tiefen verwundenen Rechteceksflügels. 

Schubert. DVL: Über eine Integrodifferential- 
sleichung der Luftschraubentheorie. 

Die Vorträge der Herren Flügge, Bilharz, Col- 
latz, Neuber und Schubert sowie der von Frl. 
Ginzel werden voraussichtlich in den nächsten 
Heften dieser Zeitschrift erscheinen. Der Vortrar 
von Herrn Weber beschäftigte sich mit den Fragen. 
denen sein auf den Seiten 126-130 Heftes 
vebrachter Aufsatz gewidmet ist. 

Dresden. 


Kammerquer 
Spiezelbewegung in einem 


Moazagot! 


DVL: Die Prandtl 


Fine Verschärfune der 


Die  Auftriebsverteilung 


dieses 


Willers. 386 
Persönliches. 


Seewald wurde zum o. Prof. für 
Mathematik und >Strömungslehre an 
Hochschule Aachen ernannt. 


Dr.-Ing. F. 
anzewandte 
der Teehnischen 


Berichtigungen. 

W. Frössel, 
Tragkraft eines 

ebener Gleitbahn. 7. 
1941) S. 321 bis 340. 

In der Kontinuitätsgl. (2) muß es natürlich heißen 


Berechnung der Reibung und 
endlich breiten Gleitschuhes auf 
angew. Math. Mech. Bd. 21 


ou dv Ow 


de  dy 02 


Göttingen. W. Frössel. 374 


K.Wieghardt, Zur Theorie der Wirbelbewegung. 
7. anzew. Math. Mech. Bd. 22 (1942), S. 58 bis 60. 
Die Gleichungen (6), die aus (5) folgen, sind 
versehentlich falsch geschrieben; sie lauten richtig: 


en 


1 y ‚(n.Z)(n — MY) 
KR 


- WW +(E— 2 
(s) 


1 [ vn E)(E 
2a) (y — y?+(£ 
(8) 


Aus diesen gehen die speziellen Formeln (7) hervor. 


Göttingen. K. Wieghardt. 38 
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